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XXII. BAND VIERTES HEFT 1954 


Uber Anwendungen von ,,Planen relativer Normalbeschleunigungen“ 
in der Getriebedynamik. 


(Trigheitspol einer ebenbewegten Scheibe und Kriimmungsmittelpunkt von Koppelkurven. ) 
Von O. Tolle. 


1. Einleitung. Infolge der vielseitigen Anwendung in der Getriebedynamik ist es wiin- 
schenswert, zur Ermittlung des Tragheitspoles einer eben bewegten Scheibe mig- 
lichst einfache Konstruktionen zur Verfiigung zu haben. 

Es ist Zweck dieser Arbeit, fiir jeden in Frage kommenden Falleine 
graphische Lésung anzugeben, und zwar liegt in den folgenden Konstruk- 
tionen der Gedanke zu Grunde, fiir zwei Beschleunigungszustinde die Wirkungslinien der 
dadurch bedingten Tragheitskrafte der Scheibe zu finden, deren Schnittpunkt dann der ge- 
suchte Tragheitspol ist. Im Wesentlichen wird dabei von den ,,Planen relativer Normal- 
beschleunigungen‘‘ Gebrauch gemacht, wie sie bereits von M. Tolle! in seiner Arbeit ,,Zur 
Dynamik der eben bewegten Scheibe“‘ entwickelt wurden, um dort fiir die Lésung der dyna- 
mischen Grundaufgabe benutzt zu werden, die darin besteht, Beschleunigungen und Stiitz- 
krafte bei gegebenen eingepragten Kraften zu bestimmen. Eine weitere Anwendung der 
Plane relativer Normalbeschleunigungen“‘ machte spater (im Jahre 1938) N. Rosenauer? 
bei der Bestimmung von Beschleunigungen kinematischer Ketten. Rosenauer gibt in seiner 
Arbeit den Hinweis, daB solche Plane bereits 1909 bei N. Joukovsky zu finden seien, daB 
»aber diese Idee im Schrifttum nicht geniigend ausgeniitzt wurde‘‘. Demgegeniiber sei fest- 
gestellt, da bereits im Jahre 1932 von M. Tolle nicht nur von den,,Planen relativer Normal- 
beschleunigungen“‘ Gebrauch gemacht wird, sondern auch in der angegebenen Quelle hierfiir 
einfache Ableitungen und Beweise mitgeteilt werden. 

Die bekanntesten Verfahren zur Ermittlung des Tragheitspoles sind die von R. Beyer, 
R. Kraus und G. Gerber. 

Bei Beyer® erfolgt die zeichnerische Bestimmung mit Hilfe des Wendepoles und des Ge- 
schwindigkeitspoles. 

Kraus* findet den Tragheitspol aus den Geschwindigkeiten nach Annahme zweier Be- 
schleunigungszustande, indem er die Wirkungslinien der zugehérigen Tragheitskrafte nach 
den im Ing.-Arch. 1930 von O. Tolle® mitgeteilten Konstruktionen bestimmt. Es wird dabei 
fiir den ersten gewahliten Beschleunigungszustand die Tangentialbeschleunigung eines Punktes 
mit fester Drehachse null, im zweiten Fall eine beliebige Beschleunigung des Schwerpunktes 
der Scheibe angenommen. Die von Kraus vorgeschlagene Liésung der Aufgabe vereinfacht 
sich, wenn die vom Verfasser® im Ing.-Archiv 1951 angegebene neue Konstruktion fiir die 
Wirkungslinie des resultierenden Massenwiderstandes angewendet wird. 

Das Verfahren von Gerber? benutzt einen Geschwindigkeitsplan sowie einen beliebigen, 
mit der Bewegung der Scheibe vertraglichen Beschleunigungsplan nach Mehmke, bei welchen 
die Geschwindigkeiten und die Beschleunigungen von einem Punkt aus aufgetragen werden. 

Zum SchluB dieser Arbeit wird, ebenfalls unter Benutzung eines ,,Planes relativer Normal- 
beschleunigungen“, noch eine Konstruktion gezeigt, wie zeichnerisch der Kr iimmungs - 
mittelpunkt einer Koppelkurve bestimmt werden kann. 


1 M. Tolle, ZVDI 76 (1932), S. 799. 

2 WN. Rosenauer, Maschinenbau RM-Afg 6 (1938), S. 543. 

3 R. Beyer, Maschinenbau RM-Afg 3 (1935), S. 642. 

4 R. Kraus, Maschinenbau RM-Atg 4 (1936), S. 337. 

5 0. Tolle, Ing.-Arch. 1 (1930), 8.377; ZVDI 74 (1930), 5S. 1453. 
6 0. Tolle, Ing.-Arch. 19 (1951), S. 355. 

? G. Gerber, Maschinenbau RM-Afg 8 (1940), S. 533. 
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2. Kurze Zusammenfassung der wichtigsten Beziehungen fiir die ,,Pline relativer Normal- | 
beschleunigungen“. In der vorgenannten Quelle trennt M. Tolle die Massenkrifte, welche | 
von der Winkelgeschwindigkeit @ herrithren und durch den Geschwindigkeitszustand der 
Scheibe bedingt sind, von denen, welche infolge der Winkelbeschleunigung ¢ auftreten. 

Ist nach Abb. 1 die Winkelgeschwindigkeit w = 0, so geht bekanntlich die auf dem Pol- 
strahl $8S senkrecht stehende Wirkungslinie g, der resultierenden Tragheitskraft Peo=0) =— bse 
durch den auf den Geschwindigkeitspol ‘$ bezogenen Trigheitsmittelpunkt Ty; dieser liegt 
auf dem Polstrahl $$Sim Abstand ty= ig/@s von S aus. Hierin bedeutet ts = YJs|m den. 
Tragheitshalbmesser der Scheibe, Js das Tragheitsmoment dieser, bezogen auf die in S 
senkrecht stehende Achse, m= G/g die Scheibenmasse und Qs den Polabstand des Schwer- 
punktes. ' 


Abb. 1. Bestimmung der Wirkungslinie des Abb. 2. Ermittlung der Wirkungslinic der resultierenden Tragheitskraft 
resultierenden Massenwiderstandes fiir infolge der Winkelgeschwindigkeit w bei Winkelbeschleunigung ¢ = 0. 
Winkelbeschleunigung ¢ und Winkel- : 

geschwindigkeit w = 0. Abb, 2a. Beschleunigungsplan zu Abb. 2. 


Fir ¢=0 (Abb. 2) geht die von der Winkelgeschwindigkeit @ abhangige Tragheits- | 
kraft Poco) = — m bs» einerseits durch den Schwerpunkt S, andererseits durch den Wende- 
pol W der Scheibe. Die Wirkungslinie von P,-o) liegt in der Geraden g,, welche sich auch 
ohne Kenntnis des Wendepoles mit Hilfe eines ,,Planes relativer Normalbeschleunigungen“ 
folgenderma8en bestimmen ]aBt: 

Man geht (Abb. 2a) von einem dem Dreieck I 2 S ahnlichen Dreieck (1) (2) (S) aus, welches 
dem ersten parallel gelegen ist und bei welchem die Seite (2) (1) gleich der relativen Normal- 
beschleunigung n,, des Punktes 2 gegen I gemacht wird. Dann tragt man in den Anfangs- 
und Endpunkten (2) und (1) der relativen Normalbeschleunigung n,, = v$,/21, wobei 
Vy1 = V. — V, die Relativgeschwindigkeit des Punktes 2 gegen ] ist, die Normalbeschleunigungen 
n, = v2/0, und n, = v?/0, auf; die Senkrechten auf n, und n, in deren Endpunkten (2)’ und 
(1)’ schneiden sich im Punkt (W), in welchem die von w herrithrenden Beschleunigungen by », 
bi}, und bs, endigen. Die Parallele zu bs. durch S in Abb. 2 ist die Wirkungslinie g, von 
P, e=0)° 

Bei gleichzeitiger Beriicksichtigung der Winkelgeschwindigkeit @ und der Winkelbe- 
schleunigung ¢ (vgl. Abb. 3) ist die resultierende Tragheitskraft P—— mbs, deren Wir- 
kungslinie nach einem der bekannten Verfahren! bestimmt werden kann, wenn neben der 
Beschleunigung eines beliebigen Punktes noch die des Schwerpunktes S bekannt ist. 

Nimmt man an, daf die von ¢ herrithrende Tangentialbeschleunigung (oh bekannt ist, 
dann findet man nach Abb. 3a wiederum in einem ,,Plan relativer Normalbeschleunigungen“ 
im Punkt [1]’ den Endpunkt der Gesamtbeschleunigung des Punktes 1 entsprechend der 


* H. Alt, Z. angew. Math. Mech. 6 (1926), S. 58. — K. Federhofer, ZVDI 74 (1930), S. 234. — O. Tolle, 
Ing.-Arch. 1 (1930), 5.377; ZVDI 74 (1930), S. 1453; Ing.-Arch. 19 (1951), S. 355. 
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Vektorgleichung b, = 7,+t,,.+1t,,. Hieraus ergibt sich fiir jeden weiteren Punkt der 
Scheibe, z. B. fir Punkt 2, der Endpunkt [2]’ seiner Beschleunigung 6, auf der Senkrechten 
zu 12 mit (2) als Anfangspunkt. Der Endpunkt [S]’ der Schwerpunktsbeschleunigung bg ist 
dann der Eckpunkt des ahnlichen, um 90° gedrehten Dreiecks 1 2 3; ihr Anfangspunkt ist 
(S). Man findet [S]’ als Schnitt der Senk- 
rechten auf 1 S durch [1]’ mit der Senkrechten 
auf 2S durch [2]’. 


[sl 


L715 


1 Te Ge ' 
(1) ti —*(W) te (1 
Abb. 3. Beschleunigungszustand fiir gleichzeitig vorhandene 
Winkelgeschwindigkeit m und Winkelbeschleunigung e. Abb. 3a. Beschleunigungsplan zu Abb. 3. 


Nachdem, wie beschrieben, fiir Punkt 1 und fiir den Schwerpunkt S die Beschleunigungen 
bestimmt worden sind, findet man nach Abb. 4 hieraus in der Geraden g, die Wirkungs- 


linie der resultierenden Tragheitskraft P, = — m bs nach der Konstruktion des Verfassers!, 
indem man im Lageplan durch den Endpunkt [J] von 6, die Parallele zu 1S zieht und diese 


mit der Verbindungslinie des Endpunktes [S] von bs 
mit dem auf Punkt 1 bezogenen Tragheitsmittel- 
punktes T, in A’ zum Schnitt bringt. Die Parallele 


zu bs durch A’ ist die gesuchte Wirkungslinie von P,. 


vy) 
aS 


es 


Abb. 4a. Nachweis der Konstruktion nach Abb. 4. 


Abb. 4, Zeichnerische Bestimmung der Wirkungslinie 


der resultierenden Tragheitskraft P = — mbg, wenn : : 5 
zuvor bei bekannter Beschleunigung b, eines Punktes 1 Wiahrend die beschriebene Konstruktion in der 


mit Hilfe eines ,,Planes relativer Normalbeschleuni- . ; 7S 
gungen‘‘ nach Abb. 3a die Schwerpunktsbeschleunigung * angegebenen Quelle aus einem »Plan relativer Normal 


by gefunden worden ist. beschleunigungen“ hergeleitet ist, sei an dieser Stelle 
ein anderer einfacher Beweis angegeben (Abb. 4a). 


In der in FuBnote 1 von S. 227 bezeichneten Arbeit wurde gezeigt, daB die resultierende 
Massenkraft P durch einen auf der Verbindungslinie 1 S liegenden Punkt A hindurchgeht, 
dessen Entfernung von S aus gemessen 


t,tsl 
See a 


“bssina 


1 Q, Tolle, Ing.-Arch, 19 (1951), S. 355. 
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betrigt; hierin bedeutet t, = i3/1S die Entfernung des auf 1 bezogenen Tragheitsmittel-_ 


punktes T, und tg; die relative Tangentialbeschleunigung des Punktes S gegen I. 
Aus der Ahnlichkeit der Dreiecke S’ A’ [S] und S T, [S] folgt 


a tsy 
t,  bgsina’ 


d.h. die Richtigkeit der graphischen Lésung. 


3. Konstruktionen fiir dieWirkungslinie g, der resultierenden Massenkraft P,(.,—0), herriihrend | 
von der Winkelbeschleunigung ¢ bei Winkelgeschwindigkeit @ == 0. Wohl in den meisten 
Fallen werden von vornherein die Geschwindigkeiten v bzw. die gedrehten Geschwindig- 
keiten v fiir alle Punkte der Scheibe bekannt sein. Es ergeben sich dann folgende Mdéglich- 
keiten: 

a) Der Geschwindigkeitspol $ liegt zuganglich auf dem Zeichenblatt 
(Abb. 5). Man findet dann den auf S$ bezogenen Tragheitsmittelpunkt Ty mit Hilfe des 
Hoéhensatzes im rechtwinkligen Dreieck, indem man auf der Senkrechten zum Polstrahl $3 S| 
in S die Strecke SC = ig auftragt und in C auf der Verbindungslinie $C die Normale er- 
richtet; diese schneidet auf der Verlangerung von ‘S$ S den gesuchten Punkt Ty ab. 


Abb. 5. Konstruktion des auf den Geschwindigkeitspol ‘§ Abb. 6. Dieselbe Aufgabe wie in Abb. 5, jedoch 
bezogenen Tragheitsmittelpunktes Ty, mit Hilfe gedrehter liegt der Geschwindigkeitspol ‘$ auBerhalb des | 
Geschwindigkeiten bei zuginglichem Pol; die Wirkungs- Zeichenblattes. | 


linie von Pe(w =0 geht durch Ty; und steht auf dem 
Polstrahl 3S senkrecht. 


b) Bei nichtzuganglichem Geschwindigkeitspol (Abb. 6) ist gleich- 
falls die Anwendung des Héhensatzes im rechtwinkligen Dreieck méglich, wenn man fiir 
den Endpunkt C der Strecke SC = igs in bekannter Weise (durch Parallelenziehen) die ge- 
drehte Geschwindigkeit v¢ bestimmt; die Senkrechte auf vc in C schneidet auf der Rich- 
tun von vy den Tragheitsmittelpunkt Ty ab!. 


c) Liegt bei nichtzuganglichem Geschwindigkeitspol der Schwerpunkt S auf der geraden 
Verbindungslinie zweier Punkte 1 und 2 (Abb. 7), deren gedrehte Geschwindigkeiten be- 
kannt sind, so sucht man zunachst vg nach M., Tolle?, indem man den Endpunkt von v, mit 
J und mit dem Endpunkt von v; verbindet. Ziehe dann von S parallel zu vj bis zur ersten 
Hilfslinie, dann weiter parallel zu vj bis zur zweiten Hilfslinie und erhadlt dadurch den End- 
punkt- von vg. Man tragt nun wiederum auf der Senkrechten zu vg in S die Strecke SC = is 
auf. Es ergibt sich dann v¢ durch Parallelenziehen aus v' und vz und nunmehr in gleicher 
Weise wie in Abb. 6 der Tragheitsmittelpunkt Ty. 


d) Sind fiir Punkt 1 und Schwerpunkt S der Scheibe die wirklichen Geschwindigkeiten 0, 
und vs bekannt, so ergibt sich folgende einfache Konstruktion (Abb. 8) fiir die Wirkungs- 


* Eine andere Konstruktion des auf $$ bezogenen Tragheitsmittelpunktes Tg findet sich bei M. Tolle, 
ZVDI 76 (1932), S. 800, Abb. 7; dort wird von dem zum Punkt 1 gehérigen, auf der Strecke 1S ge- 
legenen Tragheitsmittelpunkt T, Gebrauch gemacht. . 

_ » M. Tolle, Regelung der Kraftmaschinen, 3. Auf]. (1921), $.19. Die dort angegebene Konstruktion 
ist fir wirkliche Geschwindigkeiten durchgefiihrt; sie gilt auch fiir gedrehte Geschwindigkeiten. 
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linie g, der allein von ¢ abhangigen resultierenden Tragheitskraft. Die Verbindungslinie des 
Endpunktes von vg mit dem auf 1 bezogenen Tragheitsmittelpunkt T, schneidet sich mit 
der durch den Endpunkt von 0, gezogenen Parallelen zu 1S im Punkt A’, durch welchen, 
parallel zu vs, die Gerade g, hindurchgeht. Beweis: Da die Winkelgeschwindigkeit a = 0 
angenommen wird, so sind nur senkrecht zu den Polstrahlen stehende Tangentialbeschleu- 
nigungen 6b, = @-« méglich. Die gegebenen Geschwindigkeiten der Scheibe v = 9@ sind 


Abb. 7. Lésung der gleichen Aufgabe wie in Abb. 6 fir Abb, 8. Gegeben sind die wirklichen Geschwindig- 
den Fall, daB der Schwerpunkt S auf der geraden Ver- Neesten ah diasden Pan Ltenl au ; 
bindungslinie der beiden Punkte 1 und 2 gelegen und der sian On Poe 3 Sagat a argh al ah 


Geschwindigkeitspol J§ auBerhalb des Zeichenblattes liegt. die Wirkungslinie von P, (w = 0). 


ebenso wie die Tangentialbeschleunigungen den Polabstinden @ proportional und haben mit 
diesen gleiche Richtung, so daB die Wirkungslinie g, in der oben angegebenen Weise nach 
der fiir Beschleunigungen in Abb. 4 beschriebenen Konstruktion gefunden werden kann. 
Der Nachweis dafiir, daB die fiir die Beschleunigungen b, und bs giiltigen Konstruktionen 
auch fiir die Geschwindigkeiten v, und vs an- 
gewendet werden diirfen, kann auch folgender- 
maen gefiihrt werden: 


Loe 
@x7 
if 
py 


IS & 


Abb. 9. Lésung wie in Abb. 8 fiir den Fall, daB die gedrehten Abb. 10. Sonderfall: Die gedrehten Geschwindig- 
Geschwindigkeiten oy und 25 bekannt sind. keiten v, und vg der Punkte 1 und S sind parallel. 


Aus ts = ¥, + vs; folgt durch Differentiation nach der Zeit t 


dvs _ dv, a dvs, 
ifn 9 Oe ify” 


also é ; 4 

bs = by Sh bsy ) 
d.h. die Geschwindigkeiten 0,, ¥s und vs; haben mit den Beschleunigungen gleiche Rich- 
tungen und es ergibt sich also aus v, und vg eine resultierende Tragheitskraft 


P,(y<0) = — mbs = — mb, + (— mbs:). 


Die letzte Gleichung entspricht derjenigen, auf Grund deren die bekannten Konstruktionen 
durchgefiihrt werden kénnen. 


| 
| 
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e) Sind anstelle der wirklichen Geschwindigkeiten die gedrehten Geschwindigkeiten 0, | 


und vg gegeben, so ergibt sich mit derselben Begriindung die entsprechende Konstruktion 


nach Abb. 9: 


Errichte in T, die Senkrechte auf 1 S und durch S die Senkrechte auf v,, dann schneiden > 


sich beide Linien im Punkt A’, durch welchen senkrecht auf vg die gesuchte Wirkungslinie 


von Peo) hindurchgeht. ch. , 
f) Fir den Sonderfall (Abb. 10), daB® vj || vg ist, geht die Gerade g, 1 vs unmittelbar 
durch den Schwerpunkt S der Scheibe hindurch. (tg = is/@s = 0, weil Qs = ©). 


\ 4, Anwendungen: Zeich- 
\ 2 nerische Ermittlung des 
Tragheitspoles T einer eben 
bewegten Scheibe. Wie be- 
reits erwahnt, soll derTrag- 
heitspol als Schnittpunkt 
der Wirkungslinien zweier 
. resultierender Massenkraf- 
te bestimmt werden, die 
zwei beliebig angenomme- 


standen entsprechen. 


nahmen und zwar: 

Fall iva)o=0:,De 
Wirkungslinie Pe(w=0) geht 
durch den auf den Ge- 
schwindigkeitspol {§ hbe- 
zogenen ‘Tragheitsmittel- 
punkt Tx, 

b) e=0: Die Wirkungs- 
linie Peo) ist nur vom 
Geschwindigkeitszustand 
abhangig und wird aus 
einem  ,,Plan_ relativer 
Normalbeschleunigungen“ 
nach Abb. 2a _ bestimmt, 
oder 


(1)' Wirkungslinie P,, (e=0) Wird 


Abb. lla. l Abb. 11. 


(W) wie bei Fall 1b) gefunden. 
b) e0: Mit Hilfe eines 


Abb. 11. Bestimmung des Tragheitspoles, Wendepoles und Wendekreises einer eben . 


Fall 2: a) e=0: Die} 


nen  Beschleunigungszu- | 


Die folgenden Liésungen 1 
* beruhen auf zwei von ein- | 
ander getrennten An- 


bewegten Scheibe mit Hilfe eines ,,Planes relativer Normalbeschleunigungen‘ bei ,Planes relativer Normal- | 
J 


zuganglichem Geschwindigkeitspol ‘}. 


Abb. lla. Beschleunigungsplan zu Abb. 11. (Beschleunigungen gegeniiber Abb. 11 on 
~ "in dreifacher GréBe.) zunachst nach Abb. 3a fiir 


eine beliebig angenommene 
Beschleunigung eines Punktes die zugehdrige Schwerpunktsbeschleunigung ermittelt und 


darauf aus diesen beiden Beschleunigungen die Wirkungslinie von P, nach Abb. 4. 


5. Beispiele. a) Gegeben in Abb. 11 die gedrehte Geschwindigkeit v;, des Punktes 1 
sowie die Kriimmungsmittelpunkte K, und K, der Geschwindigkeitspol §§ ist zuginglich. | 


Gesucht der Tragheitspol J, Wendepol W und Wendekreis. Die Lésung lautet: 

Aus 0; folgt in bekannter Weise durch Parallelenziehen v, und vs. Der auf % bezogene 
Tragheitsmittelpunkt Ty wird wie in Abb. 5 gefunden; die Senkrechte g, in Ty auf Pol- 
strah] S$ ist die Wirkungslinie von Pe(w=0)- Aus 0 = 0, — v, findet man n 1 = U5)/ 2a 
(Héhensatz im rechtwinkligen Dreieck), ebenso n, = v?/0, und ny = v3/o,. In Abbie liefer 


beschleunigungen“* wird | 


| 
| 
] 


| 
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der ,,Plan der relativen Normalbeschleunigungen“ hieraus, entsprechend Abb. 2a, den 
Punkt (W) und die von @ fiir ¢= 0 abhangige Schwerpunktsbeschleunigung bs,,. Die 


Parallele g, hierzu durch den Schwerpunkt S ist die Wirkungslinie von De und ihr 
Schnittpunkt mit der erst gefundenen Geraden g der Tragheitspol T der Scheibe. Da alle 
von @ herrithrenden Beschleunigungen nach dem Wendepol W hinweisen, so findet sich 
dieser in Abb. 11 als Schnittpunkt der Parallelen zu 5, durch I mit der Parallelen zu Boss 
durch 2. Die Strecke SW ist der Durchmesser des Wendekreises der Scheibe. 


Abb, 12a. Beschleunigungsplan zu 
Abb. 12. Konstruktion des Tragheitspoles bei nicht auf dem Zeichenblatt Abb. 12. (Beschleunigungen gegeniiber 
gelegenem Geschwindigkeitspol. Abb. 12 in doppelter GréBe.) 


b) Dieses Beispiel (Abb. 12) unterscheidet sich von dem vorhergehenden nur dadurch, 
daB hier der Geschwindigkeitspol nicht zuganglich ist. Die Lésung der Aufgabe ist mit Hilfe 
des ,,Planes der relativen Normalbeschleunigungen“ (Abb. 12a) durchgefihrt, indem einmal 
fiir « = 0 die durch den Schwerpunkt gehende Gerade g, || bs » (wie in Abb. 2 und 2a) und 
dann fiir ¢+0 (entsprechend Abb. 3 und 3a) fiir eine beliebig angenommene Tangential- 
beschleunigung t;, die Beschleunigungen b, und bs ermittelt worden sind. Die Wirkungs- 
P,, welche durch die beliebig gewahlte Winkelbe- 
schleunigung ¢ bedingt ist, wird dann, wie in Abb. 4 angegeben, aus b, und bs mit Hilfe des 
auf 1 bezogenen Tragheitsmittel punktes T, gefunden. In Abb. 12 ist dann noch die Wirkungs- 


linie g, der resultierenden Massenkraft P 


€? 


linie g, der resultierenden Tragheitskraft P,(.—9) eingetragen, welche ebenfalls durch den 
Tragheitspol & gehend, auf vg senkrecht steht und hierauf den zum Geschwindigkeitspol 
gehérigen Tragheitsmittelpunkt Ty abschneidet. 

c) Wird, wie in Abb. 13, die Relativgeschwindigkeit v,, klein, also auch die relative Nor- 
malbeschleunigung 7,, klein, dann kénnen in solchem Fall, wie in Abb. 13a gezeigt, auch die 
Ecken des sehr kleinen, parallel] zum Dreieck 1 2S gelegenen, ahnlichen Dreiecks (1) (2) (S), 
ohne daB dadurch eine wesentliche Ungenauigkeit eintritt, zu Anfangspunkten der von w 
herriihrenden Beschleunigungen gemacht werden, da ja im Allgemeinen die Normalbe- 
schleunigungen n, und 7, verschieden grof sein werden. Der Schnittpunkt (W) der beiden 
Senkrechten in den Endpunkten (1)’ und (2)’ auf n, und n, kommt dann meist auBerhalb 
des Zeichenblattes zu liegen. Die Bestimmung der Richtung von bs, ist, wie Abb. 13a zeigt, 
durch Zeichnen zweier ahnlicher Dreiecke méglich, d.h. es ist die geometrische Aufgabe zu 
lésen: von einem Punkt (S) aus die Verbindungslinie nach einem nicht zuganglichen Schnitt- 
punkt zweier Geraden (Biss und Boa) zu konstruieren. Praktisch wird es oft geniigend genau 
sein, wenn bei kleinen relativen Normalbeschleunigungen n,,, wie in Abb. 13b, die drei 
Punkte (1), (2) und (S) in einem Punkt (0) zusammengelegt werden und dieser zum Anfangs- 
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punkt aller Beschleunigungen gemacht wird. Die Konstruktion des Tragheitspoles z iq | 
Abb. 13 als Schnittpunkt der Geraden g, und g, ist einerseits mit Hilfe der pe na | 
Abb. 9, andererseits mit Hilfe des ,,Planes der relativen Normalbeschleunigungen“ ftir den | 
Beschleunigungszustand ¢ = 0 entsprechend Abb. 2a durchgefithrt. ea = actaail | 
d) Sind im Sonderfall (Abb. 14) die gedrehten Geschwindigkeiten vj und v, parallel, s 
wird die Relativgeschwindigkeit %,,—= 0, daher auch die relative Drapes 
iy, = 0. Tragt man also in Abb. 14a die Normalbeschleunigungen 7, und n,, die beide gleich- | 


Abb. 13. Bestimmung des Tragheitspoles fiir den Fall, daB die Relativ- Abb. tac eondertall: Die gedrehten Geschwindigkeiten 
geschwindigkeit Voy des Punktes 2 gegen 1 klein ist. v, und v, sind parallel, d. h. v,, = 0. 


Abb. 13b. (7/’ (1)! Abbe isa: 


Abb. 13a. Beschleunigungsplan zu Abb. 13. 


Abb. 13b. Beschleunigungsplan zu Abb. 13 mit der zulassigen Annahme, da bei 
kleiner Relativgeschwindigkeit v,, die drei Punkte (1), (2) und (S) in einem Punkt (0) 
zusammengelegt werden, 


Abb. 14a. Beschleunigungsplan zu Abb. 14. 
(Beschleunigungen gegeniiber Abb, 14 in 
doppelter GréBe.) 


gerichtet sind, von (0) aus ab, so erkennt man, daB® der Punkt (W) im Unendlichen liegt, © 
die Beschleunigung bs ,, des Schwerpunktes, weil ebenfalls nach (W) gerichtet, auch senk- 
recht auf n, bzw. n, und somit auch auf v, steht. Die te entgegengesetzt gerichtete re- 
sultierende Tragheitskraft Page) geht durch den Schwerpunkt S und fallt, wie in Abb. 10 


gezeigt, mit der Wirkungslinie der resultierenden Tragheitskraft Boe zusammen. Man 
benotigt daher noch eine weitere Wirkungslinie einer resultierenden Tragheitskraft, die nach 
Abb, 14a aus dem ,,Plan der relativen Normalbeschleunigungen‘ dadurch gefunden werden 
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kann, daB man, wie in Abb. 3a, eine beliebige Tangentialbeschleunigung tts wahlt. Errichtet 
man in deren Endpunkt [1]’ die Senkrechte auf 1 2, so erhalt man auf der Parallelen zu t,, 
durch (2)’ den Endpunkt [2]’ der von ¢ herriihrenden Beschleunigung b,. Die Senkrechten 
durch [1]’ zu I S und durch [2]’ zu 2 Sliefern den Endpunkt [S]’ der Sofapescheaden Schwer- 
punktsbeschleunigung bs mit ihrem Anfangspunkt in (S) = (0). Nachdem auf diese Weise 
zwei Zusammengehiorige, mit der Bewegung der Scheibe vertragliche Beschleunigungen b, 

und bs gefunden sind, ist in Abb. 14 die Gerade gy || bs als Wirkungslinie der entsprechenden 
resultierenden Massenkraft P, wie in Abb. 4 bestimmt worden. Der gesuchte Tragheits- 
pol Z liegt im Schnittpunkt der beiden Geraden g, und gy. 


e) Zum Schlu8 soll noch der Fall behandelt werden, bei welchem der Schwerpunkt S auf 
der geraden Verbindungslinie der beiden Punkte 1 und 2 gelegen ist, deren gedrehte Ge- 
schwindigkeiten v; und v, als bekannt vorausgesetzt werden (Abb. 15). Hierfiir bestimmt 
man, nachdem die gedrehte Geschwindigkeit vj des Schwerpunktes gefunden ist, zunachst 


J 
zor . 
Abb. 15. Der Schwerpunkt S liegt auf der geraden Verbindungs- 
linie zweier Punkte 1 und 2 2 _der eben bewegten Scheibe; 
gegeben y und v, Abb. 15a. Beschleunigungsplan zu Abb. 15. 


die Wirkungslinie g, der resultierenden Tragheitskraft P,.(,, 9), z.B. in Abb.15 nach der 
Konstruktion der Abb. 9 mit Hilfe des auf I bezogenen dine heiteuitvel panktes T, aus %, 
und ts 

Im Evian der relativen Normalbeschleunigungen“ nach Abb. 15a wird, nachdem zuvor 
wiederum ny, = v3,/1 2, n, = v?/0, und n, = v3/0, bestimmt worden sind, die Lage des Punktes 
(S) auf der Geraden (1) (2), etwa nach dem Strahlensatz wie folgt gefunden: Ziehe von (1) 
aus eine beliebige Gerade und trage auf dieser die Strecken (1) H = 1S und hierauf weiter 
H J = S2 auf, verbinde J mit (2), dann schneidet die Parallele zu J (2) durch H auf (1) (2) 
den gesuchten Punkt (S) ab. Man findet dann in der Strecke (S) (W) die von w herriithrende 
Schwerpunktsbeschleunigung bs,, und zu dieser parallel durch S gehend die Wirkungs- 
linie g, der resultierenden Massenkraft, welche fiir ¢ — 0 allein von w herriihrt. Die Geraden g, 
und g, schneiden sich im Tragheitspol X der Scheibe. 


6. Konstruktion des Kriimmungsmittelpunktes von Koppelkurven mit Hilfe eines ,,Planes 
relativer Normalbeschleunigungen“. 

Die meisten bekannten Methoden zur zeichnerischen Bestimmung des Kriimmungsmittel- 
punktes fir die Bahn eines Koppelpunktes setzen voraus, da8 der Geschwindigkeitspol zu- 
ginglich ist, da® ferner z. B. der auf dem Normalstrah] des Koppelpunktes gelegene Wende- 
punkt (Zickzackkonstruktion von Griibler!) oder, zur Anwendung des Bobillierschen Satzes, 
der Wendekreis, der Wendepol oder die Polbahntangente (Aronholdsche Konstruktion”) be- 
kannt ist. Das Verfahren von Hartmann® beruht auf der Verwendung der Polwechselge- 


1 M. Griibler, Getriebelehre (1917), S. 120. 
2 §. Aronhold, Grundziige der kinematischen Geometrie 1872, 
3 W. Hartmann, Maschinengetriebe 1913, 
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schwindigkeit. Von den vorstehenden Bedingungen frei ist das Verfahren von Péschl'; aus- 
gehend von der Gleichung fiir die Normalbeschleunigung n= v2/0, findet man den Kriim- 
mungshalbmesser 0, fiir die Bahn eines beliebigen Punktes 3, wenn neben der Geschwindig- 
keit 0, dessen Normalbeschleunigung n, bekannt ist. Nimmt man an, daf von einer bewegten 
Scheibe die gedrehten Geschwindigkeiten und v, zweier Punkte I und 2, sowie deren Kriim- 
mungsmittelpunkte K, und K, gegeben sind, so wird nach Péschl zunachst in bekannter 
Weise die gedrehte Geschwindigkeit v; des Koppelpunktes 3 bestimmt. Nach Wahl einer 
beliebigen Beschleunigung ‘b, wird dann die des Punktes 2 aufgesucht und aus 6, und 6, nach 


dem Satz von Burmester die Beschleunigung bs des Punktes 3 gefunden. Die Normale zur 
Bahnkurve ist durch die Richtung von v, gegeben; die in diese Richtung fallende Kompo- 


nente 7, von b, liefert den gesuchten Kriimmungshalbmesser Q, und nach Halbkreiskonstruk- 
tion daraus den Kriimmungsmittelpunkt K3. 


3" 


7 
é ké 
W1//1 
Abb. 16. Neue Konstruktion des Kriimmungsmittelpunktes von Abb. 16a. Beschleunigungsplan zu Abb. 16. 


Koppelkurven mit Hilfe eines ,,Planes relativer Normal- 
beschleunigungen“. 


In der nun folgenden neuen Konstruktion wird ebenfalls davon ausgegangen, daB von 
der Scheibe in Abb. 16 die gedrehten Geschwindigkeiten v, und v, der Punkte J und 2 be- 
kannt sind, desgleichen fiir diese die Kriimmungsmittelpunkte K, und K,. Zeichnet man 
aus Ny, = V3,/1 2, ny = vi/o, und n, = v3/0, den ,,Plan der relativen Normalbeschleunigungen“ 
in Abb. 16a, so findet man 63,, d.h. die nur von @ herriihrende Beschleunigung nach Rich- 
tung und GréBe als die Strecke (3) (W). Zerlegt man dann 63, in ihre beiden Komponenten n, 
in Richtung von v; und t3,, senkrecht dazu, so ist m, nach Grée und Richtung die Normal- 
beschleunigung des Koppelpunktes 3. In Abb. 16 ist nun der Kriimmungsmittelpunkt K, 
nach dem Héhensatz im rechtwinkligen Dreieck hieraus folgendermaSen gefunden: Trage im 
Punkt 3 die wirkliche Geschwindigkeit 1, (Endpunkt 3’) und ebenso von 3 aus die im Be- 
schleunigungsplan Abb. 16a gefundene Normalbeschleunigung n, in entgegengesetzter Rich- 
tung auf (Endpunkt 3’’). Die Senkrechte in 3’ auf die Verbindungslinie der Punkte 3’ und 
3’ schneidet auf der Richtung von v, den gesuchten Kriimmungsmittelpunkt K, ab. — Auch 
diese Liésung ist, sofern der Kriimmungsmittelpunkt tiberhaupt auf das Zeichenblatt fallt, 
in allen Fallen méglich, da, wie eingangs gezeigt, ein ,,Plan relativer Normalbeschleunigungen“ 
stets gezeichnet werden kann. 


(Eingegangen am 29, September 1953.) 
Anschrift des Verfassers: Professor Dipl.-Ing. 0. Tolle, Konstanz (Bodensee), Bahnhofplatz 10. 


1 Th. Péschl, Kinfiihrung in die ebene Getriebelehre (1932), S. 71, 
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Uber den geraden Knickstab mit begrenzter Durchbiegung. 
ji Von H. Link 


1, Einleitung und Voraussetzungen. In der vorliegenden Arbeit wird das Verhalten eines 
geraden diinnen Knickstabes untersucht, dessen Durchbiegung nach beiden Seiten durch 
ebene, starre Flachen begrenzt wird. Die zwei zueinander parallelen Begrenzungsflachen 
sind senkrecht zur Verformungsebene des Stabes angeordnet und haben beide von dem zwischen 
ihnen befindlichen Stab den gleichen Abstand f. Das Maf8 f gibt die gréBte mégliche Durch- 
biegung des Stabes an und soll sehr klein gegen dessen Lange sein. 

Der Stab wird aus einem Werkstoff mit linearem Elastizitatsgesetz angenommen und soll 
alle im Verlauf der Untersuchung vorkommenden Verformungen ertragen kénnen. Die 
Biegesteifigkeit EJ ist iitber die Stablange unveranderlich. 

Es werden der an beiden Enden gelenkig, der an beiden Enden eingespannt und der an 
einem Ende gelenkig, am anderen Ende eingespannt gelagerte Stab behandelt. In allen 
Fallen ist eine seitliche Unverschieblichkeit der Stabenden unterstellt. 

Die Grundlage aller durchgefiihrten Betrachtungen ist die Differentialgleichung des ge- 
bogenen, geraden Stabes 


EJy’=—M, (1) 


die den Zusammenhang zwischen der Durchbiegung y des Stabes und dem Biegungsmoment M 
ausdriickt. Der Einflu8 der Querkraft auf die Stabverformung ist nicht beriicksichtigt. 
Ebenfalls sind Reibungskrafte, wo solche etwa auftreten kénnen, vernachlassigt. 


2. Der an beiden Enden gelenkig gelagerte Stab. Mit den Bezeichnungen der Abb. | wird 
das Biegungsmoment im Stab an der Stelle x 
M=Py—Qx. 
Damit geht (1) iiber in 
yt y= Sx mit a= a 
und hat die allgemeine Liésung 


y = Asin ax + Beosax +x. 
l 3 . 
Mit den Grenzbedingungen y = 0 fiir x= 0 und y’ = 0 fiir x = 7 ergeben sich die Kon- 
stanten daraus zu 
as | 
fry Ghana! 
COR —— 


2 


Aus der Forderung, daB sich die Stabmitte um das Ma& f durchbiegen soll, folgt fiir die Stiitz- 
kraft Q als Funktion von a in dimensionsloser Darstellung 


(3) 
Qe ey 


SESS re ad 


57 3 


Te se 


(2) 


Der Verlauf von (2) ist in Abb. 2 dargestellt. Der Wert Q1°/8 EJ f ist zunachst negativ. 
Die Krafte Q miBten daher fiir kleine Werte o 1/2 zur Erzeugung der Durchbiegung f in der 
Stabmitte entgegengesetzt wie in Abb. 1 dargestellt wirken. Fir a 1/2 a 0 ergibt sich der 
aus der Biegelehre bekannte Wert Q//8EJf = — 3, und fiir al/2 =a7/2 wird C=(; Die Zu- 
gehiérige Stabkraft hat dann die Gréfe P= w/l- EJ der Eulerschen Knickkraft eines a 
beiden Enden gelenkig gelagerten Stabes von der Lange 1. Fir « y2=2 wird Q 13/8 EJf = 2°, 
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Die zugehérige Stabkraft ist 
W OEJ=45Es 

und entspricht der Ewlerschen Knickkraft eines an beiden Enden gelenkig gelagerten Stabes 
von der Lange J/2 oder eines an beiden Enden eingespannten Stabes von der Linge I. SchlieB- 
lich geht noch fiir «1/2 > 1,43 x die Stiitzkraft Q—> co, wobei die Stabkraft die GréBe 

m2 

P = 2,045 aa EJ 

hat. Diese Stabkraft entspricht der Eulerlast eines an einem Ende gelenkig, am anderen Ende 
eingespannt gelagerten Stabes von der Lange 1/2. 


es 
BEI 


Abb. 1, Der an beiden Enden 
gelenkig gelagerte Stab. 


Stabform fiir + < a Ze Abb. 2. 


al 
in Abhiangigkeit von ~~ ; 


Qu 
8EJf 
Die Durchbiegung des Stabes wird unter Beriicksichtigung der Grenzbedingungen und 


der Gleichung (2) 
— y sin a x 
J coll all ( al “4 rae (3) 


85 5} es 


Der durch Abb. 1 charakterisierte Belastungs- und Verformungszustand des Stabes kann 
bei wachsendem P nur solange bestehen, bis « 1/2 = a geworden ist. Aus der zweiten Ab- 
ie der Gleichung (3) folgt namlich bei Beschrinkung auf die nachste Umgebung von 
ali2=—axa 


cee fiir a 
a2 tg — al 
(Y )eny2 = en = 0 fiir rem 
ao 
<0 fiir ae me 


Kine positive Krimmung der Stabmitte, wie sie sich fiir a 1/2 > a ergibt, ist wegen der Art 
der Durchbiegungsbegrenzung nicht mdglich, da dazu der verformte Stab die in Abb. 3 dar- 
gestellte Form haben miiBte. Die Gleichungen (2) und (3) sind daher nur solange giiltig, 
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bis die Stabkriimmung an der Anlagestelle 
grenzenden Flache wird. 


Bei weiterer Steigerung von P(q« 1/2 >z) teilt sich die Stiitzkraft 2 Q in zwei Einzelkrafte Q, 
die symmetrisch zur Stabmitte angreifen, wie es in Abb. 4 dargestellt ist. Mit den dort ver- 
wendeten Bezeichnungen werden die Biegungsmomente im Stab 


Mi Fy, 02, , M, = Py,—Q%, 


gleich der Kriimmung der die Durchbiegung he- 


und die allgemeinen Lésungen von (1) lauten, wenn fiir x, und x, einfach x geschrieben wird, 


y= Aysin ax + B, cos ant 2 x 
mit a? = + 7 e 
: EJ 
Yo = A, sin ax + B, cos ax + % x | 


Daraus ergibt sich mit den Grenzbedingungen 
y,=— 0 fir +=0, 
=f fir «= ras 
Voo= jf fur 2 &, , 


/ 


v1 0 fir x == 4, 


¥2 = 0 fir x= x,, 
vo 0 fir x= 1/2 


B, = 0 und das Gleichungssystem 
i A, sin a% + & xy ; 


f = Ay sin 0x) + B, cos 0x) + ne 


e Q * 
0 = A, a cos ax + >> (*) 
0 = A, cos ax, — B, sin ax, 
Abb. 3. Nicht még- Abb. 4. Mégliche Stab- 
al eel! licheStabverformung verformung fiir 
0 = A, cos — — B, sin a wl 
2 fiir 2 >%. ye as 


Die beiden letzten Gleichungen dieses Systems kénnen nur befriedigt werden entweder durch 
a) Ao B, ==4) 
al 


oder durch b) [HUTT Om (We seal 25S, cers) 


Die Bedeutung dieser beiden Bedingungen wird nachfolgend untersucht. 


a) A,= B,=0. Mit dieser Bedingung ergeben sich aus dem Gleichungssystem (*) die Be- 
ziehungen 


apy) Si 
foo x, A, sinax% = 0, A, = be en (**) 
Wegen A,+0 folgt aus der mittleren dieser Gleichungen sinax)=0 oder ax=na 
(n = 1,2, 3,....). Mit n = 1 wird dann aus der ersten der drei Gleichungen (=) 
(2) 
Q 18 = 2 (4) 
SUBS, os 


Der Verlauf von (4) ist ebenfalls in Abb. 2 aufgetragen. 
Fiir die elastische Linie dieses Belastungszustandes gelten dann die Gleichungen 


y, = (sin ax + a2), =f. 
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Der verformte Stab und die daran angreifenden Krafte sind in Abb. 5 dargestellt. Die 
Resultierende aus P und Q am Stabende geht dabei stets durch den Stabquerschnitt im Ab- 
stand x) vom Stabende hindurch. 

In Abb. 6 ist das Verhaltnis x,):1/2 in Abhangigkeit von a1/2 nach der fiir n = 1 geltenden 
Gleichung 

Xo 7% 
2 a W/2 
dargestellt. 


b) a =nz+ ax, Mit n= 1 wird aus dieser Bedingung 


Xo TG 


1/2 Sho 


Der Verlauf dieser Funktion ist auch in Abb. 6 eingetragen. Mit Hilfe 
von a1/2 = 2% + ax, ergibt sich aus dem Gleichungssystem (*) 


ee: ; 


SHIT fal 
2 


Der Verlauf von (6) kann der Abb. 2 entnommen werden. 


al 
It 85 


Zo 


Abb. 5. Stabform fiir 
l 
ws ~ < 2% nach 
der Bedingung a). Abb. 6, x,:1/2 in Abhangigkeit von « 1/2 nach den Bedingungen a) und b). 


a 
VA 


3 


Fiir die elastische Linie des verformten Stabes gelten die beiden Gleichungen 


— i sin a x 
tal ar ol een 
m— tg cos 


yy, i 2; (7) 
al i (Jal l ol | 
Sf peng Ta tg 5 cos [a( x) | 5 sy : 
gen ean ep 


Die Stabform und die am Stab angreifenden Krafte zeigen Abb. 7a fiir 27> a1/2 > 32/2 
und Abb. 7b fiir «1/2 = 37/2. 

Der sich aus den bisherigen Ableitungen ergebende Verformungsvorgang des Stahes stellt 
sich in Abhangigkeit von der Stablast folgendermaBen dar: 

Ab a«l/2 = 2/2 wird mit wachsendem P die Kriimmung der an der ebenen Stiitzflache an- 
liegenden Stabmitte (Abb. 1) immer kleiner und verschwindet schlieBlich fiir «1/2 =m. Bei 
weiterer Steigerung von P bildet sich ein an der Stiitzflache anliegender, zur Stabmitte sym- 
metrischer, ungekriimmter und damit momentenfreier Teil des Stabes aus, an dessen beiden 
Enden die Einzelkrafte Q angreifen (Abb. 5). Die Lange dieses anliegenden Stabteiles wachst 
mit der Stablast, wodurch sich die Stiitzkrafte Q zu den Stabenden hin verschieben. 

Fir «1/2 = 2 7 erreicht der gerade Teil des Stabes die Linge 1/2. Das zugehérige P ent- 
spricht der Eulerlast eines an beiden Enden eingespannten Stabes von der Lange 1/2, so daB 
jetzt der gerade Teil des Stabes ausknicken kann. Dabei entfernen sich die Krafte Q weiter 
von der Stabmitte, so daB wahrend des Ausknickens die Knicklange wachst. Der Knickvor- 
gang kann daher mit kleiner werdender Stablast P aufrecht erhalten werden. Es tritt 
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Durchschlagen der Stabmitte ein, wenn die Last P die der Beziehung al/2 = 2 7 entsprechende 
GroBe beibehalt. Hierbei geht der Stab von der in Abb. 5 dargestellten Form sofort in die Form 
der Abb. 7b iiber. Dieser Vorgang erfolgt schlagartig und sein Ablauf wird anschaulich klar 
bei Betrachtung des mit (6) bezeichneten Kurvenzuges der Abb. 2. Hiernach nimmt beim 
Ubergang von der in Abb.5 dargestellten Stabform zu der in Abb. 7b gezeigten der Wert od/2 ab. 

Wird der nach Abb. 7b verformte und mit a 1/2 = 3 7/2 belastete Stab (Q = 0) als neue 
Anfangsform aufgefa8t, weil sie in ihren Teilen der Stabform nach Abb. 1 bei einer Belastung 
al/2= 7/2 mit ebenfalls Q = 0 véllig entspricht, so gelten bei Steigerung von P die gleichen Uber- 
legungen fiir jedes Stabdrittel, wie sie fiir den ganzen Stab durchgefiihrt worden sind. Alle 


Abb. 7a. Stabform fiir Abb. 7b. Stabform fiir Abb. 8. Eine weitere, Abb. 9. Der an beiden 
al 32 al 32 mégliche Stabform Enden eingespannt ge- 
Ee ee nach = oT nach fi al lagerte Stab. Stabform 
4 ‘ ir —= 27. l 
der Bedingung b). der Bedingung b). 2 fiir = < “ Za, 


abgeleiteten Formeln behalten daher Giiltigkeit, wenn in ihnen unter Beibehaltung der Be- 
deutung von Q als Querkraft an den Stabenden die Stablange J durch 1/3 ersetzt wird. Kurven- 
zug (2’’) in Abb. 2 veranschaulicht den Zusammenhang zwischen Q und P fiir diesen Fall im 
Bereich 37/2 <al/2 <5a2/2. Bei dem vorerwahnten Durchschlagen der Stabmitte springt 
bei der Belastung «1/2 = 2” der Wert Q/3/8 EJf von Punkt A der Kurve (4) nach Punkt B 
der Kurve (2”’). 

Noch ein Umstand soll hier erwahnt werden. Fiir al/2 = 7 ist nicht nur die hier erwahnte 
Stabform, sondern auch die in Abb. 8 dargestellte méglich. Mit al/2 = a (Q = 0) als Anfangs- 
belastung fiir diese Stabform gelten dann ebenfalls die bisher abgeleiteten Beziehungen, 
wenn in ihnen / durch 1/2 ersetzt wird. Die Kurvenziige (2’) und (4’) der Abb. 2 geben den zu- 
gehérigen Zusammenhang zwischen Q und P fiir diesen Fall an. 

3. Der an beiden Enden eingespannt gelagerte Stab. Die Untersuchung verlauft im ersten 
Teil genau wie beim Stab mit beiderseits gelenkig gelagerten Enden. Nach Uberschreitung der 
Eulerschen Knicklast stiitzt sich der Stab gegen eine der Flachen ab, wie es in Abb. 9 dargestellt 
ist. Mit den dort verwendeten Bezeichnungen ergibt sich mit Hilfe der Grenzbedingungen aus 


der allgemeinen Lisung von (1) 
L Ne 
QL 4 


(8) 


OE fe ihe ite 
reer a 
aly’ al 
GP ep es 
RE Ef |) Tw als 
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Fir die elastische Linie des verformten Stabes gilt die Gleichung 


y= —— fa [iin ae + tw "G 0 — cos) — a]. (9) 
2(eE —4) 
Die Stabkriimmung an der Anlagestelle verschwindet fiir ol/ 
geschriebenen Formeln Giltigkeit. 
Fir Werte o 1/4 >a ergeben sich mi 
ander unabhangigen Lisungen 
a) A, = B,= 0, 
b) a ne oe 


4— 2. Bis dahin haben die an- 


t den Bezeichnungen von Abb. 10 die beiden vonein- 


140 


| 
100 + { ; ze 
‘ii fa) i y, 
60 eae | ce | 
| I} y 
fa WAN 
60 = IL 
i 
L | 
40 whe : 
20 —4 | 
acl 
4 | 7 
ee / 
Abb. 10. Mégliche Wt 
Stabverformung fiir = 2 = St z = = £2 
= S kp Abb. 11. fa in Abhingigkeit von x : 


Mit der Bedingung a) folgt aus dem hier nicht angeschriebenen, mit Hilfe der Grenz- 
bedingungen aus den allgemeinen Lisungen von (1) gewonnenen Gleichungssystem ax)=22an 
(n = 1, 2,3---) und damit, wenn fiir n = 1 gesetzt wird, : 


fa 

OB Lue 

So: fae Sas (10) 
M, = 9 


n=<@x— sin «x), 


(11) 
vie 
Die Bedingung b) liefert mit n = 1 auf die gleiche Weise 


Qi 4 
Ms 


al al 
et ky Game Os 


(12) 
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ge (J cue 3! 
a ae 
Zz 1 
Bf; 1 2(ue +3 . 
ts > ( +) od t é cos a5 — teG -x)— al (13) 
ot! 
Sierra mea ‘| 


In Abb. 11 ist der Verlauf der Funktionen (8), (10) und (12), und in Abb. 12 der Verlauf von 


entsprechend der Bedingung a) und 1 =)]— am entsprechend der Bedingung b) 


a 
2 a U4 
aufgetragen. Zp 

Die Abb. 13a, 13b und 13c 
zeigen die Stabformen nach 
(13) fiir die Belastungen a 1/4 08 


a, E 
= 32/2, «l/4—a und al/4 <a. ag ala 

Die Belastung Pdesin Abb. % Hie ihe 4 
13a dargestellten Stabes ent- yl ee a) 

spricht gerade der Eulerlast ~ 

eines an beiden Enden ein- UA ip = ee ea a a ea a a ; 
gespannten Stabes von der F Le 
Lange 1/3, so daB der gerade a Ue on a 3 iE ba & 


Teil des Stabes ausknicken 
kann. Der Knickvorgang kann 
wegen der zunehmenden 
Linge des mittleren Stabteiles mit kleiner werdender Stablast aufrecht erhalten werden. 
Bei konstant gehaltener Stablast schlagt daher die Stabmitte durch. 

Fir «1/4 =z verschwindet die Durchbiegung in der Stabmitte und die Stiitzkrafte werden 
Q= 0, wie es in Abb.13b dargestellt ist. Fiir Belastungswerte «1/4 <a sind zur Aufrecht- 


l 
Abb. 12. x, :1/2 in Abhingigkeit von = nach den Bedingungen a) und b). 


Abb. 13a. Stabform Abb. 13b. Stabform Abb. 13c. Stabform Abb. 14. Der an beiden 
1 3 al al Enden eingespannt ge- 

cosa eed fir — = 47. ' fir — <7. lagerte Stab mit 
4 2 4 4 3 Zwischen-Stiitzstellen. 


7 
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erhaltung der Durchbiegung (y),—.,= f Stittzkrafte Q<0 erforderlich, wie es in Abb. 13c an- 
gedeutet ist. Wie sich noch zeigen wird, kann diese Stabform noch fiir anders angeordnete 
Stiitzkrafte existieren. 

Im Ablauf des Durchschlagvorganges zeigt sich der bisher einzige Unterschied gegeniiber 
dem Stab mit beiderseits gelenkig gelagerten Enden. Da der beiderseits eingespannte Stab 
nach dem Durchschlagen nicht als aus drei Teilstaben vom gleichen Typ bestehend aufgefaBt 
werden kann, ist anzunehmen, daB sich weitere Unterschiede im Verhalten zeigen, wenn unter 
Zugrundelegung der Stabform und -belastung nach Abb. 14 die Untersuchung noch ein Stiick 
weitergefiihrt wird. 

Fiir diesen Fall werden die Biegungsmomente im Stab 


M,= Py, — Qx,— M,, 
M, = Py, — Qx, + Qo (%, — %) — My 


und die allgemeinen Liésungen von (1) lauten 


y, = A, sin ax + B, cos ax bp o* | oe 
= Q M, Q. 
Yo = A, sin ow + B, cos ax ++ > ax f Pp Pao (* Ke)is 


Daraus wird mit den Grenzbedingungen fiir die Stabquerschnitte an den Stellen x = 0,, 
x= x), und x = 1/2 nach Elimination von A,, B,, A,, B,, das folgende Gleichungssystem ge-? 


wonnen: 
| 


0O=—=—ft+ ee (0%) — sin 0X) + i + Mo (1 — cosa), 

m4 Bhat lela) Soule — 9 1 ae 
| J 

— * 4 © (1 — cos a m) + * + MP sin ay 


adm cae) acces cae 


Dieses System kann nur befriedigt werden, wenn seine Koeffizientendeterminante verschwin-} 
det. Das ist der Fall entweder fiir 


ce) t=na2 (2 eee 


oder 


Die Bedeutung dieser beiden Lésungen wird nachfolgend naher untersucht. 


aX, : 
c) cr =na(n=1,2,3...), Mit n= 1 folgen aus dem Gleichungssystem die bereits be~! 
kannte Beziehung (10) und die Ausdriicke 


a1\ a ed 
Qu 1? ey miet is Ti. 


” ee Ts 
a2 Ry ee aad eee (14)) 
peer In ot | 
al «xl | 
eal ET (15) 
16 EJ ioe | 
re 
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Fiir die Durchbiegung des Stabes ergeben sich die beiden Gleichungen 


al 
1 ; Ce 
eae aq (o* sin « x) al _al (1—cosax)], 
gee Tha ke (16) 
f i al al 
Oe ie Tt ont OX sin ax + tg; cosax— (7p +2]. 


aE ay eae IC 


Die untere Grenze des Giiltigkeitsbereiches der vorstehenden Formeln ist gegeben durch das 
Verschwinden einer der Stiitzkrafte Q) oder Q,. Es wird nach (14) Q, = 0 fiir «1/4 = 1,46 2 
und nach (10) und (14) Q, = Q,— Q= 0 fiir a 1/4 = 32/2. 
Mithin ist die untere Grenze des Giiltigkeitsbereiches ge- 
geben durch a 1/4 = 3 77/2. 

Die obere Grenze des Giiltigkeitsbereiches wird erreicht, 
wenn die Stabkriimmung an einer der Anlagestellen ver- 
schwindet. Das ist der Fall fiir «1/4 = 2. Dort ist nam- Al| % gs 
lich, wie sich nach (16) leicht nachpriifen 14Bt, (y,'),—., y 
= (¥o'")e=x, = (Yo)x=12 = 9. Die beiden Stabformen fiir 
o1/4—= 37/2 (untere Grenze) und «l/4= 2 2% (obere Grenze) 
sind in Abb. 15a und 15b dargestellt. Die Stiitz- und 
Auflagerkrafte betragen hierbei 


; 3 3 
roe ae Ql 


—s A == = 772 
Ao Sh ETA 82H Is ; 
Q, = 0, Mo Ef, 
ee ee ats eee 2 Oe . 
ee 7. song, ° a ETF OR 


Beige MES 0 
32 BSF . S : 


My 


Wird die Stablast soweit gesteigert, dab a 1/4 > 20 ces aot iee 
wird, so bilden sich sowohl fiir die A4uBeren Stabteile von Stabformen nach Bedingung ) 
der Lange 3/8-1, als auch fiir den mittleren Stabteil von fiir ~ =>, fiir = 2m 


der Lange 1/4, ungekriimmte, an den Stiitzflachen an- 
liegende Stabbereiche aus, deren Linge mit wachsender Stablast zunimmt. 


toe 
M, 2 
EO = — if : = A, 5 j= tp nes B, 5) 
Poa OX oO Xo 122 Rory by 
2(ts = a: Al Sac ooN 
(17) 
al a Xo % Xo 
Qo if 4 Q, = Af 2 2 , 
Po l Xo ONG iP ei 1% A %y  O Xp 
20(¢—3] cn ee 20(q ye ae 
Die Biegelinie des verformten Stabes verlauft nach den Gleichungen 
v4 f [sm ax — ax + tg “5 (1 —cosaa)), 
ae pate es ele) 
85 2 


_mb(b+e—s)]] | 
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Die untere Grenze des Giiltigkeitsbereiches ergibt sich aus Q, = 0 zu a 1/4 = 2. Die Kriim-- 
mungen des Stabes an allen Anlagestellen verschwinden bei der gleichen Stablast. Hieraus er-- 
gibt sich die obere Grenze des Giiltigkeitsbereiches zu a 1/4 = 2 a. 

Die Stabformen fiir a1/4 = a (untere Grenze) und 71/4 = 22 (obere Grenze) sind in Abb. 16a1 
und 16b dargestellt. Die zugehorigen Stiitz- und Auflagerkrafte betragen 


Nae Oil leek Si Ui wer Ge COs eee Mone ot 
fir a: O05) py age Ey won eee Pe ee 
eee! Q 3f Q a4 Q iD 
ae, =? — — ——— —— 
fir py s go ee Pa 2x’ M, = 0 


y, 
ty 


ee gf eel 


i 
SS 
ae 
S 
X 
Ss 
| 
ib 


ie 


= 


= zy O13-1 
-—- 0.358-1 


Abb, l6a. Abb. 16b. Abb.17. Eine weitere, Abb. 18. Der.an einen 
Stabformen nach Bedingung d) miégliche Stabform Ende gelenkig, am ande. - 
gee ae! al 1 ren Ende eingespannt ge- - 
fiir ite Ty fiir ae 22. fiir oh rte 1,43 2. lagerte Stab. Stabform 


fir 1,432 <al < 3a. 


Zu den Bedingungen c) und d) gehért die gleiche obere Grenze der beiden Giltigkeits-} 
bereiche. Fiir diese Grenzbelastung al1/4 = 2% stimmen sowohl die Stabform als auch die} 
Stiitz- und Auflagerkrafte nach den Bedingungen c) und d) iiberein. . 

Neben den hier behandelten Stabformen ist noch die nach Abb. 17 méglich. Zu ihrer Auf-- 
rechterhaltung gehért wenigstens eine Stabkraft, die dem Wert ol/4—=1/2-1,437 entspricht.. 
Der Stab kann aus 2 Teilstaben zusammengesetzt aufgefaBt werden, die an einem Ende zal 
lenkig, am anderen Ende eingespannt gelagert sind. Auf die Untersuchung dieser Stabform) 
kann hier verzichtet werden, weil sie gegeniiber der im nachsten Abschnitt durchgefiihrten Be-- 


handlung des an einem Ende gelenkig, am anderen Ende eingespannt gelagerten Stabes nichts: 
Neues bringt. | 


4, Der an einem Ende gelenkig, am anderen Ende eingespannt gelagerte Stab. Hierfiir; 


sree sich mit den Bezeichnungen der Abb. 18 die Biegungsmomente in den beiden Stab-- 
teilen zu 


M, = Py, — Q,%, — M). 


Damit lauten die allgemeinen Lésungen von (1) 


Vy aa A, sin aX, + B, cos OX e Xy : ‘ 

mit 2 
— 31 M, . 7° | 

Vo A, sin «x, + B, cos am, + 2 xy | EJ 
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Mit Hilfe der Grenzbedingungen fiir die elastische Linie, die in diesem Fall auch die Gleichheit 
der Biegungsmomente fiir die beiden Stabteile an der Anlagestelle enthalten, ergeben sich nach 
Elimination aller iibrigen Groen die beiden homogenen Gleichungen mit Q, und Q, 
Q1(tg 0 — ah) —2.0,(te3* —S) = 0,| 
. (19) 
Q, 01, + O,(te?— 258) =0. 


Diese Gleichungen kénnen nur dann eine von der trivialen Nullésung verschiedene Lésung 
haben, wenn die sich aus dem Verschwinden der Koeffizientendeterminante ergebende Be- 


dingung : 
fae Ip 
oe toed 1 (20) 
tgal, t al, 
are} 


erfiillt wird. Mit Hilfe von (20) lassen sich bei gegebenem a und / die Teilstablangen I, und 1, 
bestimmen. Die Integrationskonstanten, Stiitz- und Auflagerkrafte ergeben sich zu 


l Q is 
i pe fieos «1, i= Wee 
teal «l,? eae Po tgol, —oal, ’ 
al, 
tg = 21 
Q. s ip ae | My fte 2 a B ( ) 
Pia ae al, al, 2° _ of te ete os a 
a 2 os ae 


Die untere Grenze des Giiltigkeitsbereiches fiir die Stabform nach Abb.18 ist durch die Be- 
dingung Q = Q,-++ Q, = 0 gegeben. Damit wird aus (19) 


aol,—tgol=—22, tg it 2h 
mit den kleinsten Wurzeln 
lp 0,5 1305 pS Wg 99 
eee ee ee 


Dieser Wert entspricht der Knickbedingung eines an einem Ende gelenkig, am anderen Ende 
eingespannt gelagerten Stabes von der Lange I. 

Die obere Grenze des Giiltigkeitsbereiches der in Abb. 18 dargestellten Stabform ergibt 
sich aus dem Verschwinden der Stabkriimmung an der Anlagestelle. Das liefert fiir die beiden 
Stabteile 


t al, 
tga l, ay 2) =) 
tgal,—al, 2 al, al, 
o3 g 
mit den kleinsten Wurzeln, 
oh=%, \yi—3q. (23) 


e1,=27, 


Die Abb. 19a und 19b zeigen die Stabform an der unteren (al = 1,43 x) und oberen Grenze 
(al = 3m) des Giiltigkeitshereiches. Die Stiitz- und Auflagerkrafte betragen 


é . Q, 18 map Q, 1 Ss i. Se La Oe ee 
ee Lana 2" 1,43° E J f aN 143° ES f ’ 
MP. 
L432 Bf 2 
# = ; Ona tee Qe 1 ee ore AGL ==) 
fair ob 37: SLES Sea. ’ CIN TS ame ele 2 PT Dg es : 


Fir Stabbelastungen «1 > 32 werden mit den Bezeichnungen der Abb. 20 die Biegungs- 


momente in den vier Stababschnitten 
M,= Py,—Q,%, M, = Py, — Q,%, +4, Q (%1— %0), 
Ma ie 0, ty ao My, 2 = Py¥,— Q.%, + Ay Q (%, — %y9) — My» 


| 
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und damit die allgemeinen Liésungen von (1): 


¥, = A; sin o% + B, cos ax, +S! my , i| 


- ets. Ls eh One 
Vy, = Ay st Oke B, cos a%, + S%, Pp (4; — X19), 


: Q M. 
¥_ = A, sin ax, + B, cos wx, + 5 % + —’> 


= a ae ke ON Oe 
¥, = A, sin ox, + B, cos az, + Sz, + ra = (%_ — Xpp) - 


Abb, 19a. Abb. 19b. Abb. 20, Mégliche Stab- 
Stabformen fiir al = 1,432, fiir al = 3 2. verformung fiir al > 32. 


Mit den Grenzbedingungen yy= 0 fire x, = 4,9. 91 = 0 fir a =), V7 — 0 far 4, ee 
¥2 = 0 fiir x, = 1, ergeben sich daraus die Gleichungen 


7 i : 
A, 4 cos 6%) — B,a sin ey 
A, acos al, == B, «sin « 1, te =o 
und 
~ eae 0 
Ay & COS Xp — B, a sin O Xing = — 
A,acosal, — Basinal, = 28 S | 
Wegen 1,0 + A,Q = Q,+ Q, kénnen diese Gleichungen nur dann befriedigt werden, wenn eine | 
der beiden Doppelbedingungen 
a 4Q@—Q,=0,| und [ A = B,=0 | 
|4.9—@=0,| |4 2 = B,=0 | 
nye Q,=0,| und }%(l,— %o) = na (n= 1, 2, 3...) | 
|\A,9—Q=090 | O(l, — Xo) = nw (n = 1, 2,3...) | 


erfiillt ist. Die Bedeutung dieser beiden Lisungen wird nachfolgend untersucht. 


a) /419—= 0 | und 


1 cesar | 
ie a . Mit dieser Bedingung ergeben sich aus | 
|, Q— Q, = 0 | | 4, = B,= 0 
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dem mit Hilfe der noch nicht benutzten Grenzbedingungen gewonnenen, aber hier nicht auf- 
gefiihrten Gleichungssystem die Beziehungen 


Be, Sea f= x, Geax Es, Ae B,=0, 
O Nyy = 270, oe. Q=La ES, ApS B, = M,= 0. 


Die Gleichungen der elastischen Linie lauten im Bereich 0 < eee 
Dal == _ (sin OX, -- AX;) und bel = i 
im Bereich 0 <x, <1, 
fae SIN AX, + OX) und be yeokd [ie 


102 0290) (a(l, — x9) = 

b))41 1 und 1 10 

Lom ol= 0) | X(lp — 9) = J 
Diese Bedingung ist gleich fiir den Fall n = 1 an- 

geschricben. Hiermit ergeben sich auf die gleiche 


Zi _ X29 


Lk 


x Xa. he : Abb, 22. Stabform Abb. 23. Stabform 
Abb. 21. i und 7% in Abhangigkeit von al, und al, nach den nach Bedingung a) nach Bedingung b) 
o Bedingungen a) und b). fiir al = 52. fir xl = 3,587 2. 


Weise wie bei a) ermittelt, fiir den Verlauf der Biegelinie die Gleichungen 


i sin Ax, ‘ 
Ns tZ AX 9 —AX yy \ cos aX oat ‘ 
f } im, Bereich 0 = 4, F. 
ies tg XX) — A219 {tE AAyq cos [a(%, — xXo)] — A249} » 
if , DB King 
v= ae Aa: sin ox, + tg 5 (1 — cos ax,) — ax, 
20  %% 
2(te 2 : im Bereich 
a aX = F OR ale 
Vo = me i as te i [1 + cos [aw (%,— xy9)]] — ay} yee ce 
2 (ts a ro 
Fur den Zusammenhang zwischen x,) und %, gilt die Gleichung 
O Hoo 
eit fe es is (24) 
tg AX19 tg asee 


Zur Ermittlung der Grenzen, zwischen denen die Bedingungen a) und b) giiltig sind, sind 
in Abb. 21 die Verhaltnisse x1,/l, baw. %9/l, in Abhangigkeit von « 1, bzw. al,, wie sie sich ent- 
sprechend den Bedingungen a) und b) ergeben, dargestellt. 


250 Link: Uber den geraden Knickstab mit begrenzter Durchbiegung. Ingenieur-Arch i 


Die untere Grenze der Bedingung a) ist gegeben durch (23). Der Stab hat dabei die | 
Abb. 19b dargestellte Form. Die obere Grenze fiir die Bedingung a) betragt 
al=20, | 

al—=32, | 

Die hierzu gehérige Stabform zeigt Abb. 22. Die Stiitz- und Auflagerkrafte sind dabei 

fe meh rcs 
2, = a5) @ Ee 0, = Ww? > an M,= 0. 

In diesem Zustand entspricht die Stablast der Eulerschen Knicktast eines beiderseits eingespann)| 
ten Stabes von der Linge 2-1/5. Der mittlere Stabteil schlagt durch. Dieser Vorgang wird) 
in seinem Ablauf durch die Bedingung b) bestimmt, deren obere Grenze mit der der Bedingung: 
a) itbereinstimmt. Die untere Grenze fiir die Bedingung b) ergibt sich aus A,Q=0 zu al=3,587a1, 
Die zugehérige Stabform zeigt Abb. 23. Die Stiitz- und Auflagerkrafte betragen hierfiir 


ae 3,587 
g, =f Sita eI Oy O20 any | Mp eee 


Si = Bas 


Ey 


Fir « 1< 3,587 a wird entweder A,Q < 0 oder an der Stelle x, = I, greift eine weitere Stiitzkraft} 
an, damit das Gleichgewicht der duBeren Krafte am Stab gewahrleistet ist. Da negative Stiitz-| 
krafte oder auch positive Stiitzkrafte an Stellen mit Durchbiegungen y <|f| ohne praktische 
Bedeutung sind, wird von einer Weiterverfolgung des Verhaltens des Stabes bei Wertem 
al< 3,587 abgesehen. 


5. Zusammenfassung. Die Untersuchung des durch ebene, starre Flachen in seiner Durch-: 
biegung begrenzten geraden Knickstabes mit Hilfe der elementaren Biegetheorie fiihrt auf eine: 
interessante Durchschlagerscheinung. Unabhangig von der Art der Lagerung der Stabendeni 
stiitzt sich nach Uberschreitung der zugehérigen Eulerschen Knicklast der Stab an der Stelles| 
seiner gréBten Durchbiegung gegen eine der Stiitzflachen ab. Die sich dabei ausbildende| 
Stiitzkraft wachst mit zunehmender Stablast an, wahrend sich gleichzeitig die Stabkriimmun | 
an der Anlagestelle verringert. 

Erreicht die Stabkriimmung an der Anlagestelle die Kriimmung der Stiitzflache, dannlj 
bildet sich bei weiterer Steigerung der Stablast ein an der Stiitzflache anliegendes, ungekriimm--, 
tes und daher momentenfreies Stabstiick aus, an dessen Enden Einzelkrafte von der Stiitzflache?) 
auf den Stab iibertragen werden. Die Lange dieses Stabteiles wachst mit der Stablast solange i 
an, bis diese die GréBe der Eulerschen Knicklast des als an beiden Enden eingespannt an-- 
zusehenden geraden Stabteiles erreicht. Der nun einsetzende Knickvorgang des bis dahin)) 
an der Stiitzflache anliegenden Stabteiles kann wegen der zunehmenden Knicklange mit sick }) 
verringernder Stablast aufrecht erhalten werden. Bei unveranderter Stablast schlagt daher} 
der mittlere Stabtei] durch und kommt an der gegeniiberliegenden Stiitzflache zur Anlage. 

Bei weiterer Steigerung der Stablast kann der Stab immer als aus Teilstaben zusammen: - 
gesetzt aufgefaBt werden, deren Lagerung den hier untersuchten drei Grundfallen entspricht.. 
Es ist daher nicht zu erwarten, daB eine Weiterfiihrung der Untersuchung noch zu grundsatz- - 
lich neuen Erkenntnissen fiihrt. 


(Eingegangen am 6, Oktober 1953.) 
Anschrift des Verfassers: Heinz Link, Oberhausen (Rheinland), Goebenstr. 153. 
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Zur Analyse der Daimpfung freier Schwingungen. 
Von H. Kauderer, 


1. Einleitung. Bei mechanischen Schwingungen fester Kérper setzt sich die gesamte 
Dampfung meistens aus mehreren, durch verschiedene physikalische Vorginge hervorgerufenen 
Kraftwirkungen zusammen, die an dem schwingenden System gleichzeitig angreifen. Die 
wichtigsten hiervon sind bekanntlich die ihrem Betrag nach als konstant angenommene Kraft 
der trockenen Reibung, die der Geschwindigkeit proportionale sog. Fliissigkeitsdampfung und 
die bei Schwingungen in gasférmigen Medien auftretende, dem Quadrat der Geschwindigkeit 
proportionale sog. Luftdimpfung. Es besteht nun bei technisch ausgefiihrten Schwingern 
haufig der Wunsch, aus dem (etwa durch ein Oszillogramm aufgenommenen) Ablauf des freien 
Schwingungsvorgangs zu erschlieBen, welcher Teilbetrag an der gesamten Daimpfungskraft 
auf jede der einzelnen Kraftwirkungen entfallt, d.h. also mit Hilfe des Bewegungsablaufs den 
ganzen Komplex der dampfenden Krifte zu analysieren. Zur Lésung dieses Problems hat 
G. Plato! eine Formel entwickelt, welche in allgemeiner Gestalt die Dampfung als Funktion 
der Geschwindigkeit mit Hilfe eines uneigentlichen Integrals naherungsweise wiedergibt. 
Da jedoch die numerische Auswertung von Versuchsergebnissen mit Hilfe dieser Formel 
schwierig ist, soll hier die Aufgabe erneut angegriffen und auf einem anderen, in den meisten 
praktischen Fallen bequem begehbaren Weg behandelt werden. 

Wir beschranken uns auf einen einfachen freien Schwinger mit gerader Kennlinie, von dem 
wir voraussetzen, dafi sein Ausschlag q(t) der Differential gleichung 


q+ w*q+ €g(q) = 0 (1) 


mit der Konstanten w geniige. Das ,,Dampfungsglied“ ¢ g(q) bestehe aus der dimensionslosen, 
gegen Kins kleinen Gidf8e ¢ (d.h. wir setzen eine ,,schwache‘“+ Dampfung voraus) und einer 
allein von der Geschwindigkeit g abhangigen Funktion g(q), die sich durch eine Reihe von der 
Form 
(4) = cong + ed + eldld + esld?G+ ++» = > 4 lal" (2) 
v=0 

mit noch unbekannten Koeffizienten c, darstellen lasse. Wir wollen zeigen, wie man aus der 
Folge der Amplituden eines Schwingungsvorgangs und aus seiner Kreisfrequenz die Koeffi- 
zienten C,= ec,(y = 0,1,2...) naherungsweise ermitteln und somit die Dampfungskraft 
analysieren kann. 

Um dieses Ziel zu verfolgen, leiten wir zunichst (Ziff. 2) mit Hilfe eines dem Verfahren von 
N. Kryloff und N. Bogoljuboff? zugrundeliegenden Ansatzes die Grundgleichung her, die im 
Sinn einer Naherung erster Ordnung fiir schwache Dampfung einen Zusammenhang herstellit 
zwischen den Dampfungskoeffizienten c, und den Werten benachbarter Amplituden. (Hine 
exakte Untersuchung dieses Zusammenhangs ist bekanntlich nur in einigen Sonderfallen még- 
lich.) Aus der Grundgleichung wird (Ziff. 3) eine Methode zur Berechnung der Dampfungs- 
koeffizienten entwickelt, die auch noch auf den Fall erweitert wird, da8 man nicht unmittelbar 
aufeinanderfolgende Amplituden, sondern nur etwa jede m-te gemessen hat. An einem Zah- 
lenbeispiel (Ziff. 4) wird das Verfahren erlautert. SchlieBlich zeigen wir noch (Ziff. 5), wie man 
direkt aus den gemessenen Amplitudenwerten eines Schwingungsvorgangs ohne Rechenauf- 
wand ein naherungsweise richtiges Diagramm fiir den Verlauf der Dampfungskraft als Funktion 
der Geschwindigkeit gewinnen kann. 


2. Die Herleitung der Grundgleichung. Wir betrachten zunachst den Bewegungsvorgang 
wahrend einer einzigen, zwischen den maximalen positiven Ausschlagen g, und q,,1 verlaufen- 
den Schwingungsperiode, die, von einer beliebigen Periode an gezahlt, als k-te bezeichnet sei. 


1G. Plato, Z. angew. Math. Mech. 25/27 (1947), 5S. 93. 
2 N. Kryloff und N. Bogoljuboff, Introduction to non-linear mechanics, S. 8, Princeton 1947. 
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Fiir diese setzen wir den Ausschlag q(t) mit der Konstanten 


Ap = 5 (qe + Qu41) (3) | 


und den noch unbekannten Funktionen Q(t) und g;,(t) in der Gestalt 
q(t) = [Ap + € Q:(t)] cos [wt + € qx (t)] (4) 


an, wobei wir zugleich zum Ausdruck bringen, daB sich die zu erwartende Schwingungsform 
wegen der (infolge des Faktors ¢) geringen Dampfung auch nur wenig von derjenigen freien 
Schwingung q(t) = A, cos mt unterscheiden soll, deren Amplitude gleich dem Mittelwert aus 
Anfangs- und Endamplitude der k-ten Periode der gedimpften Schwingung ist. 

Da wir z wei willkirliche Funktionen Q, und gy, gewahlt haben, so diirfen wir fiir diese auBer 
der Forderung, da der Ansatz (4) die Bewegungsgleichung (1) befriedigen mu, noch eine 
zweite stellen. Sie bestehe darin, da wir auch noch vorschreiben, wie sich die Geschwindigkeit 
q(t) darstellen lassen soll. Wir verlangen namlich, daB 


q(t) = — w[ A, + € Qy(t)] sin [wt + € px (4)] (5) 


werde, wobei wir also wiederum eine nur geringe Abweichung vom Geschwindigkeitsverlauf | 


der ungedimpften Schwingung voraussetzen. 
Wir kénnen nun aber auch andererseits aus (4) durch Ableitung nach der Zeit einen Aus- 


druck fiir q(t) gewinnen. Der Vergleich mit (5) wird dann eine Beziehung liefern, die zwischen | 


den Funktionen Q; und g, (bzw. ihren Ableitungen) bestehen muB. Wir vernachlassigen hierbei 


simtliche Glieder, die in ¢ von hoherer als der ersten Ordnung klein sind, gehen also nicht — 


direkt von (4), sondern von einem Ansatz aus, der sich aus(4) ergibt, wenn man die rechte Seite | 
nach Potenzen von ¢ entwickelt und sich dabei auf die in ¢ linearen Glieder beschrankt. Dieser 
Ansatz wird 


q(t) = A, cosmt + e [Q,(t) cos wt — Ay yy (t) sin wt]. (6) 
Leitet man ihn nach der Zeit ab, so ergibt sich 
q(t) = — A,w sin wt + «€ (Q, cos wt — w Q, sin wt — A, G, sin wt— A,wg,cosmt). (7) 


Um diese Darstellung mit (5) vergleichen zu kénnen, miissen wir auch (5) in ¢linearisieren. Wir 
erhalten so 


q(t) = — A, @ sin wt — ew(Q,sinwt+ A,g, cost). (8) 

Damit sich nun die Gleichungen (7) und (8) in EKinklang bringen lassen, ist zu verlangen, daB 
zwischen Q;(t) und ;(t) die Beziehung 

Q:(t) cos wt — Ay Gi (t) sin wot = 0 (9) 


besteht, die fiir die in ¢ linearisierte Theorie der Forderung, da8 q(t) die Form (5) annehmen 
soll, gleichwertig ist. 


Mit der aus (8) gebildeten zweiten Ableitung 
q(t) = — A,w? cos wt — ew (Q, sin wt + wQ,cosmt-+ A, Pecos Wt — Aj, © x sin @ t) 
und mit q(t) aus (6) geht jetzt die Bewegungsgleichung (1) tiber in 
Q, sin ot + A, GP, cos Wt = ~ g(@). 


Lassen wir — wieder im Sinn einer linearisierten Theorie — die mit ¢, e2,... behafteten Glieder 
weg, die sich in g(q) ergeben, wenn wir in (2) mit dem Ansatz (8) eingehen, so diirfen wir g(q) 
durch g(— A; @ sin w t) ersetzen und erhalten die Beziehung 


Qx(t) sin wot + Ay Px (t) cosmt= = g(— A, @ sin wt), (10) 
die zusammen mit (9) von den beiden Funktionen Q,(t) und @,(t) unseres Lésungsansatzes be- 
friedigt werden mu. Durch Multiplikation von (9) mit cos @t und von (10) mit sin mt und 


Addition erhalten wir 


&(—A,@ sin wt) sin wt, (11) 


e[+ 


Ox(t) = 


} 
| 
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ud durch Multiplikation von (9) mit — sin wt und von (10) mit cos wt und Addition ebenso 
: 1 
AL gelt) = es &(—A, @ sin wt) cos mt, (12) 
wobei wir jetzt, wenn wir das Vorzeichen des Arguments von g im Bereich 0 < wt < 2 a be- 
riicksichtigen, gemaB (2) 

— (¢o + ¢, Axo sinwt-+ c, Aj w? sin? wt + cz Apo sindwt + -- 

fiir 2 <t<(2y4+ 1) = 
&(—A,o@ sin wt) =4 : ; 
+ Cy — ¢, Ay wsin wt + c, Aj w* sin? wt — c, Ap w sin? wt + — 
fiir (2 w+ 1) 7 <t<2(w+1) 7 | 


mit w= 0,+ 1, + 2,... zu setzen haben. Es werden also sowohl On als auch gy; periodische 
Funktionen von ¢ mit der Periode T= 2 a/w. 

Jetzt ist die Linge T;, der k-ten Periode der gedimpften Schwingung zu berechnen. Wer- 
den ihr zeitlicher Anfangs- und Endpunkt mit 4 und t,1 bezeichnet, so mub 


qk = Q(te)> We+1 = Q(tey1) und (th) = q(tey1) = 0 
werden. Aus der letzteren Forderung folgt wegen des Ansatzes (5), daB die Werte von 


wt + eg,(t) fir t= t und t= t,,1 ganzzahlige Vielfache von 2 sein miissen, die sich um 2 7 
unterscheiden, d. h. dai die Periodenlange 


)] 
| oe, 


thay 
1 20 : 
Th = tpi — t= — {2 %— € [pe (ters) — gu (te)]} = =~——= | palta (14) 
"e 
werden mu. Setzt man den Wert von q;(t) aus (12) in das Integral ein: 
+1 kt +1 
: 1 : 1 : ; 
| pxit) dt a | g(—A;, @ sin wt)cosat dt = aie | g(—A, sinw t) d(sinw t), 
th th, A=ipe 


so erkennt man aus der Form von g(— A; @ sin w t) nach (13), daB® dieses Integral verschwindet, 
wenn t und t,,; zwei benachbarte Vielfache von 2 z/m sind, so daB T, = 2 2/w wird. Mit 
diesem Wert von TJ), wird nun aber auch die Gleichung (14) befriedigt, und dies bedeutet, daB 
bei der hier vorausgesetzten nur von der Geschwindigkeit abhangigen und von erster Ordnung 
kleinen Dampfung die Kreisfrequenz der Schwingungen den Wert  beibehAlt, den sie bei feh- 
lender Dampfung angenommen hitte. 
Wir denken uns jetzt den Beginn der Zeitmessung so gewahlt, da 
t= 2kalw, tyy1 = 2 (k+1) alo 
wird und berechnen fiir diese Werte die Differenz der beiden benachbarten Amplituden q, und 
qk41- Da nach (4) wegen des Verschwindens von q(t.) und q(t, +1) 
2(k+1) n/w 
Ge — MU+1 = € [Qu (te) — Ox (tey-1)] = — € Q;,(t) dt 
2k alo 
werden muB, so ergibt sich nach (11) mit (13) fir »—k 
2(k+1)2/o 
g (— A, wsin @ t) sin w t dt 
2kxn/w 
(2k+1)2/o 
(cysinwt+ c, A,wsin?wt+ c, Agw*sin'wt + -- -) dt — 
2ka/w 
2(k+1) 2/0 
(cy sin wt — ¢, Ay @ sin? wt-+ c, Afw* sin’ wt— + ---) dt 
(2k+1)2/o 
n/2 ow 
(cy sin wt + c,A,w sin? wt + c, Ajw? sin? wt + ---) dt 
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oder schlieBlich, wenn wir die Integrationen ausfithren, die Grundgleichung 


; uf 2 
Gk — Wke+1 = 4 C9 + 2 6 Ayo +3 % vy sie = 


1.3. 
2 


| 
| 


5-++(2»—1) HEA 2v—1 9y—I1 2480 2v A2v 2 ether 
SWRUTURS) OL ee” r 5.5 fs--(29 il) ae ee i : 


die gemeinsam mit (3) den Zusammenhang zwischen den Amplituden und den Dampfungs- 
koeffizienten herstellt. 

3. Die Berechnung der Dampfungskoeffizienten. Um mit Hilfe der Grundgleichung (15) 
die Dimpfungskoeffizienten zu berechnen, fithren wir in ihr als Abkiirzungen die Summe 
zweier Nachbaramplituden 


Xp = 2A, = Qe + G41 (16) 


und ihre Differenz 
Vk = Wk — Wk+1 (17) 
sowie durch die Gleichungen 


(2 ps Asa) ERA Tiers | 
wo? ag a 1 a ons OO) (a) 18) 
Ely = Cy= 70> Bes see 2\9—2 “ — 
fils ppm abhony carer! hg 0} eee) c 
t 2-AeOses y (a) 


(15) | 


die reduzierten Daimpfungsfaktoren 7, y,, Yo,--- ein. Die Grundgleichung (15) 1a8t sich | 


dann in der Form schreiben 


Ye = f(x) » (19) 


wobei 


f(x) =r net yeti Dye (20) | 


eine Potenzreihe bedeutet, deren Koeffizienten die reduzierten Dampfungsfaktoren sind. 

Die praktische Aufgabe besteht nun darin, da man sich die experimentell gewonnenen 
Amplituden qo, q,,--+> qk, --- sowie die Kreisfrequenz @ vorschreibt, hieraus die Summen 4% 
nach (16) und die Differenzen y;, nach (17) berechnet und nun die Koeffizienten der Reihe 
fiir f(x) so bestimmt, daB die Gleichung (19) fiir simtliche Wertepaare (x, y,) erfiillt wird. 


Wir werden demnach bei N + 2 Versuchswerten qo, q,,---4nvs, und somit N-+ 1 Werte- | 


paaren (%, Yo), (%1 ¥1)>+++> (%n, ¥w) die Bedingung (19) schon mit einer bei y= N ab- 
brechenden Reihe (20) fiir alle (x;, y,) befriedigen kénnen. Die reduzierten Dampfungsfak- 


toren y, kann man hierbei etwa mit Hilfe einer der bekannten Interpolationsformeln von : 


Newton oder von Lagrange berechnen. 


Meistens wird man aber von vornherein voraussetzen, da sich die Dampfungsfunktion | 
g (q) durch ein Polymon mit einer bestimmten, endlich groBen Zahl von n + 1 Koeffizienten — 


Co: Cy, +++ > Cy Wiedergeben JaBt. (Oft wird man sich dabei auf n = 2, also den Ausdruck 
s@ =o Gt ait alala 
beschranken.) Liefert nun der Versuch mehr Wertepaare (xp, y,) als die Anzahl (n+ 1) 


der gesuchten Dampfungskoeffizienten betragt, so werden wir im allgemeinen die Gleichung (19) 
nicht mehr fiir simtliche Versuchswerte streng befriedigen kénnen. Wir suchen dann wenigstens 
die Faktoren y, so zu bestimmen, daB die Gleichung (19) fiir alle (%,, yx) méglichst gut er- 
fillt wird. Dies geschieht am zweckmaBigsten mittels der Methode der kleinsten Quadrate. 
Verlangen wir, da die Summe der Quadrate der Abweichungen der N-+ 1 Funktions- 
werte f(x.) (k = 0,1,...,.N) von den experimentell vorgeschriebenen Werten y, ein 
Minimum werde, so miissen bekanntlich die Dampfungsfaktoren y, das folgende System von 
linearen Gleichungen befriedigen: 


n N N - 
= by me) Y= 2 (xt yn) (w= 90,1,2,...n), (21) 


wobei die Summen iiber k sich auf simtliche aus den beobachteten Amplituden berechneten 
N + 1 Wertepaare (xo, Yo), .- +, (xn, Yn) (N> n) erstrecken. Es sei noch darauf hingewiesen, 
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dai diese Wertepaare auch von verschiedenen Amplitudenfolgen stammen diirfen, die zu 
verschiedenen Anfangswerten q) gehéren. Man kann so durch eine hinreichend gro®e Anzahl 
von MeBreihen eine beliebig dichte Folge von MeSpunkten im (x, y)-Diagramm herstellen. 

Bisweilen ist es, besonders bei sehr schwach gedimpften Schwingungen, aus meBtech- 
nischen Griinden nicht méglich, unmittelbar aufeinanderfolgende Amplituden so genau zu 
beobachten, daf ihre Differenzwerte y, noch zuverlissig sind. In diesem Fall beniitzen wir 
nur die Amplituden nach je m Vollschwingungen (m> 1): 


; o> Imo Fom> > ++ > Wk ms Wk+1) m Cr) 
setzen jetzt aber in (4) fiir die Konstante A, nicht den Mittelwert der Nachbaramplituden, 
sondern denjenigen hintereinander gemessener Amplituden q, , und q@+41)m ein und nehmen 
sinngemaf} an, da® dieser Ansatz zwischen der m k-ten und der m(k + 1)-ten Vollschwingung 
bestehen soll. Die (in diesem Fall formal gleich lautende) Gleichung (15) denken wir uns 
nun m-mal, und zwar fiir die Differenzen gm — Qmk+1> Imk--1 — Umk+20+++> Ym(k+1)—1 — Im(k+-1) 
angeschrieben. Die rechte Seite wird hierbei stets dieselbe Form haben wie in Gleichung (15). 
Addieren wir diese m Gleichungen, so erhalten wir links die Differenz x, — Xm(k-+ 1) und 
rechts das m-fache der rechten Seite von (15). Wenn wir jetzt in Abanderung von (16) 


und (17) 
1 
tp = 2A = Gm + Im(e+1)> Yk = — (Fmk — Im(e+1)) (22) 


setzen, so behalten die Beziehungen (18), (19) und (20) ihre Gestalt unverandert bei; wir 
kénnen sie also nach wie vor zur Berechnung der Dimpfungskoeffizienten beniitzen. Die 
Genauigkeit des Verfahrens nimmt natiirlich im allgemeinen mit wachsendem m ab. 


4, Ein Zahlenbeispiel . An einem Zahlenbeispiel wollen wir die Anwendung des Verfahrens 
zeigen und gleichzeitig die allgemeinen Formeln fiir die praktische Rechnung im Fall n = 2 
wiedergeben. 

Durch Spiegelablesung wurden an einem Pendel mit der Schwingungsdauer von 1 sek 
(also @ = 2 a sek) nach je einem Hin- und Hergang die in der zweiten Spalte der Tabelle 
im BogenmafS§ wiedergegebenen 11 Amplituden gemessen. Wir berechnen nun zuerst die 
Summen x = q, + q,41 und die Differenzen y, = q, — qx413 sie stehen in der dritten und 
der vierten Spalte. Sodann bilden wir folgende Summen, jeweils von k = 0 bis k = N mit 


ay == "9: 


- Le x} = N+ 1a) k q % ¥ 
g 7) k hb 
hee iO AON. gi |g Seales wailed hee 

a) XE = 12002 - 10-3, 0 | 0,0287 | 0,0514 | 0,0060 

x2 = 2,7788 - 10-4, 11 0.0227 | 0.0410 | 0,0044 

oe. eer 2} 0,0183 | 0,0330 | 0,0036 

2, Xp = 1,1828 - 10%; ie 31 0.0147 | 0,0266 | 0,0028 

Ye = Wo — Inqi = 90,0272, 4} 0.0119 | 0,0214 | 0,0024 

>) Xe Ve = (Go + an+1) (Go —9N+1) = 5 | 0.0095 | 0.0170 | 0.0020 

== $2144 - 10-4, 6 | 0,0075 | 0,0132 | 0,0018 

ee Axe 71 0.0057 | 0,0099 | 0,0015 

2 *f Ye = 3,1221 - 10-°. 8} 0,0042 | 0,0070 | 0,0014 

Diese Werte setzen wir in das folgende Gleichungssystem iar 0,004 0,0018 


ein, das sich aus (21) fiir n = 2 ergibt 
(N+ 1)-¥o+ 2 te Yt Xk Yo = 10 — Inti» 
Xn Yo & Xe Ya + Ys Xb Y2 = (% + IN+1) (Go — IN+1) » 
Dh Yo Ke Yat th Yn = 2s Hh Yee 


Lésen wir das System nach den drei Unbekannten yo, y, und y, auf, so erhalten wir 


¥o= 0,001113, y,=0,03307, y,=—1,19%6. (23) 
Hiermit berechnen wir die Dampfungskoeffizienten nach (18) zu 
2 2 3 


C= ZY» C,= ea ag? Ss Vie 


also in Zahlen mit w = 2 7 sek"?: 


C, = 0,010 99 sek, C, = 0,1323 sek"*, C, = 1,7964. 
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Zusammen mit den Punkten (xz, y,) ist die Parabel 
¥=Yot 1% + Yo x 
in Abb. 1 dargestellt. 
Um dieses Ergebnis mit demjenigen vergleichen zu kénnen, das man erhalten hatte, 
wenn nur nach je m= 3 Vollschwingungen die Amplitude abgelesen worden ware, setzen 
wir jetzt voraus, es seien nur die Werte q, 43, 4g. 7g Vorgegeben. Wir bilden dann nach (22) 


Xo = I + G3 = 09,0434 , 


g(G) 
sek? X= 93 + Ig = 90,0222 , 
Q06+- 9.006 %> = de + IW = 0,0105; 
1 
OD ih (qo — 43) = 0,00467 , 
Q05|- 0,005 : 
is ae (43 — qs) = 0,00240, 
O04 0,004 1 
Yo = = (Ie — %) = 9,00157. 
G03 003 Diese Werte miissen die folgenden drei 
Gleichungen befriedigen 
4 4002 Oa Yo ame % + Y2 *0 
Y= Vot MU a+ 24> 
Yo = Yo t V1 %2 + Yo ¥3> 
gor} oor | 
aus denen sich die reduzierten Dampfungs- _ 
faktoren der Reihe nach gemah 
J 0 5 ae: 1 YO Sahi Bima 6 
" Xp —%X_ \Xo— X, X— Xo 
0 505 410 14sek" ey 
Abb. 1. Pendelversuch. MeBpunkte mit Niherungsparabel und pala ee care (Xo 8 %,) ? 


Naherungsdiagramm der Dimpfungskraft. 
0 = Yo Vi XM — V2 % 
zu Yo = 1,1277, y, = 0,03310, yp = 0,001109 ergeben. Der gréBte Fehler (rund 6%) gegen- 
iiber den Werten (23) tritt somit bei y, auf. 

5. Direkte Aufnahme des Diagramms der Dimpfungskraft. Haufig diirfte es geniigen, 
fiir den Verlauf der Dampfungskraft als Funktion der Geschwindigkeit q eine graphische 
Darstellung in einem Diagramm zu erhalten. Wir kénnen ein solches Diagramm auf ein- 
fachstem Weg direkt aus dem (x,, ¥;,)-Diagramm gewinnen, wobei wir allerdings einen prin- 
zipiellen, aber im allgemeinen unbedeutenden Fehler in Kauf nehmen miissen. 

Wir gehen hierbei von der folgenden Uberlegung aus. Schreibt man die sich aus (2) er- 
gebende, fiir gq > 0 giiltige Gleichung der Dampfungskraft 


£8(9) = Ga a Gig | .C,g" > Cag caer: 
in der Form 
Ae q q\ q\* 
ef 89) = Yo tM Cre + Y ae. “eye Cs as aie ays 
so werden nach (18) die Faktoren 
Oy = 8/r = 2,546, a,=76=—2,449, og = 7128/3 x = 2,388 , an 
oy = 730 = 2,340, o,= Y1024/5n—=2,306, a6 = ¥140 = 2,279,..., 


und es 1aBt sich aus der Konvergenz des Wallisschen unendlichen Produktes fiir 2/2 schlieBen, 
daB die a, fiir y—> co gegen den Wert 2 streben miissen. Brechen wir die Reihe fiir g(q) 
nach dem Glied mit q? ab, so treten nur die Faktoren a, und a, auf. Wir ersetzen diese nahe- 


rungsweise durch den einfachen Zwischenwert 5/2 und erhalten so — allerdings nicht mehr 
streng richtig — mit den zu q bzw. ¢ g(q) proportionalen Veranderlichen 
yp 4 ; 
Secret i=) © BAG) (25) 
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die Darstellung 


m() ~ Yo + 71 € + Yo & (26) 


fiir die Dampfungskraft ¢g(q) =? 7/4. Die strenge Gleichung fiir die in & und 7 nach (25) 
dargestellte Dampfungskraft ware 


(5) = Yo + = M1 + oe Y2 5°. (27) 


Stehen uns jetzt hinreichend viele experimentell ermittelte Punkte (x, yx) im (x, y)-Dia- 
gramm zur Verfiigung, um zwischen ihnen eine Parabel hindurchlegen zu kénnen, so laBt 
sich diese Parabel bei geeigneter Wahl des Mafstabes der Achsen sofort in guter Naherung 
als Bild des Verlaufs der Dampfungskraft ¢ g(q) fiir ¢ >0 verwenden. (Den Verlauf fiir gq < 0 
erhalt man durch Spiegelung des Parabelstiicks am Ursprung.) Da aus den in (24) wieder- 
gegebenen Werten der Konstanten 4, a),... hervorgeht, daB diese ziemlich langsam kleiner 
werden, kénnen wir sogar auch noch ohne allzu grobe Fehler eine geeignete (natiirlich még- 
lichst einfache) Kurve von héherer als zweiter Ordnung als Naherungskurve nach Augen- 
maf} zwischen den Punkten (2%, y;,) hindurchlegen und erhalten auch so noch durch einfache 
Anderung der Mafstibe der beiden Achsen direkt eine brauchbare Naherung fiir das Dia- 
gramm der von q abhangigen Dampfungskraft, wobei wir aus der Streuung der MeBpunkte 
unmittelbar auf die Zuverlassigkeit des Ergebnisses schlieBen kénnen. 

Fiir das Zahlenbeispiel in Ziff. 4 weicht die Naherungsparabel fiir 7(&) nach (26), die ja 
mit der eingezeichneten nach Einfiihrung eines (£,7)-Systems identisch ist, nur um wenig 
mehr als Strichstarke von der durch die genaue Gleichung (27) bestimmten Parabel ab, auf 
deren Wiedergabe daher verzichtet wurde. 


(Eingegangen am 24. Oktober 1953.) 
Anschrift des Verfassers: Dozent Dr. H. Kauderer, Stuttgart-N, Viergiebelweg 11. 
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Gesetzmifigkeiten von Torsionsschwingungszahlen und Aufstellung 
einfacher Grenzwertformeln, 


Von H. Neuber 


1. Einfiihrung. Die kritischen Drehzahlen von mit Masse belegten biegesteifen Wellen — 
sind bereits eingehend untersucht worden, wobei sich einfache Formeln zur Bestimmung 
oberer und unterer Grenzwerte ergeben haben!; ahnliche Gedankenginge werden in vor- 
liegender Arbeit fiir das Torsionsschwingungsproblem durchgefiihrt und zur Aufstellung ein- 
facher, fiir die praktische Anwendung geeigneter Rechenverfahren herangezogen. Bei sehr 
vielen Massen ist es von grofBem Vorteil, fiir die Kigenfrequenzen obere und untere Grenz- 
werte zu kennen, um bei der numerischen Berechnung Kontrollméglichkeiten zur Verfii- 
gung zu haben. Aus der fiir das Torsionsschwingungsproblem mafgeblichen Determinante 
werden nachstehend fiir diese Grenzwerte durch besondere Umformungen einfache Formeln 
hergeleitet. 


2. Ausgangsgleichungen. Werden die in den einzelnen Wellenabschnitten (Index k) tiber- 
tragenen Torsionsmomente entsprechend Abb. 1 mit My, die dazugehérigen Federkonstanten | 
mit c,, die Traégheitsmomente mit 0, und O,,, (links bzw. rechts), die zugehérigen Dreh- 
winkel mit gy, und @z.1 gekennzeichnet, so gelten die Gleichgewichtsbedingungen 


Mi — Mi_, = Ox Pr fiir [r= iy 2 oe o7p (1) 
My=0 i Me Mn-2 | Mn-1 =0 aaa: —_ 
Cy | & fn-2 I nt My 4 My 
0, % 2 
<— Ok 
Abb. 1. Bezeichnungen beim Drehschwingungssystem mit n Massen. 


Hierzu kommen als Formainderungs- bzw. Elastizitatsgleichungen 


Mi = ce (Pe+1 — Pr) - (2) 
Wird (1) durch ©, dividiert, so folgt 


. M;, — Mr—1 
— 3 
Pk O, (3) 
bzw. mit k+ 1 statt k 
#2 Mi+1— M, 
UD hater ae k P (4) 
Wird (3) subtrahiert, so kommt 
ven Lee ee ea 
Pkt1— Pk Ok Ms (5, aa T  On+1° (5) 
Andererseits liefert (2) nach Division durch c, und Differentiation nach der Zeit 
Te My . (6) 
Aus (5) und (6) folgt als Differential-Differenzgleichung des Drehschwingungsproblems: 
Mr-1 M, 1 1 Mk+1 
we M a 
O, ck (a, Gena) 4 Or+1 ee (7) 


Da fiir die Abhangigkeit von der Zeit nur trigonometrische Funktionen in Betracht kommen, 
d.h. 


M, SS M, sin (@ t) (8) 


1 Biezeno-Grammel, Technische Dynamik Bd. 2, Kap. X, 2. Aufl. Berlin, Gottingen, Heidelberg 1953. 
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(w Eigenfrequenz), so gilt 
ih, = 
und (7) geht in die Differenzengleichung 


Mr—-1 a2 1 1 
O, les oO, ou aie 


mit den Anfangs- bzw. Endbedingungen M, = M, = 0 iiber?. 


— w? M;,, , 


(9) 


Mr+1 


Or+1 


(10) 


3. Die Frequenzdeterminante. Fiir die n—1 Unbekannten M, gelten daher n—1 lineare 
homogene Gleichungen, welche nur dann auf nicht-triviale Lésungen fiihren, wenn die Koeffi- 
zientendeterminante verschwindet, d.h. es muB 


Dee i ; i 
Cy 0; O, 0, : y 
il w” il 1 1 
== (() 
1 w? ] 1 1 
0 — 0 0 
0, c¢, O, O, O, (11) 
1 w? 1 1 
0 
: “ c Over en=1 On-1 Go, 


sein. Wird noch das Vorzeichen jeder geradzahligen Spalte und jeder ungeradzahligen Zeile 
gewechselt, so stehen vor allen Tragheitsmomenten positive Vorzeichen: 


1 1 wo? | il 
tea 3 lal: oeeo+ oe ; 
1 1 1 w? | 1 
aie eae ge er — 0 0 0 
0, 2 ‘ 3 Co Os; 
1 1 1 w? 1 
0 moll eee MiB ee ae — 0 0 (12) 
0; J O, C3 O, 
it 1 i w? 
0 0 0 o ne Onat. Oni a Qa Ga 


Die gestrichelten Abgrenzungen zeigen die fiir n = 2, 3, 4 usw. jeweils in Betracht kommenden 
Determinanten. 


4, Die Auflésung der Determinante. Bei der Auflésung der Determinante ist es vorteilhaft, 
zunachst die fiir w = 0 verbleibende Restdeterminante naher zu untersuchen. Diese wird 


- ape na | e emnes 
fine yi— sar D, 6, be me Oe: 
1 I 1 
% (5, + ) O, eae 0, + 9, + 9; 
ee: Dace) Shy We re 00,0,” 
0, G a) 
1 J 1 
Se ls —— 0 
rts 
1 ene 1 _64+0,4+0,+ 6, ; 
fir n=A: Die" oo; ogi) 6, Fo ah eee. 
I 1 if 
° 0; fae a) 


Ei Die Differenzgleichung (10) ist identisch mit der Grammelschen Rekursionsformel; siehe Biezeno- 
Grammel, Technische Dynamik, Bd. 2, S. 362—363. 


18 
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demnach allgemein bei n Massen: 


Dp = (14) 


k=1 
Fir die Hauptunterdeterminanten gilt Entsprechendes. Sie bestehen aus einem Pro- 
dukt von zwei Determinanten, fiir welche das Bildungsgesetz (14) maBgebend ist, und er- 


rechnen-sich aus 


=a (Sell 9.) (15) 


l=1 m=k+1 


oder mit Bezug auf (14) ! 
H, = D,( 3 @,) ( s @,) Vie @,): (16) 


Wegen der Identitat 


. 


‘SG a ee. (17)| 


1 1 
Mh, = Dy, jh i 5 
| sO, ~ 2 5H (18) 
p= =1 | 
geschrieben werden. 


Weiter treten bei Auflésung der Determinante (12) auch Teildeterminaten auf, welche? 
ebenfalls dem Bildungsgesetz (14) entsprechen, sich jedoch nur von einer bestimmten (h-ten) 
Zeile und Spalte bis zu einer anderen (k-ten) Zeile und Spalte erstrecken. Zu ihrer Kenn-; 
zeichnung mége das Symbol 


kann (16) auch in der Form 


oodles aah ya lea eee | 

dh, == oO; ee mit ee pes 1 (19) 
dienen. Vor allem treten GréBen auf, welche Produkte dieser Teildeterminanten darstellen 
jedoch mit der Ausnahme, da bei benachbarten Teildeterminanten das Produkt durch die¢ 
Teildeterminante des so entstehenden gesamten Teil gebietes zu ersetzen ist. Zur abgekiirzten 
Schreibweise seien folgende Symbole eingefiihrt : 


Sa = hoi» (264) 
d d ae i} 
fan bt) fir k> ae (21) 
dines fir kesh 1) i) 
dy n+1 Ue, n+ U,141 fiir k> h+1;1> k+1, 
dh ne 141 fiir k=h+1; LSS h+2, ] 
Sari = i (22) 
dh n41 Ak, +2 fiir k> h+ 1; [pe k+1, 4) 
Upes fir k=htl1;l=h+2, | 
dh, h+1 Ak, k41 Opty i aaa fiir k> h+1, tes k+1, m>I+1, 
di, 5.2 Gat Gig ey 4 fiir k= hele I> k+1, m> I+ 1, 
dh naa Ge, p02 Gnome 1 fiir kh, t=k-+ 1, ms Tsoi 
Sakim = ) %,n41 Uk, e+1 Uh, 142 fiir k>htl,is>k1, m=]=1-Pr (235) 
deena doen: fir k=A+1; 1=h+2, m>14+1, | 
din+1 Ak, n+3 fiir k> hl lS kh me ee 
dhn+4 fiir k=h+) l= h-2 0m has 
usw. 
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Mit Hilfe dieser Abkiirzungen fiihrt die Aufl lésung der Determinante (12) auf folgende 
Gleichung (nm — 1)-ten Grades fiir w?: 


n—l 
1 n—l 
eet) LT 1? 2 — gy 2-4 3 Cf, + @? *™-$ > Sy Hee sin 
p=l “p h=1k=h+1 


n—3 n—2 — 
= orn = Ds x Ch Ck C1 fiat (24) 
hal) Roh 1h 
hy 


Sigg? est sy "5 "Ss tho em farm — + °*-|-+ Dy = 0. 


h—1 k=h+1 l=k+1 m=I141 


Dividiert man diese Gleichung durch w?"—*, so folgt 


D si SF 5 ee fe ve 
ze 1 1 . Sh ist 1 “ . Shk 1 Hess + = Ehkl 
@M2n—2 @2n—4 OS Cc, w2 n—6 ps %  @pDst ae = me c 


Cc, Cc Cc, ¢ 
b= Hl knee (=1 k=hs Lee 


m4; na 3 fo al + a 
Sees > > ae ee be] Caps pio, 


Cc; ¢ 
i Phat teow Dimi Re hove 


Hierbei wurden noch folgende weiteren Abkiirzungen verwendet: 
h n n 
-(36)( 3 o}f(39) 
p=1 q=h+1 re 
h k n n 
sx = ( 3 0,)( va 6,)( a 6,) (36), (26) 
p=1 s 


q=h+l1 r=k+1 =] 
h k l n n 
ou =(2 6) 2 0) 2 6)( 2 o)f/(Se) 


usw. Zu demselben Ergebnis ware man auch durch Betrachtung der Potenzreihe fiir den 
Kehrwert des Fiequenzquadrates gelangt. 


5. Sdtze iiber die Eigenfrequenzen. Wird weiter beachtet, da die Determinante (12) 
infolge der Symmetrie reelle Wurzelquadrate hat und da die Potenzentwicklungen (24) 
und (25) alternierende Vorzeichen aufweisen, so folgt die bekannte Tatsache, daB bei n Massen 
sdmtliche n — 1 Eigenfrequenzen reel] sind. Schreibt man entsprechend dem Fundamental- 
satz der Algebra die Potenzentwicklungen (24) und (25) in Produktform, so gilt, wenn die 
Eigenfrequenzen mit w, (Grundfrequenz) <@, <@3 <::+: <@n—», allgemein mit @, ge- 
-kennzeichnet werden, 


(co? — w?) (w? — w2) se (w? — w?_3) = 0 (27) 


iaSet AY Lak aed 1 aes 
Oo Racoon 
Durch Koeffizientenvergleich erhalt man aus (24) und (27) 


bzw. 


n—l n—l n—2 n—l n—2 n—l 


Pa (he rips ba) pes, Ono => > Ch Chlbns 
h=1 h=1 h=1 k=h+41 h=1 k=h+1 29 
n— n—3 n— n—l ( ) 
"3 = ox Wh Ok OL = 2 = oS Ch Ck er fk 
mail eel fieey een h=1 k=h41 l=k+1 
usw., bzw. aus (25) und (28): 
utes 1 Begs n—2 n—l 1 be Sy eS Lvs 
= Oe ee Ce > so 2 m2 Cuiey 
esi to ie FDS (Oe h=1k=h+1 ah 
n—3 n—2 —!l n—3 n—l ( ) 


For SS ae 
ye Oh Pk I Ch °k | 


h=lk=h+1 I=k+ =1k=h+1 1=k+1 
usw. 
Es stellt sich heraus, daB die in (29) und (30) rechtsstehenden Ausdriicke Frequenzprodukte 
bzw. deren Kehrwerte darstellen, und zwar von Ersatzsystemen, die einerseits durch Weg- 
Leis 
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nahme von Steifigkeiten andererseits durch Erstarrung einzelner Wellenteile entstehen. Zu- 
nichst erkennt man in der ersten Gleichung (29) als rechtsstehenden Summanden die Figen- 
frequenz eines Zweimassensystems, gebildet von den Tragheitsmomenten Or, und Oni1> 
mit dem dazwischen befindlichen Wellenteil (Federkonstante c,, Abb. 2). Wird das zuge- 
horige Frequenzquadrat mit a3,, gekennzeichnet, so gilt mit Bezug auf die Definitionen 


(19) und (20) 


Se were pes teat Ain = Orth = Ch hry (31) 
| | | %h | | | In der zweiten Gleichung (29) treten gemaB 
73 Ons7 der Vorschrift (21) einerseits Produkte solcher 

Abb. 2. Zweimassensystem, entstanden durch Wegnahme Frequenzen von Zweimassensystemen auf, 
aller Steifigkeiten, bis auf die Stelle h. andererseits aber auch (bei benachbarter Lage 

der Stellen h und k) Produkte von Frequenz-: 

aot | | : | | | |  quadraten von Dreimassensystemen, wie sie: 
ee eos ieee al) Abb. 3 zeigt. Die nahere Betrachtung zeigt | 
i zF 4 1 9 ‘sofort im AnschluB an die Definition (19), da8B} 
O,nO7 Oso in der Tat fiir das Produkt der Quadrate der: 


Abb. 3. Dreimassensystem, entstanden durch Wegnahme aller beiden EKigenfre quenzen des in Abb. 3 dar- 
Steifigkeiten, bis auf die benachbarten Stellen h, h+1, h+2. gestellten Dreimassensystems die Beziehung gilt; 


02, hO3,h41 = Ch Ch4i Sh hoa = Ch Ch41 di,nz2+ (32)) 


Hae Ch__| Cher) Che eh Ebenso stellt der Summand der rechten Seite 
der dritten Gleichung (29) das Produkt dreier’ 

‘ Led | , Ones Oper ngs | me Fre quenzquadrate dar, wobei es sich ent-- 
: sprechend der Vorschrift (22) entweder um die 
Steifigkeiten bie wel die benachbarten Sellen te het hex, Eigenfrequenzquadrate von drei selbstindigen} 
ite Zweimassensystemen oder um ein Dreimassen-- 

und ein Zweimassensystem, oder um ein Vier-- 

massensystem (Abb.4) handelt. Die jeweils auftretenden Kigenfrequenzen erhalten als zweiten) 
Index einen der Indices der in Betracht kommenden Federkonstanten. Bei benachbarter Lage > 
gilt fiir das Produkt der Frequenzquadrate des Viermassensystems 


2 2 2 = — 
a3, h X3, h+l X3, h+2 rs Ch Chi4 Chi oth he hes a Ch Chat. Chie ding . (33)) 


Allgemein gilt fiir das Produkt der Frequenzquadrate eines aus q + 1 Massen, also q Wellen-; 
abschnitten bestehenden Teilsystems die Beziehung 


q x 
2 == } 
if QXg, ht+p—1 = IT Ch+p—l di hate : (34) 
1 1 


Analog lassen sich die Summanden der rechten Seiten von (30) deuten. In der ersten Glei- 
chung (30) erkennt man mit Riicksicht auf die Definition (26) als Summanden den Kehrwert, 
des Frequenzquadrates eines durch Erstarrung, 

Ch aller Wellenteile mit Ausnahme der Stelle 

entstandenen Ersatzsystems (Abb.5). Alslinke 

Masse gilt die Summe aller Massen von der 


0, 6) G; O, 4 : ; 
7 é h +7 Gn-1 Stelle 1 bis zur Stelle h, als rechte Masse die 
Abb. 5. Zweimassensystem, entstanden durch Erstarrung aller : ai 
Willcaréilo iat? Ace dale sdee Stace Summe aller Massen von der Stelle h+-1 bis zur 


Stelle n, wobei alle Federkonstanten auBer cy) 
unendlich sind. Im Hinblick auf die erste Definitionsgleichung (26) erkennt man, da in de 
Tat fiir das Frequenzquadrat dieses Systems die Beziehung 


y= ¢,°°° Ch—1 Ch41 Chang == **' = hy = OW, 
ee 35 
Th gh ( ) 


gilt. In der zweiten Gleichung (30) handelt es sich um Produkte der Kehrwerte der Fre quenz” 
quadrate von Dreimassensystemen, wie sie durch Erstarrung aller Wellenteile mit Ausnahme: 
von zweien entstehen (Abb. 6). Der Vergleich des Bildungsgesetzes der zweiten Gleichung (26) 
mit den GesetzmaBigkeiten beim Dreimassensystem bestatigt unmittelbar diesen Zusammen-: 
hang. Fiir das Produkt der Frequenzquadrate des in Abb. 6 gezeichneten Ersatzsystems er-} 
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halt man daher 


2 3 ae ch Sle 
B3,n 2, x bh, k : (36) 


In der gleichen Weise ]aBt sich der in der dritten Gleichung (30) rechts auftretende Summand 
als Kehrwert des Produktes der drei Eigenfrequenzquadrate eines durch teilweise Erstarrun 
nach Art von Abb. 7 entstehenden Viermassensystems deuten. Man erhilt fiir das Produkt Ts 
drei Frequenzquadrate des in Abb. 7 ersicht- 

lichen Systems mit Bezug auf die dritte C3 

Gleichung (26) den Ausdruck 4 - 


Cy Cy, C 
2. pe ee MAGES () 
B3,n° B3,4° B3,1 = semi On 7 Gres Oe Sher ®n 
hkl a 
Abb. 6. Dreimassensystem, entstanden durch Erstarrung aller 
Wellenteile, mit Ausnahme der Stellen h und k. 


usw. 

Die weitere Fortsetzung ergibt sich entsprechend dem leicht erkennbaren Bildungsgesetz. 
Die so gewonnenen Deutungen bilden den Inhalt wichtiger Eigenschaften der Eigen- 

frequenzquadrate, die sich in besonderen Satzen zusammenfassen lassen. Die erste Gleichung 


(29) fiihrt im Hinblick auf (31) zu folgendem 
O Gj Ch 


G, On Oat Ox Ore1 OQ: Ons On 


Abb. 7. Viermassensysteme, entstanden durch Erstarrung aller Wellenteile, mit Ausnahme der Stellen h, k ni is 


Satz 1: Die Summe der Eigenfrequenzquadrate ist gleich der Summe der 
Frequenzquadrate aller Zweimassensysteme, die sich durch Wegnahme aller 
Steifigkeiten bis auf eine herstellen lassen. 

Ebenso fiihrt die erste Gleichung (30) im Zusammenhang mit (35) auf folgenden 

Satz 2: Die Summe der Kehrwerte der Frequenzquadrate ist gleich der 
Summe der Kehrwerte der Quadrate der Frequenzen aller Zweimassensysteme, 
die durch Erstarrung aller Wellenteile auBer einem herstellbar sind. 

Die weiteren Gleichungen (29) und (30) enthalten ahnliche Aussagen tiber die Summen der 
Produkte der Quadrate der Eigenfrequenzen bzw. deren Kehrwerte, die hier jedoch nicht ein- 
zeln aufgefiihrt werden sollen. 


6. Aufstellung oberer und unterer Schranken fiir die Eigenfrequenzen. Nachdem nunmehr 
die Koeffizienten der Gleichung (n—1)-ten Grades in einfacher Weise als Produkte von Fre- 
quenzquadraten bzw. deren Kehrwerte gedeutet werden kénnen, ist es leicht méglich, obere 
und untere Grenzwerte fiir die Eigenfrequenzen anzugeben. Zunachst liefert die erste Glei- 
chung (29) einen oberen Grenzwert fiir die héchste Eigenfrequenz 


n—l 
On he ak b On th (38) 
Kol 


Ebenso liefert die erste Gleichung (30) einen unteren Grenzwert fiir die niedrigste Eigenfre- 
quenz 

1 
<j Pi,n 


Weitere Grenzwerte erhalt man unter Beachtung der bekannten Gesetzmafigkeit, daB die 
Eigenfrequenzen bei Erhéhung der Steifigkeiten zunehmen, bei Erhéhung der Massen jedoch 
abnehmen. Die sinngemaBe Anwendung auf n-Massen Systeme fiihrt zu folgendem 

Satz 3: Die n—1 Eigenfrequenzen eines n-Massensystems nehmen samtlich 
bei Erhéhung einer Federkonstanten des Systems zu und bei Erhéhung einer 
Masse des Systems ab, und umgekehrt. 

Nimmt man an einem System Steifigkeitserhéhungen oder Massenverringerungen belie- 
biger Art vor, so erhalt man fiir das neue System neue Eigenfrequenzen, die wieder der Gré8e 


| 
| 
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nach zu ordnen sind. Fiir die Frequenzen des Ausgangssystems galten die Ungleichungen 
Dy < Wz <.W3 <**  On_ - (40). 


Werden die Frequenzen des neuen Systems durch einen Querstrich gekennzeichnet und in der-: 
selben Weise der GréBe nach geordnet, so gilt 


Or Wee (41), 
Der Inhalt des vorstehend ausgesprochenen Satzes besagt dann fiir Streifigkeitserhéhungen | 
bzw. Massenverringerungen folgende Ungleichungen 


OQ, >; ®,>@,, O,>03, **' O,_1 > @n_1 + (42) | 


Der Satz ist offenbar in dieser Allgemeinheit noch nicht bewiesen worden, so daB® der mathe-» 
matische Beweis nachgeholt werden mu8. Fiir ein Zweimassensystem mit den Massen 0,, 0), 


4 Cy bp 


0, OQ; QO, Op O; 
Abb. 8. Zweimassensystem an der Stelle 1 mit den Abb.9. Dreimassensystem, gebildet von den Wellen- 
unmittelbar benachbarten Schwungmassen. teilen 1 und 2 mit den Schwungmassen 0, 0,, O3. 


und der Federkonstanten c, (Abb. 8) geht dieser Sachverhalt direkt aus der Forme] fiir das} 


Frequenzquadrat hervor: 


1 1 
On ala + a): (43) } 
Fir cin Dreimassensystem mit den Massen 0,, @,, @, und den Federkonstanten ¢,, ¢, (Abb. 9)) 
gelten mit Bezug auf (29) folgende Beziehungen: | 
2 ge (eels liked ) 
Oy ae OU ie ar é;) 1 VP) (5, ir 0,)” (44) 
On Ox Os 


CONS Sh Bre. o. (45) 


Man erkennt hierbei zwar direkt aus dem Aufbau der rechten Seiten dieser beiden Gleichungen, 
daB sowohl bei Zunahme einer Federkonstanten, als auch bei Abnahme einer der drei Massen | 
die Summe und das Produkt der beiden Frequenzquadrate zunehmen. Daraus folgt jedoch | 
nicht, daS jedes Frequenzquadrat einzeln zunimmt. Zur Durchfiihrung des Beweises ist es 
zweckmafig, zunichst aus den Gleichungen (44) und (45) einige Ungleichungen herzuleiten. | 
Quadriert man (44) und subtrahiert die mit 4 multiplizierte Gleichung (45), so folgt 


EM AN iio die ip PS \e Over O54, 1 

(2 — w})? = E (5, ar 6, + C (5, ate al — 4 ¢ 0, O, O, (46) 

bzw. nach Umformung 

pe ognee lee See a6 
(w3 4) = E (3; i é;) Cy (5; ae all “{ tore (47) | 

Demnach besteht die Ungleichung | 
ee tien a | 

wt — 0? > aa(a, é;] 1 (5, a, (48) 
Mit Riicksicht auf (44) folgt entweder 1 
1 gabe er | 
a> alo, +5)> ala +a)> ot ae 


oder | 
Lgauel Lante | 

o8> alg +g) > ala, + @)> oF ee 
Wird nochmals (45) herangezogen, indem das Produkt der Fre quenzquadrate durch c, é - a 


bzw. aa(g, + — 7) Z) dividiert wird, so folgen im AnschluB an die Ungleichungen (49) und (50) nau | 
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die weiteren Ungleichungen 


. 1 ie fis 1 
si = ( 0, — 0, g,| es : G 0, 1 g,)>ot a) 
oder 
; 1 1 1 1 
We = (5, 0, ac 0; = Cy (5, 0, de 0, > wt % (52) 


Zu demselben Ergebnis ware man gelangt, wenn man von den Kehrwerten der Fre quenz- 
quadrate ausgegangen ware und (30) zugrunde gelegt hatte. 

SchlieBlich wird noch eine weitere Umformung der Gleichung (47) bendtigt. Man erkennt 
leicht, da® folgende Darstellung mit (47) identisch ist: 


1 il il cr (c, + ¢,)? c2\]2 , 4c? c2[e, O ) 
©? — wo)? = biog: | ef 1 2 2 | 1&2 |¢; 3 Co Oy \ 
( 2 1) (cy + c,)? | coe @, ' 05 0, 0, a 0; O, Os |e, O; r c, Os ; iM 
3 


Andererseits ]a8t sich (44) auch folgendermafen darstellen: 


2 oa 1 PT ee ce (Cy =| €5)7 ca 
CIA hes ape ue (5, . each a a alk C 


Aus (53) geht hervor, da die Ungleichung 


il 1 c (c, + ¢2) c2 
| a | 1 2 2 
Cy C5 (5, T 5; (5 T 0, mi 5) 


besteht. Mit Bezug auf (54) ergeben sich daher als weitere Ungleichungen fiir die Frequenz- 
quadrate 


(55) 


: gen fl’ 1 ca (Gy rea)* yy e8 2 
pg Cy + Cy (5, a) z Cy + | 0, r 0, 0; cae! 28) 
d 
odaer A i cr ' (¢ — 3) | c Cy Co a l 1 2 5 
Ona. ei ak 0, ; al = Cy + Cy (5, 5;) sae ; ( > 


Wird wieder (45) herangezogen und das Produkt der Frequenzquadrate abwechselnd durch 
eine der beiden, in den Ungleichungen (56) und (57) auftretenden Schranken dividiert, so fol- 
gen die weiteren Ungleichungen 


0,+@ ) 1 1 

ee ae oe 2 > (e+ ¢) (5; a45,)>% (58) 

(0) a) (0) 1 ate 1 2 2 2 1 3 

: Se O; OQ, O; 
oder 

a 1 i C1 Cy (Cc, +°¢,) (9, + O, + O5) 
w8 > (44+ 4) (5 ee ay ot (59) 

2 1 3 0, 6, 6,| 1 { 1 2 i 

0, 0, 0, 


Auch diese Ungleichungen hatte man — ebenso wie die Ungleichungen (51) und (52) — aus 
den fiir die Kehrwerte der Frequenzquadrate geltenden Gleichungen (30) herleiten kénnen. 
Damit sind eine Reihe unterer Schranken fiir die Oberfrequenz und oberer Schranken fiir die 


== [ta 
GO, @, A, 6, 0, 0, 


Abb. 10. Ubergang vom Dreimassensystem der Abb. 9 in Zweimassensysteme, 
durch Erstarrung jeweils eines der beiden Wellenteile. 


Grundfrequenz gewonnen. Sie lassen sich z. T. als Frequenzen von Zweimassensystemen deu- 
ten, welche durch Wegnahme von Massen, Verschwinden von Federkonstanten oder teilweise 
Erstarrung aus dem Dreimassensystem entstehen. So sind in den Ungleichungen (49) und (50) 
auftretenden Schranken die Frequenzquadrate von Zweimassensystemen, die durch Verschwin- 
den einer der beiden Federkonstanten oder einer der beiden Endmassen aus dem Dreimassen- 
system entstanden sind. Bei den in den Ungleichungen (51) und (52) vorkommenden Schran- 
ken handelt es sich um Frequenzquadrate von Zweimassensystemen, die durch teilweise Er- 
starrung des Dreimassensystems entstanden sind (Abb. 10). Andererseits entspricht der in 
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der Ungleichung (56) links, bzw. in (57) rechts stehende Schrankenwert dem Frequenzquadrat} 
des bei Wegnahme der mittleren Masse entstehenden Zweimassensystems (Abb. 11). SchlieB-- 
lich findet sich in der Ungleichung (58) rechts, bzw. in (59) links als Schrankenwert das Fre-- 
quenzquadrat eines Zweimassensystems vor, das durch starre Verbindung der duferen Massen) 
aus dem Dreimassensystem entstanden ist (Abb. 12). 


Cy Ce 
Q; @ 
Abb. 11. Ubergang vom Dreimassen- Abb.12. Ubergang vom Dreimassensystem der 
system der Abb. 9 in ein Zweimassen- Abb. 9 in ein Zweimassensystem durch starre 
system durch Wegnahme der mittleren Verbindung der beiden auBeren Schwungmassen 
Schwungmassen, miteinander. 


Nach diesen Vorbereitungen gehen wir zuriick zu (44) und erhalten durch Bildung totaler: 
Differentiale 


1 1 1 1 e Cpe iiCs Co 
Hot) + Hod = (5 + g,) 4 + (6, I g,) dee #40, — 17-40, — G40. (60)) 


Andererseits wollen wir (45) einer logarithmischen Differentiation unterziehen; dann folgt 


d(w}) , d(w3) de, , de, 0, + @; d O, + O; dO 
wo? wo Cy =e Cy 0,(0, + O, + O;) 1 0, (0, + @,+ 0,) ~? 
al ES. (61)) 


3 (0, + O, + ) 
Durch Auflésung nach den Differentialen der beiden Frequenzquadrate ergibt sich hieraus 
maby (leo eV alee eh, ty 
a() oe O, a 0.) w} fai Co 0, a O,] w3 dey 


@, + 03 Cy Orta OG; Cy + Ce dO 
Lave, + 6, + 6,) ai| ioe xc +6,+6,) OF |°? 


oa pais a - gl 400l 7 

od apes lla Fe)ag~ a] [(a, ba)ar—e] | 
lara 0, ee ay dO, — Ex a Soa : aa cn a5 | 40. | 

os ea A ae a,y| 4s}: (3) 


Wie man an Hand der Ungleichungen (49), (50), (51), (52), (58) und (59) leicht einsieht, sind die | 
Koeffizienten der Differentiale der Federkonstanten stets positiv, jedoch jene der Differentiale } 
der Schwungmassen stets negativ. Damit ist bewiesen, da8 die Erhéhung einer der beiden. 
Federkonstanten oder eine Abminderung irgend einer der drei Schwungmassen jede der beiden | 
Kigenfrequenzen erhéht (durch zahlenmafBige Berechnung der Koeffizienten kann man sich. 
auferdem einen Uberblick iiber die Wirkung einzelner Steifigkeits- oder Massenanderungen 
auf die Eigenfrequenzen des Systems verschaffen). Damit ist der als Satz 3 bezeichnete Sach- 
verhalt auch fiir das Dreimassensystem nachgewiesen. In ahnlicher Weise ]a8t sich der Beweis | 
auch fiir Systeme mit noch héheren Massenzahlen durchfiihren. Der Verfasser beabsichtigt, , 
noch eine alJgemeinere Beweisfithrung im Anschlu8 an algebraische Hilfssatze an anderer’ 
Stelle nachzuholen. 

Nachdem nunmehr die Giiltigkeit des Satzes 3 fiir Systeme mit beliebig vielen Massen vor- 
ausgesetzt werden darf, ist es ohne weiteres méglich, GesetzmaBigkeiten fiir die Aufstellung ; 
oberer und unterer Grenzwerte der samtlichen Eigenfrequenzen anzugeben. Die nachstehenden | 
beiden Satze ergeben sich sinngem4® aus Satz 3, wenn die Steifigkeitserhéhung einerseits bis | 
zur Erstarrung, andererseits die Steifigkeitsabminderung bis zur Auflésung der elastischen | 
Bindung gesteigert wird: 
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Satz 4: Durch Erstarrung vont Wellenabschnitten entsteht ein neues Dreh- 
schwingurgssystem, fiir welches jede seiner n—t—l-Eigenfre quenzen (der Grobe 
nach geordnet) gréferist als die entsprechende Eigenfrequenz des Ausgangs- 
systems, d.h.esgilt, wenn die neuen Frequenzen wieder durch einen Querstrich 
gekennzeichnet werden, 


Oa Wis Oe W538i Ogg SS On aes On_441 = Dy—ire oy On oO. (64) 


Satz 5: Durch Wegnahme von q benachbarten oder nicht benachbarten Stei- 
figkeiten reduziert sich bzw. zerfallt das Ausgangssystem in zwei oder mehr 
Teilsysteme. Werden samtliche neuen Frequenzen wieder der Gréfe nach ge- 
ordnet, so gilt 


O, = 0 == O, = 0; Oy 41< Oo433 Daya Oyen? Op Ont (65) 


Mit Anwendung von Satz 4 ist der Frequenzrechner in der Lage, sich durch Vergleichs- 
rechnungen an durch teilweise Erstarrung gebildeten Ersatzsystemen obere Grenzwerte fiir 
die Grundfrequenz und die nachst héheren Frequenzen zu verschaffen. Die hierzu erforder- 
liche Rechnung ist infolge der Verstiimmelung des Systems wesentlich kiirzer als die des Haupt- 
systems. Andererseits liefern Vergleichsrechnungen an Ersatzsystemen, die durch Wegnahme 
einzelner Wellenteile entstanden sind, untere Grenzwerte fiir die héchste und die nachst nied - 
rigeren Frequenzen, wie in Satz 5 dargelegt ist. Da andererseits in (38) und (39) ein unterer 
Grenzwert fiir die Grundfrequenz und ein oberer Grenzwert fiir die héchste Frequenz gegeben 
waren, lassen sich Grundfrequenz und héchste Frequenz nunmehr bereits eingrenzen. Sind 
beide Frequenzen mit ausreichender Genauigkeit ermittelt, so lassen sich die zweite und die 
zweithéchste Frequenz in derse]ben Weise eingrenzen. An die Stelle der Gleichungen (38) und 
(39) treten dann die fo] genden, welche unmitte)] bar aus (29) und (30) abgelesen werden kénnen: 


wr? = wh?, — wi >wi_s, (66) 
1 i 1 1 

; 67 

or? w¥? Bie ( ) 


Damit sind ein unterer Grenzwert fiir die zweite und ein oberer Grenzwert fiir die zweithéchste 
Frequenz bekannt. Die Eingrenzung nach der anderen Seite ergibt sich wieder aus den Satzen 4 


und 5. 
Nach dieser Methode ergeben sich eine Reihe guter Naherungswerte, sowie wichtige Kon- 


trollen fiir die numerische Rechnung. 


7. Zusammenfassung. Durch Verwendung der elastischen Momente anstelle der Schwin- 
gungsausschlage als Unbekannte ergibt sich eine (n—1)-reihige Determinante, deren Auf- 
lésung auf eine Gleichung (n—1)-ten Grades fiir das Frequenzquadrat fiihrt. Die Untersu- 
chung der Koeffizienten dieser Gleichung fiihrt zur Erkenntnis besonderer Gesetzmafigkeiten 
der Eigenfrequenzen, aus denen einfache Verfahren zur Berechnung oberer und unterer Grenz- 
werte der einzelnen Frequenzen gewonnen werden. 


(Eingegangen am 24. Oktober 1953.) 
Anschrift des Verfassers: Professor Dr.-Ing. H. Neuber, Dresden-A 27, Westendstr. 23. 
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Turbulente Strémungen in divergenten Kanilen. 
(Mittlerer und starker Druckanstieg.) 


Von W. Szablewski. 


1. Einleitung. Das Verfahren von v. Kdérmdén und Pohlhausen zur Berechnung laminarer | 
Grenzschichtstrémungen arbeitet bekanntlich mit einem Polynomansatz, der die Wandbin- | 
dungen erfiillt und dem noch dauSere Randbedingungen auferlegt werden. Ein analoges Ver- | 
fahren fiir turbulente Grenzschichten ist bisher nicht entwickelt worden. Ein Arbeiten mit | 
einem Polynomansatz ist, wie man heute wei}, hier nicht méglich, da der analytische Charakter | 
der Geschwindigkeitsverteilung turbulenter Grenzschichten anders beschaffen ist (man denke | 
z. B. an das logarithmische Gesetz fiir die Geschwindigkeitsverteilung in Wandniahe). 

Es hat sich jedoch in den bisherigen, auf der Grundlage der Prandtlschen Mischungsweg- | 
hypothese basierenden Untersuchungen des Verfassers! iiber einige wichtige Modelle der 
Mechanik turbulenter Strémungen gezeigt, daB man auch hier zu brauchbaren Funktionsan- | 
sitzen fir die Approximation der Geschwindigkeitsverteilung gelangen kann. Fiir turbulente | 
Grenzschichten mit starkem oder mittlerem Druckabfall erweist sich, wie die Analyse der | 
konvergenten Kanalstroémungen ergab, das logarithmische Geschwindigkeitsgesetz als ein sehr | 
brauchbares Mittel der Approximation. Auch bei schwachem Druckabfall, wie z. B. bei Rohr- | 
strémungen, leistet dieses Gesetz bekanntlich gute Dienste 2. Man erzielt jedoch beischwachem | 
Druckabfall der Strémung bessere Ergebnisse, wenn man das in einem ersten Integrations- | 
schritt von der Wand her aus der Annahme konstanter Schubspannung zu gewinnende log- | 
arithmische Gesetz in einem zweiten Integrationsschritt erweitert. Das gab — vgl. die Unter- | 
suchungen des Verfassers? ttber die turbulente Strémung im Rohr und langs der ebenen | 
Platte — ein zusatzliches, in der Variablen des Wandabstandes lineares Glied. | 

Die vorliegende Untersuchung der turbulenten Strémungen in divergenten Kanalen, die | 
vom Typ turbulenter Grenzschichtstrémungen mit Druckanstieg sind, wird erweisen, dab 
man auch in diesem technisch besonders wichtigen und interessierenden Fall auf der Grund- | 
lage der Prandtlschen Hypothese durch analytische Entwicklung von der Wand her zu Funk- | 
tionsansatzen gelangen kann, die eine gute Approximation der Geschwindigkeitsverteilung | 
uber der Grenzschichtbreite erméglichen. 


2. Das Wandgesetz turbulenter Grenzschichtstr6mungen mit Druckanstieg. Im folgenden 
gebrauchen wir die iblichen Bezeichnungen: x Bogenlange der Kérperkontur, y Wandabstand, 
6 Grenzschichtdicke, u, v Geschwindigkeitskomponenten in x, y-Richtung, 0 Dichte, pu Zahig- | 
keit, y = y/o kinematische Ziahigkeit, p (x) statischer Druck, t Schubspannung. | 

Wir erinnern zunachst an die Herleitung des bekannten logarithmischen Geschwindigkeits- _ 
gesetzes *. Vernachlassigt man fiir eine wandnahe Zone der Grenzschicht die Tragheitskrafte 
gegentiber den Druck- und Reibungskraften, so reduziert sich die Bewegungsgleichung der | 
Grenzschichtstrémung auf 


dp , 0r — 
Segegnittay ie 0. 
Integration ergibt 
d 
tT =% + ay (t) Wandschubspannung) . (1) 


Hierbei ist 
UG = Tlam a Tturb 


mit 
Ou 
T —i | 
lam oy 
mie Z. angew. Math. Mech. 31 (1951), S. 131; ebenda 31 (1951), S. 309; Ing.-Arch. 20 | 
21, Prandtl, Ergebnisse der Aerodynamischen V: h talt Gotti 
Rat y ersuchsanstalt Gottingen, IV. Lfg. (1932), S. 18. 
4 W. Szablewski, Ing.-Arch. 20 (1952), S. 37. 
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und bei Zugrundelegung des Prandtlschen Ansatzes 


Tturh = 0 Ne (3) 


wo | den Mischungsweg bedeutet. Unter der Voraussetzung t) 0 kann man fiir hinreichend 
kleine y auch noch den Einflu8 des Druckgradienten vernachlissigen und erhilt 


TH% (tq 0). (2) 


Die Integration von (2) 1aBt sich, wenn man noch mit der Dicke 5, der laminaren Unter- 
schicht 1 = x (y —- 69) setzt, unschwer durchfihren und ergibt fiir ,,nicht zu Kleine“ y (mathe- 
matisch in der Asymptote) das logarithmische Geschwindigkeitsgesetz 

p =! ny + —1+In (4% Re,Ve) + x *% fiir ee ey ee (3) 
be v 2 Res\c 
Hierbei bedeuteny = y/6,0, = /t/o und, mit der Geschwindigkeit u, der Potentialstromung 
am Rande der Grenzschicht, p = u/u,, Re; = u, 6/v sowie cy = T/0 Tees 

Dieses Gesetz kann man bekanntlich auch in einfacher Weise direkt erhalten, wenn man in 
(2) die laminare Reibung gegeniiber der turbulenten Reibung vernachlassigt. Allerdings 
bleibt dann die Integrationskonstante unbestimmt. 

Die physikalische Existenz des logarithmischen Gesetzes ist experimentell fiir einen weiten 
Bereich turbulenter Grenzschichtstrémungen, angefangen mit Strémungen unter Druckabfall 
bis zu Strémungen mit Druckanstieg, einwandfrei sicher gestellt!. Fir turbulente Grenz- 
schichten mit starkem Druckanstieg und insbesondere in der Nahe der Ablésungsstelle wird 
man jedoch der analytischen Beschreibung der Geschwindigkeitsverteilung in Wandnahe die 
vollstandigere Gleichung (1) zugrunde legen miissen. Man ersieht das unmittelbar aus (1), wenn 
man beachtet, daB bei einer Grenzschichtstrémung gegen wachsenden Druckanstieg, der 
schlieBlich die Ablésung erzwingt, dem ansteigenden Verlauf des Druckgradienten ein monoton 
fallender Verlauf der Wandschubspannung Tt, bis zur Ablésungsstelle t, = 0 gegenibersteht. 

Die Integration von (1) mit 


0 du\2 
Tlam =f Tturb = Lb = a2 ? (55) 


laB®t sich nun leider nicht, im Gegensatz zu (2), mit elementaren Funktionen durchfihren und 
erfordert die Berechnung elliptischer Integrale. Enthalt die betrachtete wandnahe Zone jedoch 
eine Schicht, fiir die die laminare Reibung gegeniiber der turbulenten Reibung vernachlassigbar 
ist 2, so laBt sich fiir diese Schicht die Integration leicht durchfithren. Es bleibt dann allerdings, 
wie im analogen Fall des logarithmischen Gesetzes, die Integrationskonstante unbestimmt. 
Wir werden hier die Integrationskonstante approximativ aus einer einschrankenden Annahme 
ableiten. Ihre vollstandige Bestimmung behalten wir einer spateren Arbeit vor. 
Bei Vernachlassigung der laminaren Reibung ergibt (1) 


1\2 /dp\2 _ i _ = (dp/dx) 5 
hg eg ea a) i: 
Mit / = x7 kommt 
1 
ae * ey — (A/2) y 
gat | aera, (5) 
was auf 
Oo ears eR i ah 6 
= —JYe,—(A/2)y7 +14 In = =e const (6) 
eS amaccin 
fuhrt. 


Wir bestimmen nun die Integrationskonstante durch die Forderung, daB wir fir 7 —> 0 
wieder das log. Gesetz erhalten. Diese Forderung ist physikalisch so zu verstchen: Gehen wir 
bei einer Grenzschichtstrémung in den Bereich mit Druckanstieg hinein, so werden wir uns 
zunichst noch im Geltungsbereich des log. Gesetzes befinden, fiir eine wandnahe Schicht also 


1H. Ludwieg und W. Tillmann, Ing.-Arch. 17 (1949), S. 244. 
2 Die Existenz einer solchen Schicht ist nach den obigen Ausfiihrungen im Geltungsbereich des log- 


arithmischen Gesetzes experimentell nachgewiesen. 
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im Gleichgewicht der Krafte den Druckgradienten gegeniitber den Reibungskraften noch ver-- 
nachlassigen kénnen. Mit wachsendem Druckanstieg (und gleichzeitig fallender Wandschub-; 
spannung), greift nach (1) der Druckgradient immer weiter zur Wand hin bestimmend in das; 
Gleichgewicht der Krafte ein und engt die Schicht des log. Gesetzes immer mehr ein, um dieses; 
Gesetz schlieBlich aufzuheben. Die obige Forderung besagt dann, da8 wir unsere Betrachtung } 
auf nicht zu groBen Druckanstieg beschranken, bei dem das log. Gesetz noch nicht vollstandig ; 
aufgehoben ist. Fir starkeren Druckanstieg, insbesondere in der Nahe der Ablésungsstelle | 
(rt) = 0), werden wir diese Einschrankung nicht mehr treffen kénnen. 


Beachtet man, daB 


Vey —-(A/2) 4 = Vo ( = 2 Aj2 n) fiir kleine 7 , 


Cy 


so kommt : 
— 4/2 


Cy 


p volts 4 Fay 2! ) 7 + konst fir kleine 7 . 


Die Forderung, daf fiir kleine 7 wieder das log. Gesetz (3) erscheint, ergibt dann fir die Inte-) 
grationskonstante 


Vey ci il C4 eS + Oo 
konst = — 2 > ee Lee — 1 + In (4% Reg eq) + "4 : (7) | 


Damit lautet die von uns abgeleitete Form des Wandgesetzes 


2 — (4/2) n —Ve, 
p =2 (lq — OR) — Yq) + m4 4 a =o va 


C , 
+ ta | 1 + In (4 x Res Ve, ) + oad (9) 

fiir 
| | 

Deka len Wok os | 


In Prandtlschen Koordinaten 


eee und “#7 — Re, ey 


Vp Vey v 


lautet unser Ergebnis 


Ole: Oy | 


Ux ¥ 
Radda a yi ert Soares |) in4 
woo(ftret—-y+e ea 


2 Ura 
eee (9) 


xl 


mit 


1 (—A/2) __ v dp/dx 


ResVe, % x @% 

Empirische Koeffizienten dieses Gesetzes sind, nach dem derzeitigen Stand der Theorie, , 
die Koeffizienten x und v, do/v. Das sind die empirischen Koeffizienten des log. Gesetzes. Sie 
hangen von der Reynoldsschen Zahl des Strémungsmodells ab, sind aber vom Druckgradienten , 
der Strémung unabhangig. Fir gro®e Re-Zahlen nehmen sie den Charakter universeller | 
Konstanten an}, | 


: yu poe a 
Schicht > — < -* < Rese, . 


3. Verlauf der Wandschubspannungen in den Nikuradseschen Messungen. Zum Vergleich 
der Theorie mit dem Experiment ziehen wir die bekannten Messungen von Nikuradse® heran. 
Beziiglich der bei der Untersuchung der ebenen Diffusorstrémung verwendeten Symbole siehe 
Einleitung zu Abschnitt 4. | 

Hier standen wir zunachst vor einer gewissen Schwierigkeit. Die von Nikuradse angegebenen 
Wandschubspannungen weisen einen unregelmaBigen Verlauf auf, insbesondere zeigen sie auch | 
das seinerzeit viel diskutierte Anwachsen kurz vor dem Abfall auf Null. Dieses unwahrschein- | 


1 W. Szablewski, Ing.-Arch, 20 (1952), S. 37. | 


2 J. Nikuradse, Untersuchungen iiber di 6 i i 
Lansley VOCE en nae eos ie Strémung des Wassers in konvergenten und divergenten | 
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liche Verhalten der angegebenen Wandschubspannungen ist inzwischen als durch Fehler be- 
dingt erkannt worden *, Zur Erklatung hat man das Auftreten von Sekundirstrémungen, 
durch welche die Zweidimensionalitat der Strémung zerstért wird, herangezogen. Wir weisen 
hier auf folgende rein numerische Fehlerquelle hin: 

Nikuradse hat die Wandschubspannung 1, mittels des Impulssatzes 2 


1 
A 
eS tea | ot an 
0 


aus der gemessenen Geschwindigkeitsverteilung und dem gemessenen Druckgradienten be- 
rechnet. Numerisch heift das aber, wenigstens im Bereich starkeren Druckanstieges, als kleine 
Differenz groSer Zahlen, was bekanntlich immer eine miGliche Angelegenheit ist und zwanglos 
den streuenden Verlauf der angegebenen Wandschubspannungen erklart. 


In einfacher und sicherer Weise kann man nun die Wandschubspannungen aus der experi- 
mentellen Geschwindigkeitsverteilung in Wandnahe mittels unseres Wandgesetzes bestimmen, 
falls die empirischen Koeffizienten x, v, 6o/v desselben bekannt sind. Das ist hier der Fall! 
“x und v, d)/y sind, wie oben erwahnt, vom Druckgradienten der Strémung unabhangig und 
kénnen hier bei der mittleren Reynoldsschen Zahl Re = u, b/y ~ 81 500 der Messung der kon- 
vergenten Kanalstrémungen entnommen werden, die bei der nur wenig verschiedenen mittleren 
Zahl Re ~ 67500 durchgefihrt wurden. Dort hatten wir gefunden 


EKG: Ss: 


(Es ist zu bemerken, daB die Nikuradseschen Messungen bei einer verhaltnismafig kleinen 
Reynoldsschen Zahl durchgefiihrt wurden. Fiir grofBe Re-Zahlen, wie sie fir die technische 
Praxis besonders interessant sind, haben wir ungefahr ® 


er en LST ea 
Vy 


zu erwarten.) 


Unter Zugrundelegung des Wandgesetzes in der Form (10) haben wir bei Benutzung der von 
Nikuradse gemessenen Druckgradienten mit folgenden Werten der Wandschubspannung ein 
Zusammenfallen der gemessenen Geschwindigkeitsverteilung mit der Kurve des Wandgesetzes 
erreicht (Abb. 1): 


o Cy 
is 0,00137 
2° 0,00102 
3° 0,00075 
4° 0,00042 


Die so ermittelten Wandschubspannungen zeigen einen regelmaBigen, mit wachsendem Druck- 
anstieg monoton fallenden Verlauf (Abb. 2). Sie werden bei Zugrundelegung des vollstandigen 
Wandgesetzes noch etwas korrigiert werden miissen. 


Auffallend ist, iiber ein wie weites Intervall die experimentellen Geschwindigkeitsverteilungen 
durch das Wandgesetz brauchbar approximiert werden. Gestrichelt eingetragen ist das log- 
arithmische Gesetz. Man sieht sehr schén, wie die Schicht desselben mit wachsendem Druck- 
anstieg immer weiter eingeengt wird. 


4, Berechnung der ebenen Diffusorstrémung in erster Naherung. Man kann aus Abb. 1 
schon die Vermutung entnehmen, daB wir brauchbare Ergebnisse erhalten werden, wenn wir 
die Geschwindigkeitsverteilung tiber der gesamten Breite der Grenzschichtstrémung durch 
das Wandgesetz approximieren. a 

Es bedeuten im folgenden: b Kanalbreite, « halber Offnungswinkel des ebenen Diffusors, 
u, Geschwindigkeit in Kanalachse. 


1Vgl. H. Ludwieg und W. Tillmann, a. a. O. 
2 W. Szablewski, Ing.-Arch. 20 (1952), S. 37. 
3 W. Szablewski, Z. angew. Math. Mech. 31 (1951), 5. 309. 
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F ae __ — (dp/dx b) 
1G iat Trae TCP , 


lautet dann nach (8) das Wandgesetz 


a ns (A/4) 2 Vey 1) | Me ; Ha (A/4) )n—Vey 
x g (—A4) Ves — A) 0 + Vex 


Nias (10) 


De 


v 


Nach Abschn. 3 haben wir fiir Re = 81500 x = 0,6, “* Dimer 


Es verbleiben dann in der Funktion (10) noch die beiden unbestimmten Parameter c, und /. 
Zu ihrer Bestimmung ziehen wir den Impulssatz und die duBere Randbedingung py = 1 fiir 
7 = 1 heran. | 


= 
Wandschubspannung 


10 ore) 


Geschwindigke/tsverte/lung as aes 
i, “il 
ee 


Messung Nikuradse 


a Messung Nikuradse o Lmpulssarz 
ue Soe Tai os © Wanggese/z 
——— log Geserz cma hae ate! 
° oe 
erlg 7 
0 7 ? 


Abb. 1. Abb. 2. 


a) Randbedingung y=1 fir 4 =1. Die Randbedingung ergibt nach (10) 


1 =p (SE X yee In 4 | Cy Jc 1—A/4 —Ve, 
m % m = (aera 


As Ven E 1+ In(2% Repe,) + x = a (11) 
Wir kénnen dann fiir m auch schreiben 
pg =1 spy lees __ Vey Sew) ie (ss 4) 7 —Veq 'c, —A/4 + Ve. S (12) 
a cis a \ Va 4) + Vey Voy —A/4—qy 


b) Impulssatz. Der Impulssatz lautet 


1 
A 
Cy mig tga | gtdn. (13) 
Fur gp? erhalten wir aus (12) : 


y? —] gi ee (A/4) 4 =a poe Ing:+++> + 4(la=e (A/4) 4 pos’) 


x 


ay 4 Ve jas (A/4) n male Cy a Jing oe te “11,2 ese (14) 
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Bei der Integration vernachlassigen wir das Glied ¢,/x? In2 <:++>, das fiir groBe Re-Zahlen 


von hoherer GréBenordnung klein ist 1. Mit den Abkiirzungen B Beery. — ats ergibt 
dann die Integration m ‘ 
1 
fg? dy =1 +1 + yy? — p41 pat yy 
0 
1 
+ (06 + 7) — B15 5 — 266" + 7)" + 16? + 7)? — BS 
_4(82 + yy 442 4) — 292-2 Sig» +e +8 
(B ate y) Be 4 (B “ts y) 3 B 3 y In A p? (B2 ab y)1/2 — p° (15) 


Gleichung (11) schreibt sich mit den obigen Abkiirzungen 


& oe \ ie pepo ce 3) eB Vy 6 
1 2 [(B? + y) B] eae ay, (p2 + yi? +8 + B}—1 + In (2 48 Ref) + % *—!. (16) 

c) Druckanstieg und Wandschubspannung. Die beiden Gleichungen 
(13) und (16) lassen sich leicht approximativ lésen. Wir sind so vorgegangen, da wir aus (16) 
die zu vorgegebenen J baw. y zugehdrigen c, bzw. B mittels graphischer Methoden bestimmten; 
das gibt eine Kurve c, iiber 2. Aus (15) kann man dann unter Benutzung der obigen Ermittlung 


1 
die Kurve tg a tiber A nach der Formel (13) tga = (c, —ajay/ f gy? dy berechnen. Wir erhielten 
fiir den Druckanstieg und die Wandschubspannung 0 


a A theor A exp a Cy theor 


Ci exp 
is —0,047 —0,043 igs 0,00149 0,00137 
ie —0,092 —0,085 2° 0,00109 0,00102 
a —0,133 =) PAL oe 0,00082 0,00075 
4° —0,170 —0,140, 4° 0,00062 0,00042 


Vgl. ferner Abb.5 und Abb. 6. Die Ubereinstimmung ist als brauchbar anzusehen. Mit 
wachsendem Druckanstieg werden die Abweichungen der theoretischen von den experimentel- 
len Werten gréBer. Hierzu ist zu bemerken, dafs die von uns abgeleitete Form des Wandge- 
setzes ausdriicklich auf nicht zu starken P 
Druckanstieg beschrankt ist (vgl. Abschn. 2), 


d) Geschwindigkeitsvertei- Geschwindigke/tsverteilung seh 
lung und Verdrangungsdicke. ae 
Die experimentellen Geschwindigkeitsvertei- 
lungen werden fir mittleren Druckanstieg 
ebenfalls brauchbar wiedergegeben (Abb. 3). 
Fiir starkeren Druckanstieg wird die Wieder- 
gabe schlechter (Abb. 8). Dabei wurden die 
theoretischen Verteilungen nach (12) mit : 
den oben ermittelten theoretischen Werten Messung Nikuradse 
fiir c, und A berechnet. 0 a=7° 
e fae 


Entsprechend liegen die Dinge hinsicht- 
lich der Verdrangungsdicke. Fiir 


1. Naherung 


(Wandgeserz) 
1 


ace (1 — ¢) dy 
0 0 02 04 06 08 10 

kommt mit (12) “ Abb. 3. 

oy 4 “ Ye, —A/4 i es au Vicy —A/4 , x Cy (te AA =) 

5/2 a 3 (—A/A) ( H 2 shee ot oe (—A/4) x? ot *}° ih 


hs 


1 Aus (11) folgt fiir groBe Re \c,/x ~ €, ¢ kleine Zahl (~ GréBenordnung). Gleichung (13) ergibt 
—1/4 ~ tg o fiir groBe Re. 
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«(Bar — (Blase 


i 0,138 0,151 
aS 0,179 0,198 
3° 0,210 0,245 
4° 0,234 0,309 


Vgl. ferner Abb. 10. 


5. Der Koeffizient u*/x der Mischungswegverteilung. Um zu einer besseren Approximation 
fiir starken Druckanstieg zu gelangen, werden wir die Integration der Grenzschichtdifferential- 
gleichung in einem zweiten Schritt weiter ins Innere der Grenzschicht hinein fortsetzen. In 
diesem zweiten Schritt werden wir die bisher vernachlassigten Tragheitskrafte in erster Nahe- 
rung erfassen, was aus Griinden der GréBenordnung wiederum erfordert, den Vorgang der tur- 
bulenten Mischung nunmehr auch in zweiter Naherung zu beriicksichtigen gemaif§ dem Ansatz 


l * 
pp wen tat = an(l +a). (18) 


Es wird also das zweite Glied in der Taylorentwicklung der Mischungswegfunktion beriick- 
sichtigt2. Nikuradse® hat die Mischungswegverteilung mittels der Prandtlschen Hypothese 
t = ol? (du/dy)? aus der experimentellen Geschwindigkeitsverteilung und der (mit Hilfe des 
Impulssatzes gewonnenen) Schubspannungsverteilung bestimmt, wobei graphische Differen- 
tiation zur Anwendung kam. Dieses Verfahren ist, insbesondere in Wandnahe, mit Unsicher- 


heiten behaftet. Hinzu kommt die schon erérterte Schwierigkeit der Ermittlung der Schub- | 


spannung in Wandnahe (vgl. Abschn. 3). 


Wir haben die von Nikuradse berechneten Werte der Mischungswegverteilung nach der 


Formel 


I 
BY AEN eg ea eek ee 
7 (: a} scares ) 
mit x = 0,6 ausgewertet (Abb. 4). Die Werte in Wandniahe streuen erwartungsgemaf} sehr 
stark und wurden aus der graphischen Darstellung fortgelassen. Es ergibt sich einigermafen 
deutlich die wichtige Feststellung, daB uw*/x vom Druckgradienten der Strémung unabhangig 
ist: 


* 
a 3 EE bob 
wx 


Wenn man die Ergebnisse der vom Verfasser* durchgefiihrten Analysen anderer Strémungs- 
modelle beachtet, so kommt man zu der Annahme, dab y*/x sogar als universelle Konstante 
anzusehen ist. Die Untersuchung der von Nikuradse®> durchgefihrten Messungen turbulenter 
Rohrstrémungen, die einen weiten Bereich der Reynoldsschen Modellzahl bedecken, ergab aus 
einer von Prandtl aufgestellten Interpolationsformel fiir die Mischungswegverteilung u*/~ = —1,1, 
Die Analyse der von Schultz-Grunow® gemessenen turbulenten Strémung langs der ebenen 
Platte bei sehr groBen Werten der Reynoldsschen Modellzahl ergab durch Angleichung der 


theoretischen Geschwindigkeitsverteilung an die experimentellen Werte in Wandnahe py*/x = — 


— 1,33. 


Es mége noch bemerkt werden, daB folgende Hypothese nahe licgt: ,,Setzt man 1/6 = x f(y), 


so besitzt die Funktion f(7) universellen Charakter. Jedenfalls kann man aus Abb. 4 ent- 
nehmen, da f(7) vom Druckgradienten der Strémung unabhangig ist. 


6. Berechnung der ebenen Diffusorstrémung in zweiter Naherung. Wir gehen nun daran, 
das Wandgesetz (10) weiter ins Innere der Grenzschichtstrémung hinein fortzusetzen (vgl. 
Abschn. 5). Die ttber den Wandabstand integrierte Grenzschichtdifferentialgleichung der 


’ Die hier angegebenen Werte wurden von uns aus den von Nikuradse gemessenen Geschwindigkeits- 
verteilungen graphisch ermittelt. Die von Nikuradse angegebenen Werte sind z, T. fehlerhaft. 

> Wir verwenden hier die Bezeichnung u*, um Verwechslungen mit der Zahigkeit  auszuschlieBen. 

3 J. Nikuradse, a. a. O. 

4 W. Szablewski, a. a. O. 

5 J. Nikuradse, GesetzmaBigkeiten der turbulenten 6 i = - 
Reaeciieses g Strémung in glatten Rohren, VDI-Forschungs 

6 F. Schultz-Grunow, Luftfahrtforschung 17 (1940), S. 239. 
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ebenen turbulenten Kanalstrémungen lautet! bei Vernachlassigung der laminaren Reibung 


n 


A 1 \2 /dp\2 
an +(r] (zr) = ¢ —tga | otdy, (20) 
0 


Formale Integration ergibt 


Ve, — (A/4) n — tg af p” dy 
a ES dy}. (21) 


b/2 


ee il : 
Fir aD nehmen wir an 


l * 
5a ~ n(t +n). 


Unter Beschrankung auf kleine? 7» wollen 
wir folgende Linearisierungen zulassen: 


1 1 * 
7 ~ =e ( aa (22) 
% ni ae aa n) 


; nt oe gl 
A A 
. Ge ee [ev dy ~ V4 errs 
0 


2 
ae J gp? dy 
L—-- 5 (23) 
a ae U Abb. 4. 


Zu der letzteren Linearisierung ist zu bemerken, daf fir Druckanstieg 2 < 0 ist. Aus (20) 
folgt 


n 
as = tga | gtdn, 
0 


so daB die Linearisierung in der obigen Form sinnvoll ist’. 
Damit erhalten wir aus (21) 


ee recur tga fan ( # 
ei “1 er Aa A i) ar 


oder 


mit 


1 W. Szablewski, Ing.-Arch. 20 (1952), S. 37. 

2 Als untere Grenze der Schicht des Wandgesetzes haben wir in (8) die Bedingung 7 > 1/(% Re \cy) an- 
gegeben. Diese Zahl konvergiert fiir Re + oo gegen Null, wie man aus (11) schlieBen kann. 

3 Die beiden vorgenommenen Linearisierungen sowie die noch weiter unten zu treffende Vernach- 
lassigung eines Bestandteils des Tragheitsintegrals, die fiir kleine 7 exakte Operationen sind, haben fiir 
gréBere 7 die Bedeutung einer Mittelwertbildung, da durch sie sowohl Zahler wie Nenner des Integranden 
in (21) vergréBert werden. 


19 
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ie 2) Yq— (ian \1 * 


I ({'9* an) z "ay (25) 


” n 
( f an) 
~ a z dy fur kleine 7. 


2 9 n Vey — (A/4) 0 


Dabei haben wir im Tragheitsintegral einen Bestandteil vernachlassigt, der fir kleine 7 von 
hdherer GréB®enordnung klein ist!. Die Integrationskonstante wird so zu wahlen sein, da wir — 
den Anschlu8 an das im ersten Integrationsschritt (Abschn. 2) erhaltene Wandgesetz erlangen. 


Im folgenden gebrauchen wir wieder die Abkiirzungen B = Ye,/x, y = — (A/4)/ x? . 


a) Wandgesetz,. Integration ergibt bei passender Wahl der Integrationskonstanten 
das Wandgesetz (10) 


fy =2 (6 + yup —p) + pina BE VOTE + pl—1 +n (2a? Rep) + x A2, 


seam a) eed 


(26) 
b) Turbulente Mischung in zweiter Naherung. Man erhalt 
zo fir! Al 
pl ee ee a ee (27) 


c) Tragheitskrafte in erster Naherung. Wir erfassen fir kleine 7 die 
n 
Tragheitswirkungen in erster Naherung, indem wir in das Integral | gy? dy mit dem Wandgesetz | 


(10) bzw. (26) eingehen. Das Wandgesetz gibt ja fiir eine Wendaahe Zone (vgl. Abschn. 4) die | 
Geschwindigkeitsverteilung exakt wieder?. 
Fir f? haben wir nach (26), wenn wir noch die Abkiirzung { } = {— 1 + In (2 x? Ref) + | 

%U, Oo/v} einfiihren. 
PY PON Iig 6 VE 4 ae 2 1/2 __ 9 p2f.-1] Sg Pa? Ser ears 
+ 4[(B? +7)? — Be + 4B B+ n)® — Blinc: + Bt <- >. (28) 


Hier vernachlassigen wir wieder wie in Abschn. 4 das numerisch belanglose letzte Glied*®, das | 
fiir groBe Re von héherer GréBenordnung klein ist. Es folgt 


[ote =—4 2} — 2) tym — B+ SEK} — 2 + ye? — BF 


2 FH p2 2 pa | Ge A =o PW'e ea eke " 82 (B24 y AB | 
+ FUG + mt 81 + BC} 3-8) + G]n + 26} — 2) yn Es 


8 : , 2 (2 nes 
seo EN Utrera oa (29) 


Mit (29) haben wir nun in das Triagheitsintegral 


” n 
(tes) 

oe tj 
J nVe—alan 
einzugehen. 


1 Vgl. FuBnote 3 von S. 275. 

* Der Fehler, den wir dadurch begehen, daB wir im Integral die Funktion (26) bis zur Wand 7 = 0 zur | 
Anwendung bringen, ist vernachlassigbar klein. Das folgt aus der Bemerkung, da die untere Grenze der | 
Schicht des Wandgesetzes fiir Re > oo gegen Null konvergiert; vgl. FuBnote 2 von S, 275. 

3 Unterwirft man zum Zweck der Abschatzung f, nach (26) der Randbedingung f, = 1 fiir 7 = 1, so | 
folgt, daB fiir groBe Re B[] ~1 (~ GroBenordnung) sowie B ~ «, ¢ kleine Zahl. Der Impulssatz (13) 
ergibt dann y ~ tga. 
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Wir geben hier nur das Endergebnis der elementaren, wenn auch etwas mihsamen Rech- 
nung an: 


eo 5 
E> fs = ret Fp) +g] + — 8 
+ 2g + rt — n+ Sori) +n — Wadd 20k 
2 ee O02 ok nae ane Le. pp ey ay) — P| 
BIB ({-} )|(B 99)) ot rite yafin 4 (2+ y7 ae Ble (30)? 
d) Wir erhalten also in analytischer Weiterentwicklung des Wandgesetzes die Funktion 


ee eh, (31) 


mit f; nach (26), f, nach (27) und f; nach (30). 

Der weitere Rechengang ist wie in Abschn. 4. Wir approximieren durch die Funktion (31) 
die Geschwindigkeitsverteilung uber der ganzen Breite der Grenzschichtstrémung. Die Para- 
meter c, und A bestimmen wir wieder mittels der duBeren Randbedingung 


und des Impulssatzes [vgl. (13)] 
1 
A 
ete [ oan. (33) 


0 

Die Gleichung (32) der auBeren Randbedingung l48t sich nach (31) analytisch formulieren. 
Auf eine analytische Formulierung der Gleichung (33) des Impulssatzes haben wir in Anbe- 
tracht der Kompliziertheit verzichtet. 


Wir haben die beiden Gleichungen (32) und (33) approximativ in folgender Weise gelést?: 

1. Fir ein vorgegebenes « haben wir aus (32) fiir eine Reihe von j-Werten die zugehérigen 
c,- Werte graphisch ermittelt; das gibt fiir ein vorgegebenes « eine Kurve c, itber A bzw. c, — A/4 
iiber A (Kurve 1). 


2. Fiir das in 1) vorgegebene « haben wir dann zu gen in 1) berechneten J und c¢, (A)-Werten 
die Funktion a nach (31) berechnet und die Integrale jv dy graphisch ermittelt. Das gibt eine 
Kurve tga jv dy aber A (Kurve 2). 


3. Der Be iirepuni der beiden Kurven | und 2 bestimmt dann nach (33) das zu dem vor- 
gegebenen « zugehdrige A und ¢,. 


e)Ergebnisse:1.Druckanstieg und Wandschubspannung. 


a Atheor hese Oo C\ theor Ci exp 

i —0,044 —0,043 is 0,00134 0,00137 
rae —0,085 —0,085 2° 0,00097 0,00102 
3° —0,119 —0,120 3 0,00068 0,00075 
4° —0,149 —0,140, 4° 0,00046, 0,00042 


Vgl. hierzu Abb. 5 und Abb. 6, aus denen man auch den Gewinn der zweiten Naherung ersehen 
kann. Die Ubereinstimmung der theoretischen mit den experimentellen Werten ist im ganzen 
untersuchten Bereich als befriedigend anzusehen. 


1 Bei der Integration haben wir den Teil des Tragheitsintegrals, den der letzte Bestandteil von (29) 
ergibt, niherungsweise mit 


n 
8 B poe = yay 
£8 [mst (Bey nip” 


0 
berechnet unter Ppa Sean des numerisch belanglosen Gliedes 


n 
an lees ee 2 a dn}. 


(Be eis) a eB ryan) 


2 Die Durchfiihrung ier numerischen Rechnungen edenus ich meinem Mitarbeiter Herrn Herbert 
Moch. 
19* 
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2, Geschwindigkeitsverteilung und Verdrangungsdicke. Die 
mit den obigen theoretischen Werten von c, und A berechneten Geschwindigkeitsverteilungen 
(31) geben die experimentellen Verteilungen nun auch fiir starken Druckanstieg brauchbar 
wieder; vgl. Abb. 7. Die Verteilung fiir a = 4°, die etwas aus dem Rahmen fallt, haben wir in 
Abb. 8 besonders dargestellt und dazu noch die Verteilungsfunktion der ersten Naherung ein- 


getragen. 


A Gy 


Wandschubspannung 


Druckgrad 


+ 


e Messung Wkuradse © = Messing Mkuraase 


=== /heorie Theorle 


h< 
a 0,0010 
x 


—020 : 0 
0 fo Ae 3° ye je 
Abb. 5. Abb. 6. 


¥ a 


i Geschwindigkeitsverteiung % 
es oe 


05 : 
ena ie @. Naherung 
, a: / fa Messung Nikuradse 
eS Ze : * pet 
oN g —— /heorie 


i Theorie, é.Naherung 
Biel 
0 02 O4 G6 G8 10 0 Be OY G6 
ds 7 a 
Abb. 7. Abb. 8 


Abb. 9 enthalt die Geschwindigkeitsverteilung in den Prandilschen Koordinaten 


u 


u 0 Y 3 Reler 
oman an 5 (34) 


a 
OPS Vc ; v 


Als experimentelle Werte der Wandschubs 
: annun i di j i 
Wandgesetzes ermittelten Werte (Abschn. 3) ie Fe ee ae 
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Man kann der Abbildung entnehmen, wie mit groéBer werdendem a, also wachsendem 
Druckanstieg und gleichzeitig fallender Wandschubspannung, der Ubergang des Wandgesetzes 
in das logarithmische Gesetz sich immer weiter zur Wand hin verschiebt, d. h. aber, daS der 
Druckgradient immer weiter zur Wand hin sich bemerkbar macht gem4B (1). Inwieweit aller- 
dings bei starkem Druckanstieg ttberhaupt noch eine Schicht mit logarithmischem Geschwin- 


digkeitsgesetz existiert, kann erst bei Vorliegen des vollstandigen Wandgesetzes beurteilt 
werden. 


#|& 


500 


Geschwindigheitsverteilung 


Messung Nikuradse 


30,0 


eoecee log Leselz 
——— Wandgeserz > Theorie 
2. Naherung 


600 


70,0 


Abb. 9. 


In Abb. 7 ist noch die theoretische Verteilung fiir 1 = 0° angegeben. Die Verteilungsfunk- 
tion, die man analog wie bei der turbulenten Rohrstrémung! gewinnt, lautet fiir die turbulente 
Strémung im ebenen Kanal mit parallelen Wanden 


p =! [inn —(L +n] +! [14 nx Req) + 024. (35) 
Der Impulssatz ergibt hier c, = //4. Mittels der Randbedingung y = 1 fir 7 = 1 folgt die 
Widerstandsformel : 


lieu Bays Sure 3 pt v, 6 
pe le (Realty [—g—F tex) +09". (36) 
Die graphische Lésung der Gleichung (36) ergibt hier ¢ytheor = 90,0019, waihrend man aus dem 


Vergleich (Abb. 1) der in Wandnahe gelegenen experimentellen Werte der Geschwindigkeit 
mit dem logarithmischen Gesetz c, ,.) = 09,0018 erhalt. 


Allgemeiner erhailt man fir turbulente Strémungen in konvergenten oder divergenten 
Kanalen mit schwachem Druckabfall folgende Entwicklung, wenn man das Trag- 


1 W. Szablewski, Z. angew. Math. Mech. 31 (1951), S. 131. 
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heitsintegral mittels des logarithmischen Gesetzes berechnet?: 


te, ) 
p =! {Inn + An + By lng + Cylntn| +!|—1 + In (x Re fey) + 22% 
mit 
ee ie it 7 1 tga 
t 
B= —=F({]—2), 
= 1 tga 
CS Ge 


Zur Prifung dieser Entwicklung haben wir den Winkel a* berechnet, bei dem wir den 
Druckgradienten Null fiir eine ebene De ee zu erwarten haben. In guter Uberein- 
stimmung mit Nikurades (vgl. Abb. 5) ergibt sich a* = 0,13°. 


Entsprechend den Geschwindigkeits- 


i verteilungen liegen die Dinge beziiglich der 
05 Verdrangungsdicke, die wir aus der theo- 
Verdrangungsdicke retischen Geschwindigkeitsverteilung (31) 
graphisch ermittelt haben (Abb. 10) 
af | 5 5 
e §«Messung Wkuragse a ( i ee ( =D ine 
ee ee : 0° 0,099 0,095 
ie 0,166 0,151 
‘ 2 0,213 0,198 
ts 2.Maherung eS ein Z° 0,253 0,245 
: {Miherung 4° 0,287 0,309 
Bei der Konstruktion der approximieren- 
Ue ey, 1 | den Funktionen haben wir die 4u8ere Rand- 
bedingung dgy/dyj = 0 fiir 7 = 1 aufer Acht 
< gelassen. In Anbetracht dessenist die Wieder- 
0 0 78 30 ye 5° gabe der experimentellen Verteilungen iber- 


raschend gut, wenn auch nach der Mitte der 
Strémung zu die Approximation schlechter 
wird. Wir wollen jedoch die Méglichkeit einer Verbesserung durch Bericksichtgung der ge- 
nannten duferen Randbedingung hier nicht weiter erértern. 


Abb. 10. 


7. Zusammenfassung. Die theoretische Behandlung turbulenter Grenzschichtstrémungen 
mit mittlerem und starkem Druckanstieg erfordert die Beriicksichtigung des Einflusses des 
Druckgradienten bis zur Wand hin. Unter Zugrundelegung des Prandilschen Schubspannungs- 
ansatzes wird als Verallgemeinerung des logarithmischen Geschwindigkeitsgesetzes ein Wand- 
gesetz abgeleitet, das diesem EinfluB in Wandnahe Rechnung tragt. Die empirischen Koeffi- 
zienten des Wandgesetzes sind die beiden empirischen Koeffizienten des logarithmischen Ge- 
setzes und vom Druckgradienten der Strémung unabhangig. Das Wandgesetz kann unter Um- 
standen dazu dienen, die Wandschubspannung aus der experimentellen Geschwindigkeitsver- 
teilung in Wandniahe zu bestimmen. 


Der theoretischen Untersuchung turbulenter Strémungen in ebenen Diffusoren werden die 
bekannten Messungen Nikuradses zugrunde gelegt. Es zeigt sich, daB bereits mittels des Wand- 


1 Die (23) entsprechende Entwicklung lautet hier (A > 0) mit 


y 


/ 
A ” ile Yi leetera 
Os iets vO fed, fe (l— S171 —> : gy? dn fiir kleine 7, 
0 \ 8 Cy 2 Cy 


0 


n 


fs iY) 1 tga 
Pg eee ee gy? dy | + konst. 
0 


Cy 
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gesetzes Druckgradient, Wandschubspannung und Geschwindigkeitsverteilung fiir mittleren 
Druckanstieg in brauchbarer Naherung berechnet werden kénnen. Methodisch wird dabei so 
vorgegangen, da die Geschwindigkeitsverteilung tiber der gesamten Breite der Kanalstré- 
mung durch das Wandgesetz approximiert wird. Wandschubspannung und Druckgradient 
werden dann mittels des Impulssatzes und einer 4uBeren Randbedingung bestimmt. 

Eine brauchbare Wiedergabe der experimentellen Werte auch fiir starken Druckanstieg er- 
zielt man, wenn man als approximierende Funktion eine analytische Weiterentwicklung des 
Wandgesetzes ins Innere der Grenzschichtstrémung hinein benutzt, Diese analytische Weiter- 
entwicklung ist charakterisiert durch Beriicksichtigung der Tragheitskrafte in erster Nahe- 
rung und Beriicksichtigung des Vorganges der turbulenten Mischung in zweiter Naherung, 
d.h. Einbeziehung des zweiten Gliedes der Taylorentwicklung der Mischungswegfunktion 
in die Rechnung. Der dabei auftretende dritte empirische Koeffizient erweist sich als vom 
Druckgradienten unabhiangig und ist dariiber hinaus als universelle Konstante anzusehen. 


(Eingegangen am 24, Oktober 1953.) 


Anschrift des Verfassers: Dr. W. Szablewski, Deutsche Akademie der Wis-enschaften, 
Forschungsinstitut fiir Mathematik, Abt. Angewandte Mathematik, 
Berlin W 8, Jagerstr. 22/23. 
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Strémung in einem 90°-Knie*. 
Von D. Haase. 


1. Einleitung. Die Strémung in einem 90 °-Knie, das aus ebenen Wanden zusammengesetzt 
ist, soll zunachst als ebene Potentialstrémung einer inkompressiblen Flissigkeit betrachtet 
werden. Dabei wird die Annahme gemacht, da sich, ausgehend von der scharfen inneren Knie- 
kante (Abb. 1), eine Diskontinuitatsflache bildet, die von dem von der Flissigkeit durch- 

strémten Raum ein (meist als Totwasser bezeichnetes) Ge- 
biet abtrennt, in dem die Geschwindigkeit Null und der 
Druck konstant sind. 
Man kann zeigen, dafs derartige Diskontinuitatsflachen 
nicht stabil sind. Es liegt daher nahe, stromabwarts eine 
A Vermischung des strémenden Mediums mit dem Tot- 
ea wasser anzunehmen. Ist dort die Umlenkung fast beendet, 
so laBt sich ein Verlust an Strémungsenergie von der Art 
eines Carnotschen StoBverlustes berechnen. 
(ase Meee Der hiermit beschriebene Vorgang stellt eine weit- 
gehende Idealisierung der Strémung einer realen Flissig- 
keit dar. Man kann aber beobachten, daB es ein Gebiet 
gibt, das etwa dem ,,Totwasser® entspricht; allerdings herrscht dort eine recht lebhafte und 
ungeordnete Bewegung. Immerhin ist eine recht scharfe Trennung zwischen dem realen 
»»Totwasser“‘ und der ,,gesunden“ Strémung zu erkennen. Es kann daher interessant sein, 
festzustellen, wie weit der reale Strémungsvorgang mit dem oben beschriebenen, idealisierten 
ibereinstimmt. 

Bei einem Versuch wurde angestrebt, zunadchst einmal die Geschwindigkeits- und Druck- 
verteilung der Potentialstrémung in einem Querschnitt nicht zu weit stromaufwarts vor der 
Umlenkung experimentell zu verifizieren. Dann wurden Geschwindigkeits- und Druckmes- 
sungen stromabwarts von dieser Stelle gemacht und Vergleiche mit Rechnungen durchgefihrt. 

Bei der Strémung der realen Flissigkeit kann sich in der auBeren Ecke langere Zeit ein 
»Eckwirbel halten, dessen Zentrum sich wenig bewegt. Diese Erscheinung wird bei der 
gemachten Idealisierung nicht bericksichtigt. Jedoch kann man aus der Kenntnis des Druck- 
verlaufs an der 4uBeren Wand mit Hilfe der Grenzschichttheorie voraussagen, daf sich die 
Stroémung vor Erreichen der duBeren Ecke von der Wand ablést. Dieser Ablésungspunkt 
wurde mit Hilfe eines Naherungsverfahrens von Gruschwitz } vorausbestimmt. 

Bei der idealisierten Strémung wurde die Schwerkraft nicht beriicksichtigt. Die Versuche 
wurden mit Luft durchgefiihrt. Bei den gegebenen Abmessungen konnte von der Wirkung der 
Schwerkraft auf den Strémungsverlauf abgesehen werden. 


/ 
Totwasser 7 


ere tl Eckwirbel 


—— 


aubere Wand 
Abb. 1. 90°-Knie, Bezeichnungen, 


2. Potentialstromung im Knie, insbesondere Randstromlinien. Die Potentialstrémung in 
einem Knie, die durch die Bildung einer Diskontinuitatsflache gekennzeichnet ist, ist ein Son- 
derfall des Ausflusses aus der seitlichen Wand eines Behalters. Dieses Problem hat v. Mises 2 
in allgemeinerer Form behandelt. Er wendete dabei die Helmholtz-Kirchhoffsche Hodographen- 
methode an und ermittelte das komplexe Potential 


W 2OLGYy 


(@ = Geschwindigkeitspotential, Y= Stromfunktion) als Funktion der komplexen Ge- 
schwindigkeit 
w= Uz —vy, 


a®D O@ ,. 
ay = ay die x- bzw. y-Komponente der Geschwindigkeit sind. 


wobei v, = 
* Dissertation Technische Universitat Berlin-Charlottenb 1953 ki i : 

Dr.-Ing. R. Wille und Professor W. Kucharski. me abel ay a8 
1 KE. Gruschwitz, Ing.-Arch. 2 (1931), S. 321. 


2 P. Frank und R. v. Mises, Die Diff. tialglei ° i i 
Figg Sit Bessie sae ifferentialgleichungen und Integralgleichungen der Mechanik und 
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Fir die Strémung in einem 90°-Knie ergibt sich (Abb. 2) 
w \* 1\/w\ 
Saeiia Ge Gales | 
Af w\? : : 
Pees) 
v7 


wobei UY, v, die Geschwindigkeiten im Unendlichen sind. Fiir die Strahlkontraktion ¢ er- 
halt man mit Hilfe des Impulssatzes 


(2,1) 


l—e l-+fe 
Ta = In SH 2 ge. (2,2) 


Daraus folgt ¢ = 0,52555 . 


Es werde nun W(z) gebildet (s = x +i y), indem 
aus (2,1) w mittels der Gleichung 


dw 
w= @e (2,3) 
eliminiert wird. Lést man (2,1) nach w? auf, so wird Aled. Potcutialembmube iia Kase, 
iW n W | 
1 v? e+] e +1 1 — «2\2 
rg es ie i ed a na (a pees? 
ef % + il 
sodafi wegen (2,3) 
/ nw 
2 C Uo 
aug eee (ead ; (2,4) 


vy 2W i nW 
ee 8) 
€ é é 


Fuhrt man in (2,4) eine neue Veranderliche « ein durch die Substitution 


2 
wan 2 E*Ysinta —], a 
so wird 
ae i278 Wooese a) cos a ane (2,6) 
aes) — cos? a 


Die Integration von (2,6) ergibt, wenn man die Integrationskonstante Null setzt und die Vor- 
zeichen passend wahlt, 


1 + & 1— é 
Ob irae 2 see: ! eee ee a 26 ee | (2,7) 
at rs te 1 eee, 2€ ae 4 ‘ 
= 5 sin a 1 — Bi tga 


Das komplexe Potential W an einem Orte z erhalt man nun, indem in (2,5) und (2,7) jeweils 
der gleiche Parameter « eingesetzt wird. Das Auffinden der zu wahlenden Werte « wird erleich- 
tert, wenn man in einigen Fallen eine Konstante g durch die Gleichung 


e=tg (2,8) 


ro|-S 


einfiihrt, so dai 
ie 2k t Pe Pid eet 
f= er L+e to 


1 Fiir die geraden Wande allein kann man W(z) einfacher bestimmen, wenn man bereits in (2,1) 
vx = 0 oder vy = 0 setzt. Hier soll jedoch die allgemeine Lésung gesucht werden. 
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Damit kann man (2,5) und (2,7) auch schreiben: 


W =5In a (fre) 72): (2,5a) 
sin ~ | 

sin « tg a 

¢ | De fy Se en “tse | 2a (2,7a) 
samme sin oy __ sine tg p 1 184 ‘ 2 
sin tg 
Mit 

OO =O, +a, (2,9) 


(x,, X, reell) wird wegen 
sin x = sin a, Gof a, + i cos a, Sit a» 
schlieBlich 
W =*in d 
ee sin? p 


= (sin? «, Col? a, — cos?a, Sin? a, + 27 sin a, cosa, Sina, Gvj «,) — 1]. (2,5b) 


Durch Trennung von Realteil und Imaginarteil erhalt man daraus das Potential 


Fass re 
sin4 


Te a, Sof? a, — cosa, Sin? a, — sin’)? + 4 sin? a, cos* a, Sin? x, C0]2.x5]! (2,10) 


J 


und die Stromfunktion 
2 sin a, cos x, Gin a, Cof «, CU 


pated = 
Wires = are tg : = are tgy, (2,11) 


sin? x, Coj? a, — cos? a, Gin? a, — sin* 


wenn man fiir das Argument des Arcustangens y schreibt. Der Wert y = 0 entspricht den 
Stromlinien ¥ = 0, +c,...3 y = 0 ist etwa durch ao, = 0 zu befriedigen. Dann wird fir 
Oo, = @ die Variable z unendlich grof8. Daher sei zunachst 


a) &% =0, a, =0...qg. Aus (2,7a) folgt dann 


c= 0 
| SPAS: 1 Urs sin 0, = ae Gq Bo teat (2,12) 
|= ele : cat tas Sieg tm (%,.—> p: y—> &) 
und aus (2,10) | 
yal sin? a, i eter. y — into (2 13) 
n sin? I sin? g }° “s 
Aus (2,12) und (2,13) folgt 
_d@ da, —_% Oy | 
an ap oe ea ce | 
Eliminiert man mittels (2,14) «, aus (2,12), so ergibt sich | 
ie == (0) 5 
Pepa einen ee + 2earetg eT). | (2,15) 
b) 0, = 0 0, =... + . Hier folgt entsprechend 
ig = — ¢5 | 
ae: 1 sin » 2,16 
Poa slap tt sin a tigg Mtg 2 +a] ae 
mit 
OC aa Pp ¥ => (o,@) 
a a A a Sees 2¢€ 4 
Cram, =| Fe Oe a+ 
1 
=e be +9 or Sg e+ n| und wegen (2,2) 


es 1 +22 1 — ee 
sae rai nl =, 
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fg al he 
es TS (int : 1), (2,17) 
v o 
ty = — tg (3): (2,18) 
Me) = Ov, 
2G 
|y =~ 5 [Me tg + tae Ty et + 2earcege |. GY) 
si) 7) 1) 
oS > % = 0... 00 
2c Sin a 1 Tq a 
x= — — arc t 2+ — arc tg —= 4 4 x On = 0 3x = 0 
m f Ssing ' tgp eee 7) (2,20) 
y =0 (a> co : > — ov) 
__ € % Sin? a, 
=f = Fain sin? p +1), or 
v ow 
v1, = — = Tq a ’ (2,22) 
| x Sree Meo yg et earctge— +2 Ur tq— 
é é Dy Bars v, |” (2,23) 
ly = 
05 1 15 2 25 3 y 
“Z| 10000 | tose | 1weo | 12898 | 14501 7679 r 
70000 | 9902 | 077% | osm | geno | gseo4 | 95566 
35 3,65 4 45 § 6 oo 
“F| 19462 | 19984 | 27154 [a2aae | 2,4499 eG = 
95487 | 05429 | 05378 | 05330 | 05307 | Q5272 | 052555 
Totwasser 
i Diskontinuitatslinie 
8 | 
S 
Ss 
(ae eS ee x 
Abb. 3. Potentialsstrémung im Knie; Diskontinuitatsflache, 
d) o, =—, & =0...c (Gleichung der Diskontinuitatslinie). 
1 2 2 & 
Pass peak | plea z aretg 8% 1 2 4,], 
14 sing tg Pp tg Pp 2.24 
eh ITY PM ots a, . (2,24) 
y x © | sin Y 8 Coja, a 
mit 
1 2 
A= 0:4 = —C y= Se] team age +a 
und unter Beriicksichtigung von (2,2): 
° Med 
OX, —> 00: X—> — CO 
yoce. 


Wegen (2,11) und o, = 7/2 stellt (2,24) eine Stromlinie dar, und zwar die Diskontinuitats- 
inie (Abb. 3). Am Rande sei vermerkt, da8 aus (2,24) folgt 


ees | (2,25) 
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3. StoBverlust. Nimmt man an, daf nach erfolgter Umlenkung eine Vermischung desi 
stromenden Mediums mit dem Totwasser stattfindet und da® diese Vermischung plétzlich mit 
groBer Heftigkeit einsetzt, so 1aBt sich mit Hilfe des Impulssatzes ein StoBverlust errechnen, 
Setzt man stationare Strémung voraus und verwendet man die eindimensionale Stromfaden~ 


theorie, so ergibt der Impulssatz, angewendet auf eine Kontrollflache nach Abb. 4, wenn mani 
die Wandschubspannung zwischen den Stellen 2 und 3 vernachlassigt, | 


o F we (ws — w,) = (pp — ps) F. (3,1) 


Dabei bedeuten w;, p; die Geschwindigkeit und den Druck im Querschnitt 1; F ist der Quer-: 
schnitt des Rohres, 9 die Dichte des strémenden Mediums. | 
Mit Hilfe der Kontinuitatsgleichung und durch Ansetzen der Bernoullischen Gleichung) 
zwischen den Stellen 1 und 2 (Abb. 4) kann man aus (3,1) einen Druckverlust berechnen: 
1 1 ase | 

Pi — Ps =z Ol, — 4)? => 0 i= oe 1 -  (3,2)) 


Definiert man einen dimensionslosen Widerstandsbeiwerti 
2. Pie Pa ote | 

2° | 

so ergibt der Vergleich von (3,2) und (3,3) 


tie & a 1). (3,4)| 


Kontrollfliche i‘ 
Abb. 4. StoBverl i 6 i : 
Ee ee Mit ele= 0.52555 hatman ale Zablenwert ci 0. ie 


Nach Messungen am Hermann Foéttinger-Institut ! betragt der hier berechnete Wert ¢ etwas 
zwei Drittel eines experimentell ermittelten Faktors, der den durch die Umlenkung zusatzlichh 
bedingten Verlust an Strémungsenergie reprasentiert. Dieser Faktor wurde durch Messung 
des Druckabfalls zwischen zwei Stellen 1 und 3 (Abb. 4) bestimmt, die so gewahlt waren, daB in 


Ablosung. 10 
(berechnet) 


a 
| Stromungsrichtung 
20 18 10 05 Bares 10 15 20 25 Jae 
z ey 
Lcke 


Abb. 5. Druckverlauf langs der AuBenwinde. berechnet; — — — gemessen. 


beiden die Verteilung der Geschwindigkeit ititber den Querschnitt nahezu der der voll aus-- 


gebildeten turbulenten Strémung entsprach. Daher ist eine volle Ubereinstimmung mit dem) 
oben berechneten ¢ nicht zu erwarten. 


' A, Druck- und Geschwindigkeitsverlauf lings der Winde bei Potentialstrémung. Sind v, 
die Geschwindigkeit bzw. der Druck an einer beliebigen Stelle des Knies, v9, Po die Werte fir 


1 R. Wille und D. Haase, Allgemeine Warmetechnik 4 (953) eons 
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das dem Knie zustromende Medium in der Entfernung Unendlich, so hat man durch Anwen- 
dung der Bernoullischen Gleichung und der Kontinuitatsgleichung 


bien ] l v\? 

pal (es )]: a) 
Damit kann man bei bekanntem Geschwindigkeitsverlauf auch den Druckverlauf bei Poten- 
tialstrémung berechnen. Dieser ist fiir die Wande in Abb. 5 und 6 eingetragen. 


# 


7 | Stromungsrichtung 


20 15 


Ole 


Kone 
Abb. 6. Druckverlauf langs der Innenwiande. 


berechnet; — — — gemessen. 


5. Die Versuchsanordnung (Abb. 7). Das Knie war aus einem Rohr mit quadratischem 
|Querschnitt angefertigt; die Seitenlange des Quadrats betrug 


c = 265 mm 


(ce wird bisweilen auch mit ,,Rohrdurchmesser“ bezeichnet). Die Wande bestanden aus gut 
gehobelten Kiefernholzbrettern. Die Luft wurde von einem Radialgeblase in einen Kasten von 


44-144? 


—+| 265-265mm? i= 


G 
i 10, 7a: 18a: DruckmeBstellen in halber Cc S 
Kanalhohe — Le 
I, 0 Mebebenen,in ae Staurohre eingebracht al oa 
werden konnen 
1a 
Slog 
® 20 
cL c. cL 
ame Calne SaaS 
} At Pe 
Hi It Bi | hat 
fa “¥ 
om 
& 6a 
7a oD 
Oise 
78a a ba | ba Wa Ba tae tie wa do od c 
6 
uA A g CC ee Oe 
peri net Rete aay cee ea 
6 5c aa : aC 


Abb. 7. Schema der Versuchsanlage. 


1,44 x 1,44 x 2 m3 Inhalt gefordert. Von dort trat sie iber eine gut ausgerundete Diise in das 
Knie und weiterhin tuber eine gerade Rohrstrecke wieder in den Versuchsraum. 

In der Symmetriecebene des Knies waren an den Innen- und Au®Benwanden Druckstutzen 
aus Messing angebracht. AuBerdem konnten in die Querschnitte I, und II, die je einen Durch- 
messer von der inneren Kniekante entfernt waren, verschiebbar Staurohre eingebracht werden. 
Es wurden Prandtlsche Rohre mit einem Durchmesser von 3 mm verwendet. Die Messung der 
Driicke erfolgte mit verstellbaren Schragrohr-Manometern mit Alkoholfillung. 
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6. Druckverteilung lings der Wande (experimentell), In einer Voruntersuchung wurde mit 
Hilfe eines Staurohres der statische Druck in der Mitte des Querschnitts I (Abb.7), einen Durch; 
messer vor der Kniekante, bestimmt und mit den Driicken an den Wanden im gleichen Quer: 
schnitt verglichen. Der Druck in der Mitte des Querschnitts I hatte die gleiche GréBe wie dex 
Mittelwert der in der Symmetrie-Ebene (des Knies) gemessenen Wanddricke, die von dem 
Druck in der Querschnittsmitte um etwa + 1% baw. — 1%, abwichen (wenn man den Druc 
im Totwasser als Bezugsdruck Null annahm). 

Dies stimmt gut tiberein mit den Driicken, die bei der in den vorigen Kapiteln behandelte 
Potentialstrémung an der Stelle y/e = 2 auftreten. Es ist nach (4,1) usw. 


P/Po = 1,0094 fiir —a/c=0, Wi Caer. 
P/Po = 0,9900 fiir — xc a 1 A yle = o : 


Mit Hilfe des Druckes py in der Mitte des Querschnitts I wurde die mittlere Geschwindig 
keit v,, bestimmt mittels der Formel 


1p A 1 
=—|/2 , 6,1) 
Ye y/ Q 1 + Cpr ( if 


Dabei ist pg der Druck vor dem Eintritt in den Kanal; Cg, ist ein Diisenbeiwert, der mit 0,02 
angenommen wurde. Die Luftgeschwindigkeit v,, betrug 23 m/s, die kinematische Zihigkeit 
1,5 - 10-5 m?/s. | 

Um nun einen Vergleich der gemessenen Druckverteilung mit der der Potentialstrémun 
durchfiihren zu kénnen, werde die Geschwindigkeit v,, im Querschnitt I gleich der Geschwin 
digkeit v) der Potentialstrémung im Unendlichen gesetzt. Zu vg gehért ein Druck py derart/ 


daf bei der Potentialstrémung im Totwasser der Druck Null herrscht, namlich 


i 


1 | 
ges He at 1 ve. (6,2) 
Zur Auswertung der Messung der Druckverteilung lings der Wande wird also der Druck i 
Querschnitt I nach (6,2) und (6,1) 


OVC tory | See ery \ar Eee . 
PI -2(3 1 Um = & 1) ] ae ges (6,3 
gesetzt. 


Der Druckverlauf ist in den Bildern 5 und 6 eingetragen. Auf Abb. 5 erkennt man deutlic 
den EinfluB des Eckwirbels. Davon abgesehen ist eine gute Ubereinstimmung zwischen de 
gemessenen Druckverlauf und dem der Potentialstrémung bis zwei Rohrdurchmesser hinte 
der auBeren Ecke ersichtlich. Dort soll nun die Geschwindigkeitsverteilung langs des Durch4 
messers in der halben Kanalhéhe untersucht werden. 


4 
4 


7, Geschwindigkeitsverteilung 1 Durchmesser stromabwirts von der inneren Kniekante., 


a) Bei Potentialstrémung. Die x-Komponente der Geschwindigkeit (Abb. 2) 
betragt 

_ 0 _8@ ax, , a Aa, 
0% O00, O * Oa, Ox (7,1) 


Vx 


Zu ihrer Berechnung sei zunachst x und y als Funktion von «,, %, angegeben. Durch Zerlegung 
von (2,6a) in Realteil und Imaginarteil hat man unter Beachtung von (2,9) 


Seta] a 2sin @ cos a, Sin a, 
c x |sing 8 sin? p+ cos*a, — Coj2a, 
1 Sige Cie (7,2) 
are t 2 
as tg p 8 tg? p(cos 2a4,+ Coj 2a,) + cos 20,—Cojf 2a, ae x, i 
Woe eres | ped in (sin g + sin a, Cof a.) + cos? a, Gin? a, 
c 2n\sing (sin p—sin a, Coj «,)? + cos? a, Sin? a, | 
(7,3) 


1 [sin 2a,+ tg p(cos 2 a,+ Coj 2 «,)]?+ Sin? 2« 
1 1 1 2 2 
nares i [sin 2x, —tg p(cos 2a, + Coj 2a,)]? + Sin? 2a, are a : | 


Um nun die Wertepaare (a,, a) zu finden, die zu — x/c = 2 gehéren, beachte man, da 
nach Gleichung (7,2) und (7,3) fiir groBe a, offenbar — x/c proportional «, (und unabhangigi 
von a) und y/c proportional «, (und unabhangig von «,) werden. Nun ist bei der Diskontinui- 
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| tatslinie der Koordinate — x/c = 2 der Wert a, = 3,65 zugeordnet. Da damit Gina, und 
| Coj a, groB wegen 1 werden, sind die oben erwahnten Proportionalitaten fiir — x/c = 2 nahezu 
erfillt, so da8 man an dieser Stelle, ausgehend von der Diskontinuitatslinie, die Parameter 
) 4, X leicht ermitteln kann. 

Die Berechnung von v, an der Stelle — x/e = 2 nach (7,1) kann man nun wesentlich verein- 
fachen, wenn man beachtet, daB fiir groBe a, 


a a Oo, Oo, 
ee eas ice 


ist. Die erste Ungleichung sagt aus, daf fiir groBe a, oder nach obigen Uberlegungen fir groBe 
— x/e die Geschwindigkeit unabhangig von y, d.h. konstant iiber den durchflossenen Quer- 
schnitt ist; die zweite sagt aus, daf a, fiir groBe — x/c unabhangig von x ist; dies folgt ebenfalls 
_aus den genannten Uberlegungen. Daher wird fiir — x/c = 2 v, durch Division der Ausdriicke 
0@/dx, und Ox/da, berechnet und anschlieBend die GréBe des ersten Produktes 
| Gleichung (7,1) abgeschatzt. (Die trotz der Vereinfachungen noch umstandlichen Rechnungen, 
die in der Dissertation ausgefiihrt sind, sollen hier fortgelassen werden.) Die Geschwindig- 
keitsverteilung ist in Abb. 8 eingetragen. 


b) Experimentelle Bestimmung des Geschwindigkeitsver- 
laufs. Die experimentelle Bestimmung der Geschwindigkeitsverteilung eines Durchmessers 
hinter der Knieecke wurde in der halben Kanalhéhe vorgenommen. Die Strémungsrichtung 
an der MeSstelle wurde jeweils mit einer Fadensonde kontrolliert. Dabei wurde festgestellt, 
daB die Geschwindigkeitsrichtung auferhalb des Totwassers nur sehr wenig von der Richtung 
} der Kanalachse abwich. Die Luftgeschwindigkeit betrug wie bei der Bestimmung der Druck- 
) verteilung v,, = 23 m/s bei einer kinematischen Zahigkeit von y = 1,5- 10—° m?/s. Der ge- 
messene Geschwindigkeitsverlauf ist ebenfalls in Abb. 8 aufgetragen. 


_gemessen o4 


= a1 


igkeit 
Were Geschwindigkeit 


F S in gto RIC TON OS SO CRED TANGO UNOS WOE OS SOZ. VOT. 0 
Geschwindi 


| Abb. 8. Geschwindigkeitsverlauf im Mittelschnitt, ein Durchmesser hinter Knie — Abb. 9. Zur Berechnung des Ablésungs- 
Innenkante. berechnet; — — — gemessen. punktes A. 


8. Bestimmung des Ablésungspunktes vor der auferen Knieecke nach Gruschwitz. Zur Be- 
'rechnung der Lage des Ablésungspunktes A (Abb.9) vor der duBeren Knieecke werde die 

Strémung lings der 4uBeren Wand des Knies (y = 0 in Abb. 9) mit der Strémung langs einer 
Platte verglichen. Fiir diese gilt nach Messungen von Hansen’, daB die Grenzschicht laminar 
ist von der Vorderkante (x = 0) bis zu einer Stelle x = x,, fiir die die mit der Lauflange x und 
| der Geschwindigkeit U auBerhalb der Grenzschicht gebildete Reynoldssche Zahl den Wert 


Re,, = = 3,2: 10° (8.1 


} annimmt. Bei » = x, wird die Grenzschicht turbulent?. Es sind zwar fir Re, Werte von 108 
und hoher erreicht worden, jedoch nur bei besonders stérungsfreier Strémung, eine Voraus- 
setzung, die bei unseren Versuchen als nichterfillt angesehen wurde. Mit U = 23 m/s, 
oe — 1,5 -10-*> m?/s wird nach (8,1) 
we 0.21 mi 0,8 e% (8,2) 


| Bei unserer Anordnung wurde x, nur langs der geraden Kanalwand abgetragen; die kurze 
| Einlaufdiise blieb unberiicksichtigt. In Abb. 9 ist zur Berechnung des Ablésungspunktes ein 


1 M. Hansen, Z. angew. Math. Mech. 8 (1928), S. 185. 
_ 2 Der Buchstabe u als Index bedeutet Umschlagpunkt laminar-turbulent. 
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neues Koordinatensystem gewahlt worden, dessen Nullpunkt an der Stelle liegt, die bei unsere 
Vergleich der Vorderkante der Platte entsprechen soll. Die Koordinaten werden wieder mit 
x, y bezeichnet. Mit diesen Festlegungen liegt dann der Umschlagpunkt zwei Durchmesser 
vor der 4uBeren Knieecke. 

Wir nehmen an, daf in der Grenzschicht in y-Richtung der Druck konstant ist, und be; 
rechnen diesen aus der Geschwindigkeit U der Wandstromlinie y = 0 bei Potentialstrémung; 
Diese Geschwindigkeit U entspricht bei der Platte der Geschwindigkeit auBerhalb der Rei+ 


bungsschicht. Es ist nach (2,15) mit 
2,8 —x/c anstatt y/c, 


” me a (8,3) 
Mf 2 —X, J 
Dh eR pat eg 7 + et Mr Tg et + 2earceges |. (8,4) 
c Lf Vo U% Vo 


Nach Ziff. 6 ist an der Stelle x = x, der Druck gegeniiber dem Druck im Unendlichen um 
etwa 1% angestiegen. Nun ist x, berechnet worden unter der Annahme, daf der Druck langs 
der Wand konstant bleibt, da® also der geringe Druckanstieg die Lage des Umschlagpunktes 
wenig beeinfluBt. Dafiir ist jedoch zumindest zu fordern, daB nicht bereits vor dem errechneten 
Umschlagpunkt die laminare Grenzschicht sich von der Wand ablést. Um dies zu tiberprifen,, 
kann man nach Prandil! folgendermaBen verfahren. Es sei angenommen, da bis zu der be-| 
trachteten Stelle hin das Profil der laminaren Grenzschicht sich unter Wirkung des Druck-, 
apstieges so weit entwickelt hat, daB nur wenig stromabwarts eine Ablésung méglich ist.) 
Prandtl stellt nun die Frage, wie die Druckverteilung beschaffen sein muB, damit eine solche 
Ablésung gerade vermieden wird, d. h. welchen Druckanstieg die so weit entwickelte laminare | 
Grenzschicht zu tiitberwinden vermag. Prandtl findet an Hand des Kdrmdn- Pohlhausenschen 
Naherungsverfahrens als Kriterium, da, solange der Parameter 


U Uz 
= Sy > ll (8,5) 


(U’, U” Ableitungen nach x) ist, keine Ablésung stattfindet. Man erhalt nach (7,4) 
2 2 
ey ee | 
= 0 Vo \ 
i (8,6) 
Pe eG EGP 9G), 
Vo OF Vp 


und numerisch fiir drei Werte von we die Werte der Tabelle 
i) 


VY 0,99 0,98 0,95 
0 

x 
Fs 0,747 0,97 1,275 
oO 22,9 11,45 4,61 


Z | 
Demnach ist vor dem errechneten Umschlagspunkt von der laminaren in die turbulente Grenz-; 
schicht (x = x, = 0,8 c) keine Ablésung der laminaren Grenzschicht zu erwarten. | 


Zur Bestimmung des Ablésungspunktes der turbulenten Grenzschicht nimmt nun Grusch- 
witz an, daf} sich die Form des Geschwindigkeitsprofils der Grenzschicht durch einen einzigen 
Parameter kennzeichnen laBt, und bestimmt aus Versuchen den Wert des ,.f ormparameters“, 
fiir den Ablésung eintritt. Ist u die Geschwindigkeit in der Grenzschicht, so bezeichnet, wenn 
h eine Stelle auBerhalb der Grenzschicht bedeutet, 


=F | u(U —u)dy (8,7) 


Dia 


1 L. Prandtl in W. Durand, Aerodynamic Theory III, S. 112, Berlin 1935. 
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a 


die ,, Impulsverlustdicke. Gruschwitz definiert unter Beriicksichtigung seines Versuchsmaterials 
einen Formparameter 


u* (3) 
a Doin: 5 os (8,8) 
und findet, daB Ablésung eintritt fiir 
7 = 0,8 — Ablésung. (8,9) 


Um den Verlauf von 7 langs x zu bestimmen, stellt Gruschwitz aus Dimensionsbetrachtungen 
eine Differentialgleichung auf: 


d(U? n) 


o + A(U27) = BU?. (8,10) 


A, B sind Konstanten, die aus Versuchen ermittelt wurden. In der am weitesten vereinfachten 
Form des Verfahrens erhalt man den Verlauf der Impulsverlustdicke nach E. Truckenbrodt! aus . 


6 ) = 0, = ! 


«= xy 


{ Ut dx + cy (811) 


Ausgangspunkt zur Ableitung dieser Gleichung ist der Impulssatz der Grenzschicht. C be- 
stimmt sich aus dem Wert der Impulsverlustdicke am Umschlagspunkt x,, wenn man an- 
nimmt, daB 


o = § 


ist. Fir die langsangestrémte Platte hat man, wenn man den Druckanstieg wieder vernach- 
lassigt, 


“turbulent “laminar (8,12) 


Y Xu 
oe 0,664 |/ = (8,13) 


Damit la8t sich der Verlauf von # bestimmen und nach (8,10) dann auch der von 7. Als An- 
fangswert am Umschlagspunkt kann man 


Mu =1 (8,14) 
setzen”, Fiir die Konstanten der Gleichung (8,10) gibt Gruschwitz an 
A = 0,00894, B = 0,00461, (8,15) 
fiir die Reynoldsschen Zahlen ; 
Rey =" = 6-102 -- 5-108, (8,16) 


Wie sich spater zeigt, liegen unsere Messungen etwa in demselben Bereich der Reynoldsschen 
Zahlen. 


Die numerische Rechnung ergibt zunachst 
o, = 0,664 / 5 EOS a oes 074 an 


Weiterhin hat man an der Umschlagsstelle 


(*) — 0,9876, U, =0,9876- 23 = 22,7 m/s. 


Vo 


Damit wird nach Gleichung (8,11) 
Gea is 50 
s 


ie (1 +e[1 (22) aU 
ar) ca 


1 E. Truckenbrodt in H. Schlichting, Grenzschichttheorie, S. 430, Karlsruhe 1952. 
2 H. Schlichting, a. a. O., 5. 430. 


Nach Gleichung (8,4) ist 


dx = — 


(8,17) 


20 
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Damit wird 


ut 


v9 


+2¢ rote —(i+ pyre — 0.1]. 
Vo D 


0 


[usae = —22 (22) [oro tao + 2 Mle Be et 
WG \\ & Vo g 
Wenn man dies in (8,11) einsetzt erhalt man 9 und anschlieSend durch graphische Lésung der 
Gleichung (8,10) 7 = 0,8 (Ablésung) fir U = 16,4 m/s, U/vy = 0,713, x/c (Ablésung) = 1,966 
bzw. der Abstand von der duBeren Ecke 0,834 c. Dieser Wert ist in Abb. 5 eingetragen. 

Zur Kontrolle der Giltigkeit der Konstanten A, B, nach Gleichung (8,15), (8,16) hat man 


etwa 
an der Umschlagstelle x =0 U=22,7 m/s 0, =2,5-10-4m, 
etwa 1 Durchmesser vor auf erer Ecke x/c = 1,795 U =184m/s 9%, =15,7-10-4 m. 


Damit wird 
Re», — 379, Re», — 1925. 


Die Messungen von Gruschwitz liegen nach (7,16) im Bereich von Res = 600 — 5000. Daher 
scheint die Anwendung der Konstanten nach Gleichung (8,15) gerechtfertigt. 


9. Zusammenfassung. Die Untersuchungen zeigen, daB es unter gewissen Voraussetzungen 
moglich ist, die komplizierten Vorgange bei der Strémung in einem Knie wenigstens in der 
Umgebung der Umlenkung mit guter Naherung durch eine ebene Potentialstrémung zu be- 
schreiben, namlich wenn vor der Umlenkung die Geschwindigkeit itber den Querschnitt nahezu 
konstant ist (z. B. infolge einer kurzen Zulaufstrecke, in der sich keine dicke Grenzschicht 
bilden kann), so da Sekundarstrémungen weitgehend vermieden werden. Insbesondere 
existiert trotz der heftigen Bewegung im ,,Totwasser‘* zwischen diesem und der gesunden 
Strémung ein schmales Ubergangsgebiet, das sich in erster Naherung als Diskontinuitatsfliche 
auffassen 1aBt. Die Bestimmung des Ablésungspunktes vor der aufBeren Knie-Ecke nach 
Gruschwitz fihrt trotz der gemachten starken Vereinfachungen zu guter Ubereinstimmung mit 
dem Versuch. Der Einflu8 des Eckwirbels scheint sich auf einen engeren Bereich in der Um- 
gebung der auBeren Knie-Ecke zu beschranken. 


(Eingegangen am 4, November 1953.) 
Anschrift des Verfassers: Dr.-Ing. D. Haase, Berlin-Lichtenrade, Lichtenrader Damm 19, 


XXII. Band 1954 Schmutz: Berichtigung. 293 


Berichtigung 


zu meiner Note im Bd. XXI, S. 33 des Ingenieur-Archivs: 
+ Koeffizientenbedingungen zur Kontrolle des Dampfungsgrades bei 
Aus gleichsvorgaingen (verallgemeinerte Hurwitzbedingungen)*. 


Von O. Schmutz. 


In einer brieflichen Mitteilung machte mich Herr H. Bilharz darauf aufmerksam, da} 
meine a. a. O. in Abschnitt 10 (Anfang) und Abschnitt 11 (Ende) ausgesprochene Vermutung 
nicht zutreffe, indem er nachwies, daB fir das n-reihige Determinantenschema [Formel (21)] 
die Anzahl der notwendigen Rechenoperationen fiir den allgemeinen Fall die GréSenordnung n° 
hat, wahrend fiir das 2n-reihige Schema [Formel (35)] die Gré®enordnung n? betragt. 
Das eigentliche Ergebnis meiner Arbeit, die explizite Kontrolle des Dampfungsgrades an 
Hand von Koeffizientenbedingungen, also ohne Wurzelbestimmung, bei der die Kriterien (21) 
und (35) vom theoretischen Standpunkt aus gleichberechtigt sind, wird hiervon nicht berihrt. 


(Eingegangen am 28. Mai 1954.) 


Anschrift des Verfassers: Dr.-Ing. O. Schmutz, Karlsruhe, Kaiserstr. 12. Elektrotechnisches Institut 
der Technischen Hochschule. 
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Berichtigung 
zu meinen Aufsadatzen im Bd. XX, S.189 und im Bd. XXI, S. 396 
des Ingenieur-Archivs: ,Uber Hauptschwingungen mit 
endlichen Schwingweiten". 


Von Th. Péschl. 


Ein befreundeter Kollege machte mich vor kurzem in dankenswerter Weise darauf aufmerk- 
sam, daB die Steigung der ,,Grenzkurve“, die ich in meiner I. Mitteilung durch die Gleichungen | 
(24) baw. (25) erhalten habe, identisch sei mit der Steigung der Bahnkurve im Umkehrpunkt, | 
was im Allgemeinen nicht zu erwarten sei und speziell beim Doppelpendel nicht zutreffe. Die 
anscheinend gute Ubereinstimmung, die sich bei den von mir angestellten Versuchen gemaB | 
Abb. 4,5 und 6 mit der Rechnung ergeben hat, war der Grund fiir meine Annahme, daB die | 
Loésung der Differentialgleichung (25) auch die Grenzlagen geben miisse — ein SchluB, dessen | 
Unrichtigkeit ich der Wahrheit gemaB richtig stellen méchte. 

In allen Fallen, die auf einen ahnlichen Sachverhalt fiihren, diirfte es zweckmaBig sein, | 
zwischen einer Grenzlagenkurve und einer Grenzsteigungskurve zu| 
unterscheiden; diese fallen nur beim linearisierten Problem zusammen, sind aber sonst von- - 
einander verschieden. Die Lésung der Differentialgleichung (25) meiner I. Mitteilung gibt} 
tatsachlich nur ein Feld fiir die Grenzsteigungskurve, in das die eigentliche Grenzlagenkurve } 
eingebettet ist. Nahere Ausfiihrungen sollen in einer III. Mitteilung gegeben werden. 


(Hingegangen am 19, Juni 1954.) 
Anschrift des Verfassers: Prof. Dr. Theodor Péschl, Karlsruhe i. B., WendtstraBe 5. 
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Siebenstellize Tafeln der elementaren transzendenten 


Funktionen. Von Dr. Friedrich Liésch, Professor an der Technischen Hoch- 
schule Stuttgart. VIII, 335 Seiten 4°. 1954. Ganzleinen DM 49,80 


i inlei i Funktionen fiir + = 
Inhaltsibersicht: Einleitung. — Neunstellige Werte der elementaren transzendenten 
0 (0,0001) 0,1. — Siebenstellige Werte der elementaren transzendenten Funktionen fiir x = 0,1 (0,0005) 3,15, « = — 


3 (0,01) 10; « = 10 (0,1) 20; Werte von tg x fir x = si Werte von Urq x und Ar Cig x fir «x = 1. — Sieben- 
: td i veh 
stellige Werte elementarer Funktionen fir x = 0 (1)100. — Zwélfstellige Werte von oan fir n= 0 (1) 100. 


Siebenstellige Werte der Funktionen sin 4 %, cos = x fir « = 0 (0,001) 0,5. — Siebenstellige Werte der, Funkti- 


4 M4 ux 
Ss leekaeieee (ah Pa oe : 180 
onen 6 7 @ » Sin = x, Cof > x fir x = 0 (0,01) 2. — Siebenstellige Werte der Funktionen @ ’ 
wx 
iO Z* , Coj 3 fiir x = 0(1) 180. — Umwandlung von Bogenmaf (x) in Gradma (7). — Um- 


wandlung von GradmaB (q) in Bogenmaf (x), — Hinige oft gebrauchte Zahlenwerte. 


In den letzten zwei Jahrzehnten sind, vor allem im Rahmen des ,,Applied Mathematics Programm of the National 
Bureau of Standards‘, ausfiihrliche Tafeln fiir fast alle elementaren transzendenten Funktionen erschienen. Es fehlt 
jedoch, seit die von K, Hayashi im Jahre 1926 herausgegebenen Tafeln vergriffen sind, ein zusammenfassendes 
mehrstelliges Werk iiber diese Funktionen. Die vorliegenden Tafeln sollen zur Ausfillung dieser Liicke beitragen. 
Entsprechend der Zielsetzung schlieBen sie sich inihrer Grundanlage notwendig an die genannten Hayashischen Tateln — 
an. Das Bestreben, ein fiir viele Rechnungen ausreichendes, nicht zu umfangreiches und dabei bequem benutzbares 
Tafelwerk zu geben, fiihrte jedoch im einzelnen zu wesentlichen Abweichungen. Einerseits wurde die Stellenzahl der 
Tafelwerte fast durchweg auf 7 Dezimalen berchrankt. Andererseits wurde der Tafelschritt so weit verkleinert, daB 
bei-einem zulassigen Fehler von 1 bis 2 Einheiten der letzten Stelle im Hauptintervall fast durchweg lineare Inter- 
polation méglich ist. Zur Erleichterung der Interpolation wurden die ersten Differenzen beigegeben. 


Ein groBer Teil der Tafeln ist vollstandig neu berechnet wordén. Bei dem restlichen Teil wurde von den vorhandenen 
Quellen ausgegangen und soweit notig eine Erginzung der damit gegebenen Werte durch Interpolation vorgenommen. 
In jedem Falle wurden die auf 9 bis 10 Dezimalen ermittelten Funktionswerte durch Differenzbildung kontrolliert. 
Die Abrundung auf 7 bzw. 9 Dezimalen wurde in der iiblichen Weise vorgenommen. In Zweifelsfallem wurde die 
letzte Dezimale durch direkte Berechnung des betreffenden Wertes auf die erforderliche Stellenzahl ermittelt. 


Praktische Funktionenlehre. Von Professor Dr.-Ing. Friedrich Tolke, 
Karlsruhe. 
Erster Band: Elementare und elementare transzendente Funktionen. Zweite, stark 
erweiterte Auflage. Mit 178 Abbildungen und 50 durchgerechneten Beispielen und 
einer Ausschlagtafel. XI, 440 Seiten 4°. 1950. Ganzleinen DM 39,— 


Inhaltsiibersicht: Erster Abschnitt: Definierende Differential- und Integralgleichungen, Fundamentaleigen- 
schaften und gegenseitige Beziehungen der elementaren und elementaren transzendenten Funktionen. — Zweiter 
Abschnitt: Durch elementare und elementare transzendente Funktionen ausdrickbare Integrale. — Dritter 
Abschnitt: Funktionentafeln der elementaren Transzendenten. — Vierter Abschnitt: Anwendungen im Bereich 
der partiellen Differential- und Integralgleichungen. — Finfter Abschnitt: Summierung von Reihenentwicklungen. 
— Sechster Abschnitt: Tafel der zonalen Kugelfunktion sowie ihrer Ableitungen und Integrale. 


Z 


Grundziige der Ausgleichungsrechnung: nach der Methode der 
kleinsten Quadrate nebst Anwendung in der Geodasie. Von Dr.-Ing. Walter GroB- 
mann, o. Professor an der Technischen Hochschule Hannover. Mit 54 Abbildungen. 
VIII, 261 Seiten Gr.-8°. 1953, Ganzleinen DM 19,80 


e / 
Inhaltsiitbersicht: Uberblick tiber die Methode der kleinsten Quadrate. — I. Grundziige der Fehlerlehre: Fehler- 


| 
arten und GenauigkeitsmaBe. Das Fehlerfortpflanzungsgesetz. Berechnung des mittleren Fehlers aus scheinbaren 
Fehlern gleich genauer Beobachtungen. Berechnung des mittleren Fehlers aus Beobachtungen verschiedener Ge- 
nauigkeit. Berechnung des mittleren Fehlers aus Doppelmessungen. Das GauBsche Fehlergesetz und die Genauig- 
keitsmafe. — II. Ausgleichung von direkten Beobachtungen: Grundprinzip und Formen der Ausgleichungsaufgabe. 
Ausgleichung direkter Beobachtungen gleicher Genauigkeit (Arithmetisches Mittel). Ausgleichung direkter Beob- 
achtungen ungleicher Genauigkeit (Allgemeines arithmetisches Mittel). Beobachtungen mit Summengleichung. — 
III. Ausgleichung von vermittelnden Beobachtungen: Einfiihrung in die Methode der vermittelnden Beobachtungen. ¥ 
Aufstelien der Fehlergleichungen. Aufstellen und Auflésen der Normalgleichungen nebst Summenproben. (v v)- bs 
Proben, Schiu8probe und Rechenschema, Gewichtsreziproken und mittlere Fehler der Unbekannten. Mittlerer 
Fehler einer beobachteten GréBe. Vermittelnde Beobachtungen ungleicher Genauigkeit. Gewicht einer Funktion 
der Unbekannten, Ubersicht tiber die Ausgleichung von vermittelnden Beobachtungen. Gemeinsame Berechnung 
der Unbekannten und ihrer Gewichtsreziproken. Ausgleichung von Héhennetzen, Reduzierte Fehblergleichungen, 
Stationsausgleichungen. Trigonometrisches Einschneiden. — IV. Die Ausgleichung von bedingten Beobachtungen: 
Einfihrung in die Methode der bedingten Beobachtungen. Aufstellen der Bedingungsgleichungen. Korrelationen- 
gleichungen, Normalgleichungen und Proben. Mittlerer Fehler einer beobachteten GriBe. Bedingte Beobachtungen 
mit ungleichen Gewichten. Mittlerer Febler von Funktionen der ausgeglichenen Beobachtungen. Ubersicht uber die 
Ausgleichung von bedingten Beobachtungen. Einfache Anwendungen der Ausgleichungen nach bedingten Beob- 
achtungen, Bedingungsgleichungen in Dreiecksnetzen. Reduzierte Bedingungsgleichungen. Das Boltzsche Ent- 
wicklungsverfahren. Vermittelnde Beobachtungen mit Bedingungsgleichungen. Bedingungsgleichungen mit Unbe- 
kannten, Aquivalente Fehlergleichungen. — V. Anhang: Ausgleichung durch schrittweise Annéherung. Genidherte 
Seeger von Funktionen. Mittlere Fehler der GenauigkeitsmaSe. Schrifttum (Auswabl). — Namen- und Sach- 
verzeichnis. ; 
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